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ABSTRACT  OF  THE  DISSERTATION 


Selective  Reverberation  Cancellation  via  Adaptive  Beamforming 

by 

Dimitrios  Alexandrou 

Doctor  of  Philosophy  in  Electrical  Engineering  (  Applied  Ocean  Science  ) 
University  of  California,  San  Diego,  1985 
Professor  William  S.  Hodgkiss,  Chairman 

The  central  theme  of  this  thesis  is  the  application  of  the  adaptive  beamform¬ 
ing  concept  to  the  problem  of  selective  reverberation  cancellation,  with  an  emphasis 
on  experimental  verification.  Oceanic  reverberation  is  often  a  limiting  form  of  interfer¬ 
ence  in  echo  detection  applications.  It  may  also  be  the  "signal"  of  interest  conveying 
valuable  information  about  the  scatterers  or  physical  processes  influencing  their 
dynamical  behavior.  Because  of  stringent  limitations  imposed  by  the  oceanic  medium 
on  the  directional  characteristics  of  sonar  systems,  reverberation  is  often  composite  in 
nature  with  two  or  more  unrelated  reverberation  types  contributing  to  the  acoustic 
return.  In  that  setting  it  would  be  desirable  to  selectively  cancel  the  unwanted  rever¬ 
beration  component(s)  while  preserving  the- component  of  interest.  This  signal  process¬ 
ing  problem  whereby  both  "signal”  and  "noise"  are  constituent  components  of  the 
received  reverberation  process  is  the  focus  of  this  thesis.  The  related  problem  of 
extracting  reflector  echoes  from  a  reverberation  backround  is  also  considered 
Significant  contributions  of  this  thesis  include: 

(1)  Exposing  relevant  reverberation  properties.  It  was  shown  that  the  spatial 
correlation  characteristics  of  volume  and  boundary  reverberation  are  directly  applica¬ 
ble  in  a  reverberation  cancellation  context. 
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(2)  Critically  evaluating  a  representative  cross-section  of  adaptive  algo¬ 
rithmic  solutions  with  special  attention  to  recently  developed  deterministic  least- 
squares  structures.  It  was  determined,  based  on  theoretical  predictions  substantiated 
by  computer  simulations,  that  lattice  filters  outperform  direct  tapped-delay-line 
implementations  and  that  deterministic  least-squares  structures  are  more  suitable 
than  stochastic  approximation  solutions  for  reverberation  cancellation.  In  addition, 
the  multichannel  least-squares  lattice  was  shown  to  offer  a  significant  advantage  over 
a  cascade  of  scalar  filters. 

(3)  Demonstrating  experimentally  the  feasibility  of  selective  reverberation 
cancellation  with  real  and  simulated  data  from  a  shallow-water  experiment,  a  near- 
surface  Doppler  sonar  deployment  and  a  multibeam  bathymetric  system.  A  clear  dis¬ 
tinction  emerged  between  adaptive  noise  cancelling  with  preformed  reference  beams 
and  spatially  adaptive  beamforming  obtained  by  operating  on  array  elements  with  a 
multichannel  joint-process  filter.  It  was  shown  that  significant  signal  enhancement  as 
well  as  selective  reverberation  cancellation  can  be  achieved  through  constrained  adap- 
tive  techniques. 
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INTRODUCTION 


Oceanic  reverberation  is  the  composite  echo  from  a  large  number  of  sound 
scatterers  located  in  the  volume  and  boundaries  of  the  ocean.  In  the  context  of  signal 
detection,  reverberation  is  a  form  of  interference  which  often  presents  a  serious  limita¬ 
tion  to  the  useful  operational  range  of  an  active  sonar.  Hence,  it  would  be  desirable 
to  actively  reject  reverberation  while  preserving  the  signal  of  interest.  Although  the 
view  of  reverberation  as  a  form  of  noise  is  likely  to  persist,  a  new  outlook  has  been 
developing  casting  reverberation  as  an  information-bearing  signal.  It  possesses  infor¬ 
mation  about  the  nature  and  distribution  of  the  scatterers  and  to  the  extent  that  the 
scatterers  are  influenced  by  a  fluid  process  about  the  process  itself.  In  this  light,  the 
following  signal  processing  problem  arises:  Extract  from  a  composite  reverberation 
return  the  component  representing  the  "signal"  of  interest.  This  thesis  addresses  both 
aspects  of  the  reverberation  cancellation  problem  The  proposed  solution  involves  the 
application  of  constrained  adaptive  beamforming. 

The  following  general  approach  is  taken:  First,  the  theoretical  and  experi¬ 
mental  evidence  on  oceanic  reverberation  is  reviewed  in  order  to  expose  properties  of 
potential  relevance  in  a  reverberation  cancellation  context.  Second,  the  adaptive 
filtering  problem  is  studied  in  detail  in  order  to  delineate  the  operational  characteris¬ 
tics  of  existing  adaptive  structures  with  particular  attention  to  some  recently 
developed  algorithms.  Having  chosen  an  "  optimum”  adaptive  processor,  real  and 
simulated  reverberation  data  from  three  active  sonar  experiments  with  widely 
differing  objectives  are  processed  through  several  adaptive  beamforming  configurations 
in  order  to  establish  the  feasibility  of  the  selective  reverberation  cancellation  concept 

This  thesis  is  organized  as  follows:  In  the  first  chapter  a  review  is  given  of  the 
current  state  of  research  on  the  properties  of  reverberation.  The  physical  and  quasi- 
phenomenological  approaches  are  discussed.  Some  theoretical  predictions  of  the  point- 
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scattering  model  regarding  the  spatial  correlation  characteristics  of  reverberation  are 
shown  to  be  applicable  in  the  context  of  selective  reverberation  cancellation.  A  sonar 
system  simulation  program  (REVGEN)  is  described. 

In  the  second  chapter,  adaptive  filtering  is  set  in  the  proper  theoretical 
framework.  The  adaptive  noise  cancelling  (ANC)  concept  is  presented  as  a  special  case 
of  the  dual-channel  Wiener  filtering  solution  and  adaptive  beamforming  is  placed  in 
the  context  of  ANC.  The  fundamental  difference  between  stochastic  approximation 
solutions  and  the  "exact"  deterministic  least-squares  approach  is  exposed.  Three 
representative  joint-process  algorithms  are  tested  through  computer  simulations. 
Their  behavior  in  the  presence  of  power  disparities  in  the  reference  channel  and  during 
drastic  departures  from  stationarity  is  examined,  in  order  to  identify  the  algorithm 
best  suited  for  reverberation  cancellation.  An  additional  simulation  demonstrates  the 
advantages  of  a  true  multichannel  structure  over  a  cascade  of  scalar  filters 

In  the  third  chapter,  the  selective  reverberation  cancellation  (SRC)  concept  is 
applied  to  shallow-water  reverberation  data.  The  effect  of  adaptive  "whitening”  on 
low-Doppler  echoes  embedded  in  reverberation  is  examined.  It  is  shown  that  significant 
signal  enhancement  can  be  obtained  through  adaptive  noise  cancelling  with  fixed  (pre¬ 
formed)  reference  beams.  Certain  ANC  model  theoretical  predictions  as  well  as  the 
spatial  correlation  characteristics  of  reverberation  are  related  to  filter  performance. 
It  is  shown  that  SRC  can  be  achieved  with  constrained  adaptive  beamforming  imple- 
mented  through  a  multichannel  joint-process  filter. 

In  the  fourth  chapter,  data  from  a  near-surface  Doppler  sonar  deployment 
are  used  to  study  the  effect  of  surface  reverberation  interference  on  upper-ocean 
(volume)  velocimetry.  Based  on  the  intensity/velocity  structure  of  the  surface  return 
a  hypothesis  is  presented  involving  a  thin  bubble  layer  just  below  the  surface.  Spatial 
aliasing  effects  in  constrained  reference  beams  are  shown  to  be  deleterious  to  SRC  A 
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REVGEN  simulation  of  this  experiment  is  used  to  show  that  recovery  of  volume  velo¬ 
city  information  is  possible  if  the  spatial  aliasing  problem  is  aleviated. 

The  fifth  and  final  chapter  centers  around  a  REVGEN  simulation  of  the  Sea 
Beam  bathymetric  system.  It  is  shown  that  the  sidelobe  interference  caused  by  the 
near-specular  bottom  return  can  be  succesfully  removed  through  ANC,  pending  a 
minor  modification  of  the  Sea  Beam  data  acquisition  system. 


CHAPTER  I 


REVERBERATION  IN  THE  OCEAN 

In  the  active  mode  of  sonar,  the  regime  of  interference  which  usually  limits 
system  performance  is  the  sound  field  created  by  a  large  number  of  unwanted  echo 
returns.  The  ocean  abounds  with  objects  that  can  intercept  and  reradiate  acoustic 
energy.  Suspended  sediment,  organic  detritus,  air  bubbles,  plankton,  fish  and  minute 
discontinuities  in  the  thermal  structure,  are  all  capable  of  redirecting  sound.  Irregular¬ 
ities  of  the  sea  surface  and  the  ocean  floor  are  also  significant  contributors  to  this 
reradiation  of  sound  known  as  scattering  The  composite  echo  from  all  scatterers  is 
known  as  reverberation. 

When  the  objective  is  signal  detection,  reverberation  is  clearly  a  form  of 
noise.  This  consideration  is  largely  responsible  for  the  manner  in  which  reverberation 
has  been  treated  in  the  scientific  literature.  The  overiding  concern  has  been  with  the 
prediction  of  the  reverberation  intensity  level.  This  has  led  to  a  single-parameter 
description  of  reverberation,  the  backscattering  coefficient  which,  combined  with 
geometrical  and  beam-pattern  characteristics  gives  rise  to  the  reverberation  level  (RL) 
of  the  well-known  sonar  equations  [Urick,  1975] 

Although  this  view  of  reverberation  as  a  form  of  noise  will  always  be  valid,  a 
new  outlook  has  been  developing  casting  reverberation  as  an  information-bearing  sig¬ 
nal.  It  possesses  information  about  the  nature  and  distribution  of  the  scatterers  and 
to  the  extent  that  the  scatterers  are  influenced  by  a  fluid  process,  about  the  process 
itself.  The  information,  if  properly  extracted,  may  be  used  to  quantify  fishery  species 
of  commercial  interest  [Cushing,  1973;  Holliday,  1974],  to  examine  planktonic  commun¬ 
ities  [Greenblatt,  1980],  to  monitor  pollution  in  industrial  dumping  sites  [Orr  and  Hess, 
1978],  to  identify  ocean  bottom  types  [de  Moustier,  1985b]  or  for  remote  sensing  of 
oceanic  fluid  processes  [Proni  et.  al.,  1978;  Fisher  and  Squire,  1975;  Pinkel,  1981], 
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As  reverberation  is  being  reevaluated  in  terms  of  its  information  content,  so 
must  the  signal  processing  techniques  used  vis-a-vis  reverberation  be  reconsidered.  For 
instance,  rather  than  indiscriminate  reverberation  suppression,  the  situation  may  call 
for  the  extraction,  from  a  composite  reverberation  return,  of  the  component  created 
by  the  class  of  scatterers  associated  with  the  process  of  interest.  In  this  scenario, 
signal  and  noise  are  both  constituent  components  of  the  received  reverberation 
process.  This  problem  has  not  been  previously  addressed  and  is  the  focus  of  this  thesis. 

In  the  first  chapter,  we  have  attempted  a  brief  review  of  the  voluminous 
bibliography  on  the  physical  properties  of  oceanic  reverberation  and  the  mathematical 
methods  used  to  describe  it.  The  relative  advantages  of  the  physical  versus  the  guasi- 
phcnomcnological  approach  are  discussed  and  a  simulation  program  implementing  the 
latter  is  described.  A  direct  connection  is  established  between  certain  model  predic¬ 
tions  concerning  the  spatial  correlation  characteristics  of  reverberation  and  the  adap¬ 
tive  noise  cancelling  concept. 

1.1  Reverberation  types  /  Scatterer  identification 

Oceanic  reverberation  is  usually  classified  as  volume,  surface  and  bottom 
reverberation.  The  scatterers  responsible  for  volume  reverberation  are  mostly  biologi¬ 
cal  in  nature  [Clay  and  Medwin,  1977],  Inorganic  particles  are  insignificant  contribu¬ 
tors  and  reflections  from  sound  velocity  microstructure  are  effectively  masked  by  bio¬ 
logical  scattering  [Kaye,  1978],  Zooplankton  such  as  copepods  are  the  dominant  source 
of  volume  scattering  in  the  near-surface  region.  In  deeper  water,  biological  scatterers 
are  often  distributed  within  diffuse  Deep  Scattering  Layers  (DSL)  consisting  of  sipho- 
nophores,  copepods,  pteropods,  euphausiids  and  fish  The  diurnal  migrating  cycle  of 
the  DSL  and  its  frequency-selective  backscattering  properties  have  been  studied  inten¬ 
sively  revealing  a  complex,  multilayered  structure  [Clarke,  1970;  Kampa,  1970;  Isaaks, 
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1974;  Hersey,  1962;  Gibbs,  1970]. 

Surface  reverberation  is  generated  by  the  entire  spectrum  of  the  rough 
air/sea  interface  and  is  a  function  of  wind  speed  and  the  transmitted  frequency. 
[Urick  and  Hoover,  1956,  Chapman  and  Harris,  1962;  Garrison  et.  al  ,  I960].  Specular 
reflections  from  normally  inclined  wave  facets  and  scattering  from  an  isotropic  layer 
of  bubbles  in  the  near  surface  have  also  been  suggested  as  reverberation  sources  at 
high  and  low  grazing  angles,  respectively  [Medwin,  1966;  Clay  and  Medwin  1964 
Shadowing  effects  by  large-scale  surface  waves  may  be  significant  for  low  grazing 
angles  [Gardner,  1973], 

Bottom  reverberation  is  an  extremely  complex  phenomenon  owing  to  the 
diversity  of  ocean  floor  types,  lateral  inhomogeneity  and  potential  contribution  of  sub¬ 
bottom  layers.  The  complexity  is  reflected  by  the  disparity  among  the  reported  meas¬ 
urements  of  bottom  scattering  strength.  It  appears  that  both  particle  size  and  bot¬ 
tom  relief  are  important  factors  [Me  Kinney  and  Anderson,  1964;  Buckley  and  Urick 
1968] 

1.2  The  theory  of  scattering  from  rough  surfaces 

The  theoretical  treatment  of  reverberation  poses  a  formidable  analytical 
problem.  In  principle,  the  wave  equation  could  be  solved  if  the  boundary  conditions, 
describing  the  regions  where  density  and  sound  velocity  change  discontinuously  can  be 
derived.  In  the  real  ocean  this  is  an  all  but  impossible  task.  In  fact,  even  for  the  sim¬ 
plest  of  surfaces,  the  solution  is  far  from  trivial.  The  theoretical  solutions  to  the  prob¬ 
lem  of  scattering  from  rough  surfaces  can  be  divided  into  two  groups  according  to 
whether  deterministic  or  stochastic  surfaces  are  considered.  The  scattering  of  sound 
by  a  sinusoidal  surface  at  normal  incidence  was  first  treated  by  Lord  Rayleigh  [1945], 
His  view  of  the  sinusoidally  corrugated  surface  as  a  reflection  grating  led  to  an 
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approximate  solution  consisting  of  a  discrete  set  of  plane  waves  travelling  away  from 
the  surface  at  angles  corresponding  to  the  "orders"  of  the  grating.  An  exact  solution 
for  the  sinusoidal  surface  was  obtained  by  Uretsky  [1965]  which  was  substantiated  by 
laboratory  experiments  on  sinusoidal  cork  surfaces  [Barnard  et.al.,  1966).  Other  solu¬ 
tions  for  periodic  profiles  of  a  specific  type  have  been  obtained  [e.g  Twersky,  1951]  In 
the  absence  of  an  exact  solution,  the  Helmholtz  formula  with  the  Kirchoff  approxima¬ 
tion  [Baker  and  Copson,  1939]  has  typically  been  used.  This  approximation  is  valid  for 
"locally  flat"  surfaces  with  irregularities  whose  radius  of  curvature  is  large  compared 
to  the  wavelength.  A  good  review  paper  on  solutions  involving  deterministic  surfaces 
is  by  Lysanov  [1958]. 

A  rough  surface  is  more  appropriately  described  as  a  stochastic  process. 
Most  methods  used  to  treat  random  surfaces  are  generalizations  of  the  methods 
developed  for  deterministic  surfaces  and  are  subject  to  the  same  restrictive  assump¬ 
tions.  The  randomized  Rayleigh  approach  whereby  the  surface  is  described  in  terms  of 
random  Fourier  coefficients  was  developed  by  Rice  [1951].  Other  methods  depict  a 
rough  surface  as  a  planar  array  of  objects  the  size,  shape,  position  and  separation  of 
which  are  random  variables  or,  alternatively,  through  the  probability  distribution  of 
the  slopes  of  its  plane  facets  [Cox  and  Munk,  1954].  The  Kirchoff  method  may  be  used 
to  relate  the  statistics  of  the  surface  to  the  statistics  of  the  scattered  field  [Eckart, 
1953].  By  assuming  a  Gaussian  distribution  for  the  surface  scatterers,  it  is  possible  to 
calculate  the  mean  field  and  mean  power  scattered  in  an  arbitrary  direction.  This,  in 
turn,  makes  it  theoretically  possible  to  solve  the  inverse  problem,  that  is  to  determine 
the  autocorrelation  function  of  the  random  surface  from  measurements  of  the  scat¬ 
tered  intensity  distribution  [Horton,  Mitchel,  Barnard,  1967],  Eckart’s  high-  and  low- 
frequency  solutions  were  extended  to  include  intermediate  frequencies  by  Proud,  Beyer 
and  Tamarkin  [i960].  The  question  of  the  appropriate  form  for  the  surface  covariance 
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function  has  been  examined  by  Horton  and  Muir  [1967],  A  model  for  a  composite 
rough  surface  consisting  of  a  superposition  of  independent,  stationary  random 
processes  each  with  its  own  distribution  and  covariance  function  (e.g.  "swell",  "sea", 
"ripple")  was  suggested  by  Beckman  [1965]  and  the  effect  of  statistical  dependence 
between  the  component  roughness  processes  was  studied  by  Hayre  and  Kauffman 
[1965], 

Remarks 

This  "physical"  approach,  despite  longstanding  efforts  and  the  theoretical 
elegance  of  many  of  the  solutions,  provides  at  best  an  incomplete  description  of 
scattering  from  the  ocean  surface  (bottom),  because  of  the  restricting  assumptions 
that  had  to  be  used  in  the  face  of  prohibitive  analytical  difficulties.  Only  mono¬ 
chromatic  incident  waves  and  a  single  receiver  have  been  considered.  The  ubiquitous 
locally  flat  assumption  limits  the  range  of  transmitted  frequencies  and  surface 
roughness  scales  which  are  amenable  to  this  treatment.  In  addition,  Gaussian  statis¬ 
tics  with  convenient  correlation  functions  must  always  be  used  for  analytical 
expediency.  Thus,  the  direct  connection  of  these  solutions  to  the  physical  characteris¬ 
tics  of  the  ocean  boundaries  is  rather  tenuous  in  general.  Still,  useful  information  such 
as  roughness  and  correlation  area  may  be  derived  through  this  approach  if  these 
assumptions  are  satisfied  to  a  sufficient  extent.  No  "physical"  solution  exists  for 
volume  scattering,  however. 

1.3  The  point-scattering  reverberation  model 

An  alternative  to  the  physical  approach  is  to  represent  reverberation  as  a 
random  process  constructed  by  a  linear  superposition  of  the  individual  echoes  emanat¬ 
ing  from  a  large  number  of  point  reflectors  located  independently  in  a  homogeneous 
medium.  The  physical  equivalent  of  the  independence  assumption  is  the  first  Born 
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approximation  which  states  that  secondary  scattering  is  negligibly  small  This 
approach  has  been  developed  mainly  by  Faure  [1964],  Ol’shevskii  [1967]  and  Middleton 
[1967  a,b,  1972  a,b]  and  is  known  as  the  point-scattering  or  quasiphenomenological 
model  of  reverberation  Specifically,  the  backscattered  signal  is  represented  by 

i-4 

The  returning  echo  from  each  scatterer,  indexed  by  «,  is  delayed  by  f,- ,  the  twoway 
travel  time  between  the  source  and  the  ith  scatterer.  The  a,  and  o,  are  random  vari¬ 
ables  representing  distributions  of  acoustic  cross-section  and  Doppler  factor  of  the 
scatterers,  respectively.  For  a  monostatic  sonar, 

G{(t)-tBHn)m 

where  g  is  a  system  gain  factor,  B?  is  the  twoway  beam  pattern  sensitivity  in  the 
direction  of  the  ith  scatterer  and  f(t)  represents  the  twoway  propagation  loss.  «(f)  is 
the  transmitted  waveform 

Assuming  that  for  a  sufficiently  large  scattering  region  the  mean  scatterer 
density  is  constant,  the  number  n  of  scatterers  contributing  to  the  reverberation 
return  at  any  given  time  may  be  described  as  a  Poisson-distributed  random  variable. 
Based  on  Poisson  statistics,  a  characteristic  function  may  be  derived  for  the  rever¬ 
beration  process  and  the  various  moments  may  be  obtained  by  taking  appropriate 
derivatives.  This  approach  makes  it  possible  to  calculate  the  second  order  statistics  of 
reverberation  for  general  experimental  geometries,  beam  patterns  and  transmitted  sig¬ 
nals. 

Remarks 

Owing  to  its  general  construct,  the  point-scattering  model  allows  the  calcula¬ 
tion  of  some  fundamental  statistical  measures  of  reverberation  for  realistic  system 
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parameters.  However,  no  direct  connection  is  established  by  the  model  between  these 
measures  and  the  physical  characteristics  of  the  scattering  region.  Any  such  informa¬ 
tion  must  be  inserted  in  the  model  through  a  dynamical  impulse  response  function 
representing  the  interaction  of  the  incident  field  and  a  typical  scatterer.  Middleton 
[1967a]  suggests  that  such  knowledge  must  be  derived  from  experiment  but  he  stops 
short  of  suggesting  appropriate  experimental  procedures. 

The  simplest  approach  is  to  assume  that  the  point  scatterers  are  perfect 
point  reflectors  radiating  in  all  directions  and  reproducing  exactly  the  incident 
waveform.  Several  investigators  [e.g.  Plemons  et.  al.,  1972]  have  shown  that  even  this 
simplest  of  models  predicts  temporal  and  spatial  correlation  values  which  are  in  sub¬ 
stantial  agreement  with  experimental  estimates,  with  some  interesting  deviations. 
Such  deviations  may  be  caused,  for  example,  by  departures  from  isotropic  scattering 
due  to  "patchy"  scatterer  distributions  or  they  may  be  a  reflection  of  the 
"microroughness"  characteristics  of  a  particular  scattering  region.  From  an  extensive 
data  set,  collected  in  a  variety  of  environmental  conditions,  the  experimentally 
derived  reverberation  parameters  may  be  categorized  into  patterns  representing  dis¬ 
tinct  classes  of  scatterers.  Once  established  with  the  assistance  of  independent  physi¬ 
cal  information  such  as  bottom  photographs,  wind  and  wave-height  measurements  or 
plankton  sampling,  these  patterns  could  be  used  to  classify  scatterer  distributions 
based  on  reverberation  information  only.  This  is  the  potential  solution  to  the  "inverse 
problem"  offered  by  the  point-scattering  model,  in  this  case  an  empirical  one  Modest 
steps  have  already  been  taken  in  this  direction  [Plemons  Shooter  and  Middleton, 
1972, a, b;  Wilson,  1982;  Wilson  and  Frazer,  1983;  Jackson  and  Moravan,  1983.]  The 
results  tend  to  support  the  model  predictions  on  the  temporal  and  spatial  covariance 
properties  of  reverberation. 
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No  large-scale  effort  to  collect  long  records  of  oceanic  reverberation  data  has 
been  reported  thus  far.  Such  a  project  would  require  a  stable  ocean-going  platform,  a 
multichannel  data  acquisition  system  and  a  powerful  computing  facility  to  accomo¬ 
date  the  data  validation  and  parameter  estimation  procedures  for  the  collected  rever¬ 
beration  ensembles. 

1.4  Spatial  correlation  of  reverberation 

Spatial  correlation  expresses  mathematically  the  fact  that  returns  arriving 
from  widely  separated  angular  directions  tend  to  add  out  of  phase  at  spatially 
separated  receivers  The  point-scattering  model  of  reverberation  allows  the  calcula¬ 
tion  of  the  temporal  and  spatial  correlation  properties  of  reverberation  independently 
of  "fine  structure"  morphological  information.  Although  deviations  are  to  be  expected 
for  real  reverberation  data,  the  model  does  predict  accurately  the  basic  structure  of 
these  statistical  measures.  Based  on  the  predicted  structure,  good  decisions  can  be 
reached  on  experimental  configurations  conducive  to  particular  signal  processing  goals. 
In  this  context,  the  point  scattering  model  can  be  of  great  value.  We  will  show  here 
that  knowledge  of  the  basic  spatial  correlation  structure  of  the  backscattered  signal 
as  predicted  by  the  simplest  (point  reflector)  version  of  the  model,  can  be  used  directly 
in  the  context  of  selective  reverberation  cancellation. 

The  following  results  on  the  spatial  correlation  of  volume  and  boundary 
reverberation  are  based  the  assumptions  that  the  elementary  scattered  signals  repro¬ 
duce  the  shape  of  the  transmitted  pulses  and  that  the  reverberation  process  has  been 
appropriately  stationarized  [Ol’shevskii,  1967].  In  addition,  narrowband  (quasi¬ 
harmonic)  transmit  signals  and  uniform  scatterer  distributions  are  assumed. 

The  geometry  for  the  volume  backscattered  return  is  depicted  in  Figure  1.1a. 
The  two  receivers  /?,  and  R2  are  a  distance  »  apart  and  vertically  aligned.  Volume 
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reverberation  arrives  from  a  range  of  elevation  and  azimuthal  angles  determined  by 
the  directional  characteristics  of  the  transmitting  and  receiving  arrays.  For  omni¬ 
directional  transmission  and  reception,  the  spatial  correlation  coefficient  of  volume 
reverberation  is 

JM«>  -  ^  <2») 

where 

Thus,  the  spatial  correlation  coefficient  of  volume  reverberation  is  a  decaying  function 
and  will  tend  to  zero  as  the  transmitted  frequency  and/or  receiver  separation 
increase. 

The  geometry  for  the  boundary  return  is  described  in  Figure  1.1b.  The  two 
receivers  lie  in  the  x z  plane  and  the  line  connecting  them  forms  an  angle  4>t  with  the 
x—  axis.  In  this  case,  for  omnidirectional  transmission  and  reception  and  for  ranges 
much  larger  than  the  receiver  separation,  the  spatial  correlation  coefficient  is  given  by 

Rb(«)  =  «W*«cas^#)  (2b) 

where  J,  is  the  modified  Bessel  function  of  the  first  kind.  For  <j>c  - 0  (horizontal 

separation)  the  result  is  again  a  decaying,  oscillating  function.  However,  for  4>„-  — 

(vertical  separation)  we  have 

Therefore,  no  loss  of  spatial  correlation  is  suffered  for  vertical  separation  of  the 
receivers.  This  can  be  explained  intuitively,  in  terms  of  the  lack  of  vertical  extent  of 
the  insonified  surface  patch.  The  above  relations,  are  similar  to  some  earlier  results 
on  the  correlation  properties  of  spatially  distributed  ambient  noise  fields  [Cron  and 
Sherman,  1962;  Jacobson,  1962]. 
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The  effect  of  directional  transmission  can  be  readily  included.  Letting  MM) 
be  the  beam  pattern  of  the  transmitting  array,  we  have  for  volume  reverberation 


JT 
2  2* 


Ry(e)  ~  — o —  /  f  t}(ff,4>)cos(kssin<f>)dffcos^d^> 
A**  ,0 


(3a) 


where, 


2  2* 


An,ff  =  J  Jb}{6,<i>)decos4>d<t> 
and  for  boundary  reverberation 

**(•)  ' 


j_  0 
2 


2* 


— 7 —  f  b$(0)cos(kseo4  <j>0sin0)d0 


(3b) 


where, 


2* 


M*-J*H»)**- 

o 

Equation  (3b)  indicates  that  the  surface  return  originating  at  a  thin  surface  layer,  will 

continue  to  be  perfectly  coherent  between  vertically  separated  receivers  (<f>0  =  — ). 

2 

However,  according  to  equation  (3a),  the  volume  return  will  become  inreasingly  corre¬ 
lated  with  increasing  vertical  directivity  of  the  transmitting  array  thus  reducing  the 
correlation  disparity  between  the  two  reverberation  types. 

Experimental  evidence  exists  supporting  these  model  predictions.  Using 
explosive  reverberation  data,  Urick  and  Lund  [1964]  have  shown  varying  degrees  of 
correlation  between  vertically  separated  receivers  for  different  reverberation  types, 
with  boundary  returns  displaying  decidedly  higher  values  than  the  volume  component. 
The  same  investigators  [Urick  and  Lund,  1970a, b]  have  shown  that  the  spatial  correla¬ 
tion  of  composite  reverberation  in  a  shallow  water  environment  is  substantially  higher 
in  the  vertical  than  in  the  horizontal.  This  was  confirmed  in  a  more  recent  experiment 
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involving  a  pulsed,  high  frequency  sonar  and  directional  receivers  [Wilson  and  Frazer, 
1983). 

1.4.1  Spatial  correlation  and  selective  reverberation  cancellation 

Provided  that  the  vertical  directivities  of  the  transmitting  and  receiving 
arrays  are  sufficiently  low,  the  predicted  disparities  in  vertical  correlation  for  volume 
and  boundary  reverberation  may  be  exploited  in  order  to  achieve  selective  reverbera¬ 
tion  rejection.  It  can  be  shown  [Jackson  and  Moravan,  1983]  that  the  vertical  correla¬ 
tion  of  volume  reverberation  will  be  zero  for  receiver  separations  larger  than  the  verti¬ 
cal  dimension  of  the  transmitting  array.  If  directional  receivers  are  used  as  well,  the 
"correlation  distance"  will  extend  to  the  sum  of  the  dimensions  of  the  transmitting 
and  receiving  array.  On  the  other  hand,  surface  (bottom)  reverberation  will  be  highly 
correlated  between  the  same  vertically  separated  receivers,  at  least  for  ranges  large 
relative  to  the  distance  of  receiver  separation.  This  forms  a  natural  premise  for  the 
application  of  the  Adaptive  Noise  Cancelling  principle  (ANC)  whereby  the  correlated 
component  between  two  signals  is  cancelled.  In  this  scenario,  boundary  reverberation 
will  be  cancelled,  while  volume  reverberation  (the  uncorrelated  component)  will  be 
preserved  at  the  filter  output.  Therefore,  a  vertical  array  is  needed  for  boundary 
reverberation  rejection. 

The  ANC  concept  will  be  described  in  detail  in  Chapter  2  and  this  idea  will 
be  tested  with  real  and  simulated  reverberation  data  in  Chapter  III. 

Remarks 

i.  The  ocean  surface  will  in  general  have  some  vertical  extent  due  to  the  pres¬ 

ence  of  surface  waves.  This  will  lead  to  some  loss  of  vertical  correlation  for 
the  boundary  return  and  will  adversely  affect  the  suggested  selective  cancel¬ 
lation  procedure. 
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ii.  If  indiscriminate  reverberation  suppression  in  the  presence  of  a  coherent  echo 

return  is  desired,  a  horizontal  conventional  array  will  be  much  more  effective 
than  a  vertical  one. 

1.5  Reverberation  simulation 

The  point-scattering  model  of  reverberation  has  the  added  advantage  that  it 
can  be  readily  implemented  in  software.  Several  reverberation  simulation  software 
systems  are  available  [Princehouse,  1977;  Chamberlain  and  Galli,  1983].  REVGEN, 
(REVerberation  GENerator)  [Goddard,  1985]  is  a  direct  implementation  of  the  point- 
scattering  model;  returns  from  a  large  number  of  discrete  scatterers,  distributed  ran¬ 
domly  throughout  the  volume  and  the  boundaries,  are  summed  coherently  at  each 
receiver  to  abtain  a  reverberation  time-series.  The  REVGEN  output  is  a  dual  digital 
data  stream  representing  the  I  (in  phase)  and  Q  (quadrature)  components  of  the 
complex-basebanded  reverberation  signal.  A  backscattering  coefficient  for  each  rever¬ 
beration  type  and  the  densities  of  their  Poisson  distributions  are  specified  by  the  user. 
Scattering  layers,  random  scatterer  motion,  platform  trajectories,  attenuation  and 
reflection  losses,  arbitrary  (multiple)  transmitting  and  receiving  beam  patterns  and 
several  transmitted  signal  types  may  be  specified  through  appropriate  REVGEN 
parameters  to  create  realistic  experimental  settings.  In  such  a  controlled  environment, 
a  sonar  system  concept  can  be  tested  prior  to  actual  deployment  in  the  ocean  thus 
avoiding  the  cost  and  operational  difficulties  associated  with  oceanic  experiments. 
REVGEN  has  undergone  an  extensive  validation  process  and  has  been  found  to  be 
consistent  with  the  theoretical  model  predictions  regarding  spectra  and  correlation 
functions.  A  favorable  independent  assessment  was  made  by  Hansen  [1984]  and  we 
have  succesfully  repeated  the  original  set  of  validation  tests  during  the  installation  of 
REVGEN  on  the  MPLVAX  computer.  The  current  version  of  REVGEN  has  been 
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made  available  to  us  by  the  Applied  Physics  Laboratory  (APL)  of  the  University  of 
Washington. 

The  REVGEN  software  package  will  be  used  to  simulate  three  oceanic  exper¬ 
iments  with  widely  differing  objectives.  Comparisons  with  results  from  available  real 
data  will  be  made  in  order  to  determine  the  extent  to  which  the  performance  of  the 
adaptive  systems  in  question  can  be  accurately  assessed  through  REVGEN  simulated 
data.  Once  a  reliable  connection  is  established,  new  experimental  procedures  can  be 
suggested  on  the  basis  of  simulated  results. 

It  should  be  noted  that  such  simulations  in  no  way  eliminate  the  need  for 
further  testing  through  real  world  experiments  as  the  currently  used  model  is  only 
weakly  related  to  the  actual  physical  scattering  process.  Other  physical  phenomena 
such  as  sound  velocity  fluctuations  and  surface  waves  are  not  included  in  the  model. 
Also  not  considered  are,  necessarily,  any  other  unknown  factors  in  the  environment  of 
the  actual  experiment.  Therefore,  the  simulated  results  are  presented  only  as  a  useful 
intermediate  step  in  the  process  of  experimental  validation  of  the  adaptive  techniques. 
However,  one  may  conceive  a  process  of  continual  model  update  whereby  new  rever¬ 
beration  measurements  (e.g.  of  second  order  statistics)  are  incorporated  into  the  model 
to  eventually  produce  a  more  reliable  synthetic"  oceanic  environment. 

1.6  Chapter  summary 

We  have  presented  a  brief  survey  of  the  current  state  of  research  on  the  pro¬ 
perties  of  oceanic  reverberation.  The  fundamental  differences  between  the  physical 
and  quasiphenomenological  methods  were  discussed.  A  direct  connection  was  esta¬ 
blished  between  the  vertical  correlation  characteristics  of  volume  and  boundary  rever¬ 
beration  and  the  adaptive  noise  cancelling  concept.  A  software  implementation  of  the 
point-scattering  model  was  described. 
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Figure  1.1  Spatial  correlation  geometry  (a)  volume  reverberation,  (b)  boundary  reverberation 
(  <f>0  =  —  for  vertical  separation.) 


CHAPTER  D 


ON  OPTIMUM  ADAPTIVE  FILTERING 

As  their  name  implies,  adaptive  filters  are  capable  of  responding  to  changing 
conditions  through  a  rudimentary  "learning"  process.  The  filter  parameters  adjust 
themselves  under  the  guidance  of  an  appropriate  cost  function  and  track  the  evolving 
characteristics  of  the  input  signal(s).  These  structures  grew  out  of  the  demand  for  sys¬ 
tems  capable  of  operating  in  uncertain,  time-varying  environments  and  were  made 
possible  by  the  increasing  availability  of  computational  power.  The  general,  underly¬ 
ing  problem  is  the  one  of  designing  realizable  approximations  to  optimal  Wiener  filters 
when  the  ensemble  of  a  random  process  is  not  available.  One  is  then  faced  with  the 
problem  of  extracting  the  desired  information  from  a  single  sample  function  of  the 
process.  The  problem  may  be  approached  from  a  statistical  viewpoint  by  invoking  sta¬ 
tionary  and  ergodicity  and  performing  time-averaging,  or  through  a  deterministic 
least-squares  formulation.  Although  the  two  approaches  are  fundamentally  different, 
they  lead  to  basically  equivalent  results  when  wide  sense  stationarity  prevails.  Both 
solutions  can  be  implemented  either  in  a  block  -  processing  or  a  time  recursive  mode 
and  both  give  rise  to  efficient  lattice  structures.  Nonstationary  processes  may  be 
treated  as  locally  stationary  by  restricting  the  time  interval  over  which  optimization 
is  performed.  This  is,  essentially,  the  function  performed  by  an  adaptive  filter  of  the 
type  considered  in  this  thesis.  The  Wiener  problem,  or  its  deterministic  counterpart,  is 
solved  over  the  restricted  optimization  interval,  controlled  by  an  adaptation 
coefficient  which  provides  for  a  "fading"  of  past  information  in  favor  of  recent  values. 

Rather  than  evolving  from  this  general  mathematical  premise,  the  first  adap¬ 
tive  filters  were  designed  to  solve  some  very  specific  real-time  engineering  problems. 
The  directions  defined  by  these  initial  applications  have  created  the  various  regimes  of 
adaptive  filtering  and  are  still  evident  in  current  applications,  even  as  a  more  general 
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outlook  is  developing.  Among  them,  linear  prediction  and  high  resolution  spectral 
analysis,  adaptive  beamforming  and  adaptive  noise  cancelling  are  most  commonly 
encountered.  Central  in  the  development  of  general  adaptive  structures  has  been  the 
**  prototype”  linear  prediction  problem.  A  large  and  rapidly  expanding  body  of  litera¬ 
ture  has  been  compiled  on  linear  prediction  and  its  application  to  speech  formant 
extraction,  geophysical  signal  deconvolution  and  high-resolution  spectral  estimation. 
The  solution  to  this  problem,  obtained  through  a  whitening  filter  in  conjunction  with 
certain  modeling  schemes  for  the  process  of  interest,  has  suggested  lattice-form  solu¬ 
tions  to  the  more  general  joint-process  problem  of  which  noise-cancelling  is  a  special 
case.  Of  particular  interest  throughout  this  thesis  is  the  class  of  least-squares  lattice 
filters  with  exact  time-update  properties  and  the  potential  advantages  they  enjoy  over 
"noisy"  gradient  solutions  both  of  direct  form  and  of  the  lattice  type. 

In  this  chapter,  the  problem  of  adaptive  filtering  is  set  in  the  proper  theoreti¬ 
cal  framework.  A  brief  presentation  of  both  the  statistical  and  deterministic  formula¬ 
tions  is  given  and  the  evolution  of  the  lattice  form  is  discussed.  Next,  three  adaptive 
noise-cancelling  realizations  were  chosen  to  undergo  comparison  tests  through  com¬ 
puter  simulations.  These  are  the  Least-Mean-Square  (LMS)  algorithm,  a  direct  imple¬ 
mentation  with  a  gradient-search  mechanism,  the  Gradient  Lattice  (GRL),  a  lattice 
type  also  with  gradient-based  adaptation  and  the  Least-Squares  Lattice  (LSL),  an 
exact  time-update  realization  not  succeptible  to  "  misadjustment"  noise.  The  reported 
differences  in  performance  represent  a  significant  contribution  of  this  thesis. 
Specifically,  the  effect  of  large  eigenvalue  spread  on  the  performance  of  the  LMS  is 
shown  to  extend  to  its  noise-cancelling  configuration.  In  addition,  the  demonstrated 
detrimental  effect  of  a  "scale-change"  in  the  reference  input  to  the  performance  of 
both  gradient  methods  is  a  new  result.  Finally,  the  simulation  demonstrating  the  suc¬ 
cessful  implementation  of  the  multichannel  LSL  filter  and  its  advantages  over  cas- 
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caded  processing  of  individual  reference  channels  is  the  first  reported  to  the  best  of 
our  knowledge.  In  general,  it  is  hoped  that  future  users  will  benefit  from  the  the  expo¬ 
sition  of  the  operational  characteristics  of  these  structures  in  their  noise-cancelling 
mode. 

2.1  A  historical  overview 

The  initial  developments  on  adaptive  filtering  proceeded  independently 
within  several  groups.  Howells  and  Applebaum  [1965]  designed  and  built  a  system  for 
antenna  6idelobe  cancellation  that  used  a  reference  input  from  an  auxiliary  antenna 
and  a  simple  two-weight  filter.  Gabor  et  al.  [1969]  designed  a  "learning"  machine,  a 
non-linear  adaptive  filter  implemented  as  a  highly  specialized  analog  computer,  capa¬ 
ble  of  approximating  a  "target”  function  by  successive  adjustment  of  its  variable 
coefficients.  The  term  "adaptive  filter"  has  been  used  [Jakowatz,  Shuey  and  White, 
1960]  to  describe  a  structure  which  extracts  an  unknown  signal  from  noise,  when  the 
signal  recurs  frequently  at  random  intervals.  Davisson  [1966]  has  described  an  adap¬ 
tive  method  for  estimating  an  unknown  waveform  in  the  presence  of  white  noise  of 
unknown  variance  and  Glaser  [1960]  has  presented  an  adaptive  system  capable  of 
detecting  a  pulse  signal  of  a  fixed  but  unknown  waveform.  Adaptive  prediction  opera¬ 
tions  to  fill  telemetry  "dropouts"  were  used  by  Balakrishnan  [1962],  An  interference 
cancelling  network  using  operator-controlled  cancellation  beams  was  designed  by 
Anderson  [1968]. 

The  work  of  Widrow  et  al.  [1960]  produced  the  first  generalized  adaptive 
structure,  based  on  the  mathematical  method  of  steepest  descent,  the  LMS  algorithm. 
The  LMS  has  been  used  as  the  central  adaptive  processor  in  applications  such  as 
adaptive  array  processing  [  Widrow  et  al.,  1967;  Griffiths,  1969;  Frost,  1972  ],  adaptive 
modeling  [  Schade,  1971  ],  adaptive  channel  equilization  [  Lucky,  1965  ]  and  adaptive 
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noise  cancelling  [  Widrow,  1975  ).  The  development  of  lattice  adaptive  filters  superior 
in  terms  of  convergence  characteristics  to  direct-form  tapped-delay-line  (TDL)  imple¬ 
mentations  such  as  the  LMS,  has  led  to  investigations  of  their  performance  as  the  cen¬ 
tral  adaptive  element.  Adaptive  lattice  structures  were  used  by  Satorius  and  Alex¬ 
ander  (1979]  in  a  channel  equilization  application  and  by  Hodgkiss  and  Presley  [1981] 
in  frequency  tracking  simulations  Their  potential  advantage  over  the  direct  form 
implementation  was  clearly  demonstrated  in  both  cases.  Lattice  structures  for  noise¬ 
cancelling  have  been  suggested  [Griffiths,  1978]  but  this  regime  of  adaptive  filtering  is 
lagging  markedly  in  terms  of  reported  results  either  simulated  or  experimental.  In  par¬ 
ticular,  the  recently  developed  deterministic  least-squares  lattice  filters  [e  g.  Lee  1980 
have  not  yet  been  seriously  considered  despite  their  exact  time-update  and  superior 
convergence  properties,  achieved  with  only  a  moderate  increase  in  algorithm  complex¬ 
ity.  These  structures  will  be  examined  in  detail  in  Section  2.6. 

2.2  Optimum  filtering  and  prediction 

The  development  of  the  LMS  and  its  successful  application  in  a  plethora  of 
experimental  settings,  has  demonstrated  that  adaptive  structures  of  considerable 
diversity  can  be  realized  through  a  central  adaptive  estimator  which  is  a  practical 
implementation  of  a  Wiener  filter.  The  theory  of  linear  Minimum  Mean-Square  Error 
(MMSE)  filtering,  developed  by  Kolmogorov  [1941]  and  Wiener  [1949],  has  long  been  a 
integral  part  of  estimation  theory.  Complete  treatments  may  be  found  in  standard 
estimation  theory  texts  [  e  g.  Van  Trees  1968].  What  follows  is  a  statement  of  the  gen¬ 
eral  problem  and  a  brief  presentation  of  the  relevant  results. 
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2.2.1  The  linear  MMSE  estimation  problem 

Given  two  discrete-time,  complex,  vector  stochastic  processes  x(n),  d(n)  we 
want  to  operate  on  the  first  process  (data)  to  obtain  an  estimate  of  the  second  process 
(signal).  This  constitutes  the  general  two-channel  (joint-process)  filtering  problem  and 
is  the  focus  of  our  interest.  In  addition,  the  one-step  prediction  problem  </(„)=  *(n) 
has  played  an  important  role  in  the  development  of  lattice-type  solutions. 

We  denote  the  estimate  as  rf(n)  and  choose  a  processor  so  that,  with 
fp(n)  =  d(n )  -  d(n),  the  quantity 

is  minimized. 

Restricting  our  attention  to  the  class  of  causal  linear  processors,  we  express  the  esti¬ 
mate  as 

d(n)  =  ^)Fj,,(n)x(n  - »') 

<-o 

The  unique  optimum  linear  MMSE  filter  may  be  obtained  through  the  projection 
theorem  as  follows: 

£(«i>(n )**(*)}  =0,  n~p  <  k  <n 

equivalently, 

E{  (rf(n)  —  ^F^,)(n)x(n_,)]J./t(^)}  _o 

i—0 

or 

Ri,(n,k)  =f]F^\n)R„(n^i,k)  (la) 

where  are  the  autocorrelation  and  crosscorrelation  functions  of  the  received 

signal  and  between  the  desired  and  received  signals  respectively: 
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R„(*,k)=E{x(n)xH(Jc)},  R*(n,k)  =  E{  d{n)xH{k)}  . 

The  resulting  MMSE  is  given  by 

EA»)  =nun{J,}E{£p(n)£pff(n)} 

=  £{£,(„)«/"(»)} 

or 

EPV)  =A«(0 )  -  &¥\n)R&(n-i,n)  (lb) 

iMi 

where  Ru  is  the  autocorrelation  of  the  desired  signal. 

The  well  known  time-invariant  Wiener-Hopf  solution,  the  continuous  version  of  which 
was  derived  by  Wiener  [1949],  is  a  special  case  of  (l)  and  holds  when  wide  -sense  sta¬ 
tionary  (WSS)  prevails  and  an  infinite  observation  interval  (p  —oo)  is  assumed. 

Equation  (la)  may  be  expressed  in  matrix  notation  as 

(2) 

where 

R„  -£{x !„:„_p|X(*[:i,_f|} 

R i*  =  E{  rf(n)^:n-p|} 

with 

=  n  );X^(n  — l),...,X^(  71  — p  )  j 

and 

F/  =  [F(°)(»)40(n)v..fjpW(n)] 

The  solution  to  the  one-step  prediction  problem,  obtained  similarly,  is  given 
by 

En(n<k)  =f^Af\n)R1J(n-i,k),  n-p  <  k  <n-l 

iM 

and  the  corresponding  MMSE  is 


(3a) 
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E;(»)  =R„t o)  -  (3b) 

Relations  (3-a)  and  (3-b),  known  as  the  Yule-  Walker  equations,  may  be  combined  to 

form  the  so-called  normal  equations: 

«[E*,0,...,0]  (4) 

with 

where  Im  is  the  m  x  m  identity  matrix. 

Remarks 

i-  An  important  property  of  the  optimum  linear  processor  is  that,  when  the 

processes  involved  are  Gaussian,  it  is  also  the  best  of  any  type  (  including 
nonlinear  )  for  minimizing  the  MSE.  The  Gaussian  assumption  is  reasonable 
for  oceanic  reverberation  as  the  nature  of  the  physical  process  that  generates 
it  invites  application  of  the  central  limit  theorem. 

ii.  The  solutions  represented  by  (2)  and  (4)  are  valid  for  general  vector,  nonsta¬ 

tionary  processes.  It  is  possible,  in  principle,  to  solve  them  directly  by  numer¬ 
ical  methods  given  ensemble  representations  of  the  processes  involved.  In 
applications  of  practical  interest  where  real-time  processing  is  of  the  essence, 
it  is  often  the  case  that  only  a  single  sample  function  of  the  underlying  ran¬ 
dom  process  is  available  and  time-average  statistics  must  be  used.  When 
WSS  holds  at  least  locally  and  ergodicity  can  be  safely  invoked,  this  will  still 
lead  to  the  optimum  solution.  In  practice,  however,  ergodicity  is  difficult  to 
prove  and  is  routinely  assumed.  In  addition,  during  periods  of  transition 
(departure  from  WSS),  the  optimality  of  the  filter  will  be  highly  questionable 
from  the  probability  theory  point  of  view. 
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iii.  Here  we  are  mainly  interested  in  solving  the  filtering  problem  However,  fast 
and  efficient  (lattice-type)  solutions  for  (2)  can  only  be  obtained  in  conjunc¬ 
tion  with  the  solution  of  the  prediction  problem.  In  that  sense,  (2)  and  (4)  are 
intimately  related 

2.3  The  Adaptive  Noise  Cancelling  (ANC)  concept 

The  set  of  circumstances  which  gives  rise  to  the  noise  cancelling  concept  is 
illustrated  in  Figure  2  la.  The  primary  channel  consists  of  the  signal  corrupted  by  a 
form  of  additive  noise  n,  and  the  reference  channel  consists  of  a  process  related  in 
some  unknown  way  to  the  primary  noise.  The  key  requirement  is  that  the  signal  be 
uncorrelated  with  both  the  primary  noise  and  the  reference  process 

R„§  =0,  o. 

We  then  have 

which  leads  to  the  noise  cancelling  solution  as  a  special  case  of  the  general  filtering 
problem  solution: 

)#.,.,(«>*)»  n-p  <k  <n.  (5) 

*M) 

Thus,  solving  the  filtering  problem  in  this  setting  is  equivalent  to  producing  the  best 
MMSE  estimate  of  the  primary  noise  process.  The  ANC  output,  obtained  by  subtract¬ 
ing  this  estimate  from  the  primary  input,  will  consist  of  the  signal  component  «  plus  a 
residual  error  n,  —  ri#. 

A  more  detailed  model  for  the  noise-cancelling  structure  is  shown  in  Figure 
2.1b.  The  "mismatch”  between  the  primary  and  reference  inputs  is  represented  by  a 
linear  transfer  function  H(z).  Additional  uncorrelated  noise  components  m,  and  m, 
are  included  in  the  two  channels.  This  particular  model,  which  is  representative  of 
many  situations  of  practical  interest,  has  been  studied  in  detail  by  Widrow  et.al. 
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[1975],  It  was  determined  that  the  uncorrelated  noise  components  have  a  deleterious 
effect  on  cancellation,  expressed  in  z-transform  notation  [Oppenheim  and  Schafer, 
1975]  by  the  following: 

P...U)  (A /,(*)  +  !  )1M,(*)+1  ] 

p^(z)  =  M.(z)  +  Mt(z)Ml(z)  +  W,(z)  (5a) 

Here,  p,„{z)  is  the  signal-to-correlated  noise  density  ratio  at  the  output,  p^(z)  is  the 
signal-to-correlated  noise  density  ratio  at  the  primary  input  and  MJz),  A/,(r)  are  the 
ratios  of  the  uncorrelated  noise  spectrum  to  the  correlated  noise  spectrum  at  the  pri¬ 
mary  and  reference  channels,  respectively.  It  is  apparent  from  (5a)  that  the  ability  of 
the  ANC  system  to  reduce  noise  is  limited  by  M,(z)  and  A/,(?)  when  these  quantities 
are  large.  When  they  are  small,  other  factors  become  important  in  limiting  ANC  per¬ 
formance.  One  such  factor  is  the  presence  of  signal  components  in  the  reference  chan¬ 
nel.  Letting  pn/(z)  denote  the  signal-to-correlated  noise  density  ratio  at  the  reference 
input  and  assuming  M„(z)  =M]{z)  =0,  Widrow  et.  al.  show  that 


P.,M 


Pr,/(2) 


(5b) 


Another  relation  of  potential  practical  significance  is  the  following,  regarding  the  out¬ 
put  noise  spectrum  S„(r), 


~  (5c) 

Equation  (5b)  indicates  that,  although  the  signal  components  in  the  reference 
input  has  a  detrimental  effect  on  cancellation,  it  does  not  render  the  ANC  operation 
useless  if  pre/(z)  is  held  to  reasonably  low  levels.  Equation  (5c)  makes  the  somewhat 
surprising  statement  that  a  high  signal- to- noise  ratio  in  the  primary  channel  will 
result  in  increased  ouput  noise  power.  This  can  be  explained  intuitively  as  follows:  If 
PPri{2)  is  low,  the  filter  will  be  "trained"  to  cancell  the  noise  rather  than  the  signal 
and  therefore  output  noise  will  be  low. 
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2.3.1  Adaptive  beamforming  via  ANC 

Conventional  arrays  are  usually  designed  under  the  assumption  that 
unwanted  noise  is  spatially  disorganized.  These  fixed-weight  arrays  suffer  a  degrada¬ 
tion  of  performance  in  the  presence  of  directional  interference  possessing  some  degree 
of  spatial  coherence.  Such  interference  may  be  due  to  natural  sources,  reverberation 
returns  or  jamming  signals.  Other  factors,  such  as  array  motion,  multipaths  and  con¬ 
stantly  changing  interference  characteristics,  contribute  to  further  deterioration. 
Adaptive  arrays,  on  the  other  hand,  are  capable  of  reducing  or  eliminating  directional 
noise  components  and  of  responding  to  changing  conditions  by  adjusting  their  pattern 
response  according  to  an  appropriately  chosen  error  criterion  In  general,  adaptive 
arrays  are  in  the  form  of  a  space-time  filter  which  can  be  implemented  through  the 
general  multichannel  filtering  structures  described  in  the  previous  section. 

The  fields  of  communications,  radar,  sonar  and  geophysics  have  each  contri¬ 
buted  to  adaptive  array  development.  The  exact  manner  in  which  control  of  the 
available  spatial  and  temporal  degrees  of  freedom  is  relegated  to  the  adaptive  proces¬ 
sor  and  the  choice  of  a  performance  measure,  depend  largely  on  the  objectives  and 
priorities  of  the  specific  application.  For  instance,  Shor[l966]  adjusted  the  array 
coefficients  to  maximize  the  signal-to-noise  ratio  (SNR)  at  the  output.  Riegler  and 
Complon[l973]  and  Widrow  [1967]  minimized  the  mean  square  of  an  error  signal  equal 
to  the  difference  between  the  array  output  and  a  reference  signal.  The  total  noise 
power  at  the  output,  in  the  absence  of  the  signal,  was  minimized  by  Lacoss[l968t, 
while  optimum  probability  of  detection  for  a  fixed  false  alarm  rate  was  the  criterion 
chosen  by  Brennan  and  Reed[l973j.  The  special  case  of  weak  signals  in  the  presence  of 
strong  interference  was  examined  by  Zahm[l973].  A  good  tutorial  introduction  to  the 
adaptive  array  problem  is  given  by  Gabriel[l976|  and  a  more  complete  treatment  is  by 
Monzingo  and  Miller[l980]. 
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Here  we  are  primarily  interested  in  signal  extraction  and  MMSE  is  the  per¬ 
formance  measure  of  our  choice.  It  is  flexible  enough  to  encompass  many  experimental 
configurations  and  it  leads  to  mathematically  tractable  results.  Our  objective  is  to 
utilize  the  general  adaptive  algorithms  which  have  been  derived  under  this  criterion  in 
conjunction  with  the  ANC  scheme,  in  order  to  eliminate  reverberation  interference 
entering  through  the  sidelobes  and  possibly  the  main  lobe  of  the  receiving  array.  It 

will  be  shown  that  selective  reverberation  cancellation  may  be  achieved  through  con- 
strained  adaptive  beamfornring 

2.3.2  Adaptive  beamforming  with  constraints 

Of  crucial  importance  in  any  adaptive  array  application  is  the  utilization  of 
all  available  a  priori  information  concerning  the  nature  of  the  signal  and  the  interfer¬ 
ence,  in  order  to  ensure  signal  preservation  at  the  output.  In  the  early  applications  of 
adaptive  processing  to  antenna  sidelobe  cancellation,  the  effects  of  signals  on  the 
adapted  response  were  commonly  neglected.  This  was  justified  by  the  assumption  that 
the  signals  of  interest  were  sufficiently  low  in  power  and/or  duration  so  that  the  adap¬ 
tive  control  algorithm  would  not  respond  to  them  (e  g.  Widrow  et.  al.[l967].)  When 
such  an  assumption  is  unrealistic,  the  response  of  the  adaptive  array  must  be  expli¬ 
citly  constrained  to  avoid  signal  cancellation.  Constraints  may  be  incorporated  into 
the  adaptive  algorithm  to  maintain  a  chosen  frequency  response  in  a  desired  direction 
(Frost[l972]).  This  is  computationally  expensive  and  operationally  inflexible  "Pilot" 
signals,  simulating  actual  signals  of  interest,  have  also  been  used  to  describe  a  desired 
main"  look  direction,  through  a  twomode  adaptation  process  [Widrow  1967],  How¬ 
ever,  accurate  replicas  of  the  desired  signals  are  not  normally  available.  An  alterna¬ 
tive  is  to  impose  preadaptation  constraints  in  the  time,  space  and  frequency  domains 
[Applebaum  and  Chapman,  1976],  Restrictions  may  be  effected  in  the  time  domain  by 
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allowing  the  array  to  adapt  only  when  no  signal  is  present;  in  the  frequency  domain 
by  allowing  adaptation  only  to  energy  received  outside  the  signal  band.  Finally, 
preadaptation  spatial  filtering  may  be  used  to  remove  the  signal,  whose  direction  is 
presumed  known,  thus  protecting  it  from  cancellation.  Such  prefiltering  may  range 
from  complete  conventional  beamforming  (Figure  2.2a)  to  simple  element-to-element 
subtraction  (Figure  2.2b).  The  former  can  be  applied  only  when  the  spatial  charac¬ 
teristics  of  the  interference,  as  well  as  the  signal,  are  known  and  constant  and  requires 
an  external  steering  mechanism.  The  latter  is  more  appropriate  when  the  directional 
characteristics  of  the  interference  are  unknown  or  variable  and  results  in  (constrained) 
spatially  adaptive  cancellation  beams.  These  two  configurations  are  compatible  with 
our  objectives  and  will  both  be  used  with  attention  to  their  relative  advantages. 

It  is  interesting  to  note  that  such  prefiltering  operations,  by  eliminating  the 
signal  of  interest  from  the  reference  channel(s),  in  effect  place  constrained  adaptive 
beamforming  in  the  context  of  ANC  and,  more  generally,  of  the  (vector)  Wiener  filter¬ 
ing  problem.  This  approach  affords  us  the  flexibility  of  using  any  efficient  solution  to 
the  general  filtering  problem  as  the  adaptive  processor  controlling  the  array  element 
weights.  It  will  be  shown  that  significant  differences  in  performance  exist  between 
seemingly  equivalent  algorithmic  implementations  and  that  the  correct  choice  is  criti¬ 
cal  to  succesful  reverberation  cancellation. 

2.4  Stochastic  lattice  solutions  for  optimum  filtering 

The  optimum  linear  filters  described  by  (2)  and  (4)  can  be  determined  by 
solving  the  matrix  equation  directly.  In  the  stochastic  approach,  the  cost  is  twofold 
First,  the  needed  second  order  statistics,  which  are  almost  always  unknown  in  prob¬ 
lems  of  practical  interest,  must  be  estimated  from  the  data  and  updated  periodically 
since  the  statistics  are  in  general  nonstationary.  The  second  step  involves  the  inver- 


30 


sion  of  the  autocorrelation  matrix  R.  In  general,  the  p  x  p  correlation  matrix  can  be 
inverted  by  standard  methods  such  as  Gaussian  elimination  and  Crout  reduction 
[Stewart,  1973],  or  more  robust  methods  making  use  of  the  nonnegative  definiteness 
and/or  symmetry  properties  of  R  such  as  Householder  transformations  or  Cholesky 
decompositions  [Lawson  and  Hanson,  1974].  All  of  these  methods  require  0(p3)  compu¬ 
tations  (multiplications  of  real  numbers).  The  computational  requirements  can  be 
substantially  reduced  when  R  is  characterized  by  certain  shift-invariance  properties 
which  lead  to  efficient  recursive  solutions  of  the  lattice  or  ladder  type. 

For  the  remainder  of  this  section,  we  will  restrict  our  attention  to  scalar 
processes.  This  was  done  for  the  sake  of  notational  simplicity  and  does  not  seriously 
affect  the  generality  of  the  presentation. 

2.4.1  Levinson’s  algorithm:  The  Toeplitz  recursion 

A  faster  method  of  inversion,  and  therefore  solution  of  the  normal  equations, 
exists  when  R  is  Toeplitz,  i.e.  of  the  form 

R=lr(« -;)].  0<  «•,;<  p  (6) 

It  is  easy  to  see  that  a  Toeplitz  structure  for  R  is  tantamount  to  the  wide-sense  sta- 
tionarity  of  the  input  process.  In  this  case,  the  shift-invariance  and  symmetry  proper¬ 
ties  of  R  lead  to  an  inversion  technique  requiring  only  0(p2)  computations.  The  pro¬ 
cedure,  developed  originally  by  Levinson  [1947]  and  rediscovered  by  Durbin  [i960],  con¬ 
sists  of  using  a  partial  solution  and  a  subsequent  correction  through  a  "reversed  time" 
auxiliary  solution.  The  same  concept  has  been  used  by  Szego  [1939]  in  the  different 
but  mathematically  isomorphic  problem  of  polynomial  approximation  on  the  unit  cir¬ 
cle.  Multivariate  extensions  of  the  Levinson  algorithm  have  been  obtained  by  Whittle 
[1963]  and  Wiggins  and  Robinson  [1965],  Toeplitz  recursions  exist  for  both  the  filtering 
and  prediction  solutions.  The  auxiliary  Toeplitz  recursion  solves  for  the  predictors 


31 


of  (4)  and  the  general  Toeplitz  recursion  solves  for  the  filter  coefficients  Fjfi  of  (2). 
These  recursions  give  rise  to  the  prediction  error  and  joint  process  lattice  structures, 
respectively.  Although  Levinson,  in  his  original  work,  was  mainly  concerned  with  the 
general  filtering  problem,  it  was  the  prediction  solution  (auxiliary  recursion)  which  has 
attracted  much  attention  in  recent  years,  in  conjunction  with  certain  modeling 
schemes  for  the  input  process.  The  general  Toeplitz  recursion  can  only  be  achieved 
through  the  auxiliary  recursion  and  thus  the  two  procedures  are  intimately  related. 
The  basic  structure  for  the  prediction  lattice  filter  is  described  by  the  following  rela¬ 
tions  (Figure  2.3a  with  k-  =  k-  =  k{.) 

«•(»)  =  e.MM  +  *,(n)r,_1(n- 1)  (7a) 

*•.(«)  =  *.*(n)«.-i(n)  +  r,-j(n— l)  (7b) 

where  x(n )  is  the  input  signal  and  e,(n)  and  r,(n)  are  the  forward  and  backward  pred¬ 
iction  error  residuals  at  stage  «. 

Levinson’s  algorithm  in  its  original  form  requires  the  external  calculation  of 
autocorrelation  lags  and  is  sometimes  referred  to  as  the  "given  correlation"  algo¬ 
rithm.  Alternative  stochastic  lattice  formulations  which  utilize  the  input  data  stream 
to  directly  estimate  the  lattice  coefficients  are  more  suitable  for  real-time  processing. 
One  such  "given  data"  structure  (GRL)  will  be  considered  here  as  a  potential  candi¬ 
date  for  our  central  adaptive  processor.  The  original  Levinson’s  algorithm  is  included 
in  Appendix  A  for  the  sake  of  completeness. 

2.4.2  Lattice  properties 

Levinson’s  algorithm  has  implications  beyond  the  original  intent  of  a  fast 
matrix  inversion  technique.  The  lattice  structures  it  has  created  possess  properties, 
summarized  below,  which  can  amount  to  significant  operational  advantages  over 
direct  TDL  implementations. 
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i.  Stage-by-stage  decoupling 

Each  lattice  stage  is  independent  of  all  other  stages.  Increasing  the  predictor 
order  is  achieved  by  adding  one  more  lattice  section,  without  changing  any  of  the  pre¬ 
vious  sections.  Thus,  given  a  lattice  implementation  of  a  predictor  of  order  p  we  have 
in  fact  an  implementation  of  all  predictors  up  to  that  order.  Stated  somewhat 
differently,  the  time-variation  of  the  finite-length  filter  takes  on  a  simple  step-function 
form  in  the  lattice,  whereby  each  lattice  section  is  turned  "on"  one  at  a  time.  This 
decoupling  or  orthogonalization  property  is  directly  responsible  for  the  improved  con¬ 
vergence  properties  of  the  lattice  (see  Simulation  1.) 

ii.  Stability  "  by  inspection " 

The  HR  (inverse)  lattice  filter  is  stable  if  and  only  if  I*, I  <1.  As  a  result  sta¬ 
bility  (or  lack  thereof)  of  the  inverse  ("synthesis")  lattice  filter  is  easily  established  by 
examining  these  reflection  or  partial  correlation  (PARCOR)  coefficients 

iii.  AR  cut-off 

For  an  input  AR  process  of  order  p,  the  kp+i,  *=1,2,...  are  zero.  This  pro¬ 
perty  can  significantly  facilitate  AR  model  identification  (see  Section  2.4.3) 

iv.  Structural  modularity 

The  lattice  forms  are  built  from  a  cascade  of  sections  with  identical  struc¬ 
ture,  a  property  attractive  from  the  VLSI  implementation  point  of  view. 

v.  Desirable  numerical  properties 

There  is  evidence  indicating  that  lattice  structures  are  less  sensitive  to 
roundoff  noise  (Markel  and  Gray,  1975;  Gray  and  Markel,  1975]  and  to  quantization 
effects  [Bruton,  1976;  Fettweis,  1974]. 

vi.  Physical  interpretations 
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The  reflection  coefficients  can  be  related  to  parameters  of  certain  models 
representing  physical  processes.  In  particular,  models  of  wave  propagation  in  a 
stratified  medium  [Thomson,  1950]  and  of  the  human  vocal  tract  [Kelly  and  Loch- 
baum,  1962]  have  been  suggested  which  have  lattice-form  configurations.  The  work  of 
Wakita  [1973]  established  a  direct  connection  between  the  reflection  coefficients  and 
the  cross-sectional  area  function  of  the  vocal  tract,  showing  that  the  vocal  tract  shape 
can  be  estimated  directly  from  the  acoustic  speech  waveform  through  linear  predic¬ 
tion.  This  parameter  estimation  technique,  coupled  with  a  recent  surge  of  activity  in 
the  area  of  speech  signal  processing,  can  account  to  a  large  extent  for  the  predomi¬ 
nance  of  the  linear  prediction  solution  in  the  literature.  On  the  other  hand,  the  joint- 
process  lattice  solution  was  largely  neglected,  leaving  direct  TDL  implementations 
(e.g.  the  LMS)  unchallenged  in  filtering  applications  such  as  channel  equilization  and 
noise-cancelling. 

2.4.3  Alternative  lattice  solutions 

The  reflection  coefficients  k(  of  the  lattice  prediction  filter  can  be  interpreted 
as  a  measure  of  the  crosscorrelation  between  et_i(n)  and  rt^(n-l).  This  idea  has 
spawned  a  class  of  lattice  algorithms  in  which  the  reflection  coefficients  are  computed 
directly  from  the  data  and  differ  only  in  the  manner  that  this  quasi-correlation 
coefficient  is  normalized.  Itakura  and  Saito  [1971]  used  the  geometric-mean  method  of 
normalization  to  compute  their  PARCOR  coefficients.  Burg  [1975],  on  the  other  hand, 
chose  the  harmonxc-mean  method  which  satisfied  his  maximum- entropy  criterion.  A 
more  general  class  of  normalization  schemes,  of  which  the  above  are  special  cases,  has 
been  defined  by  Makhoul  [1977,  1978].  We  have  chosen  the  harmonic-mean  (maximum 
entropy)  method  for  our  stochastic  lattice  implementation  (GRL)  because  its  deriva¬ 
tion  was  based  on  a  well-defined  and  intuitively  pleasing  criterion.  Morover,  it  is  the 
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method  of  choice  for  some  widely  referenced  implementations  (e.g  Griffiths  [1978], 
Satorius  [1979]).  Thus,  we  define  i,  as 


M«) 


2E{  e,(n)r,‘(n-l)} 

£{kM2+lr,(n-l)|2} 


(8) 


where  the  expectation  operator  E  is  approximated  in  practice  by  a  time  average. 

Joint-process  versions  of  the  given  data  lattice  can  be  easily  derived  by  con¬ 
sidering  the  prediction  lattice  as  a  practical  implementation  of  Wold’s  decomposition 
theorem  [Wold  1939],  In  essence,  the  lattice  performs  a  successive  orthogonalization,  of 
the  Gram-Schmidt  type,  on  delayed  versions  of  the  input  signal.  The  backward  residu¬ 
als  are  orthogonal  (£{r,r,>-^  - 1 . p)  and  constitute  an  innovations 

representation  of  the  input  signal,  forming  a  natural  orthogonal  basis  which  spans  the 
signal  subspace.  The  dual-channel  filtering  problem  can  be  solved  by  estimating  the 
desired  signal  with  a  linear  combination  of  these  orthogonal  components  A  simple 
application  of  the  orthogonality  principle  [Satorius  and  Alexander  1979]  leads  to  the 
optimum  set  of  cross-correlation  coefficients  kf  thus  forming  the  joint-process  lattice 
structure  (Figure  2.3b),  suitable  for  noise-cancelling  applications: 

iliM.-i**)  <9> 

These  given  data  structures,  by  implicitly  estimating  the  second  order  statis¬ 
tics  and  by  making  the  error  signals  available  at  each  stage  of  the  lattice,  are  more 
suitable  for  real-time  processing  than  the  original  given  correlation  Levinson  algcv 
nthm.  In  addition,  adaptive  versions  may  be  easily  derived  by  introducing  an  exponen¬ 
tial  data  window  and  an  error  gradient  criterion.  It  should  be  emphasized  that  these 
structures  are  still  based  on  a  stochastic  formulation  and  are  fundamentally 
equivalent  to  Levinson’s  algorithm.  They  are  optimum  only  when  WSS  and  ergodicity 
can  be  safely  assumed.  They  are,  however,  precursors  to  purely  deterministic  least- 
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squares  structures,  to  be  examined  in  Section  6  of  this  chapter. 


2.4.4  The  question  of  modeling 

It  is  often  useful  to  develop  a  parametric  model  for  the  behavior  of  a  signal 
or  system.  The  problem  of  determining  the  parameters  of  a  hypothesized  mathemati¬ 
cal  model  given  observations  of  the  signal  to  be  modelled  is  the  subject  matter  of 
parameter  estimation  and  system  identification.  A  variety  of  estimation  techiques  have 
been  used  in  this  context,  with  least-squares  methods  occupying  a  prominent  position 
[Astrom  and  Eykhoff,  1974;  Kailath,  1974;  Willsky,  1979], 

A  very  successful  application  of  modeling  has  led  to  the  development  of  Kal¬ 
man  filters  [  Kalman  and  Bucy,  1961]  which  have  proven  to  be  very  useful  in  control 
applications  where  internal  models  (or  state-space  descriptions)  are  relevant.  Parame¬ 
ter  identification  problems  also  arise  in  several  digital  signal  processing  applications 
where  external  models  (  input/output  descriptions)  are  appropriate,  among  them 
seismic  signal  processing  and  the  analysis,  coding  and  synthesis  of  speech.  The  latter 
application,  in  particular,  has  received  a  great  deal  of  attention  in  the  past  10  years  [ 
e.g.  Markel  and  Gray,  1976,  Makhoul,  1975],  Modeling,  in  this  context,  entails  describ¬ 
ing  the  process  of  interest  in  terms  of  the  linear  (possibly  time-varying)  system  that 
would  generate  it  when  driven  with  white  noise.  The  most  general  system  of  this  type 
is  the  autoregressive-moving  average  (ARMA)  model,  described  by  the  following 
difference  equation. 


In  z-transform  notation  [Oppenheim  and  Schafer,  1975]  we  have 

A(z)X(z)=B(z)U(z)  =>  X(z)=H(z)U(z) 


(10) 


where 


(11) 
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H(z)-B(z)/A(z) 

is  a  rational  transfer  function.  If  j-0,  then  x.  is  known  as  an  autoregressive  (AR)  pro¬ 
cess  and  the  model  is  referred  to  as  "all-pole": 

■  H(z)  =  l/A(z). 

The  all-zero’  filter  which  will  reduce  a  given  AR  process  of  order  p  to  white  noise  is 
known  as  the  inverse  or  whitening  filter  and  is  identical  to  the  prediction  error  filter  of 
order  p.  In  fact,  parameter  estimation  can  always  be  regarded  as  a  mathematical 
modeling  process,  even  if  the  primary  goal  is  not  explicitly  stated  in  this  form.  There¬ 
fore,  AR  system  identification  can  be  effected  through  the  solution  of  the  one-step 
prediction  problem.  This  procedure  is  the  basis  for  most  high  resolution  spectral  esti¬ 
mation  techniques,  where  the  spectrum  of  the  input  process  is  modeled  as  an  AR 
transfer  function.  In  the  two-channel  filtering  problem,  on  the  other  hand,  it  is  the 
cross-spectrum  between  the  two  channels  which  is  modelled  as  an  all-pole  transfer 
function  and  must  be  compensated  for  by  the  all-zero  filter  represented  by  the  lower 
section  of  the  joint-process  lattice  (see  Figure  2.1b.)  It  follows  that  the  optimum  choice 
of  filter  order  is,  in  general,  quite  different  for  the  single-channel  (prediction  error)  and 
the  two-channel  (joint-process)  filters. 

2.5  Gradient-based  time-recursive  algorithms 

The  coefficients  of  the  whitening  (prediction)  filter,  in  either  its  direct  or  lat¬ 
tice  form,  can  be  obtained  by  operating  on  a  segment  of  the  input  time  seriesuntil  a 
white-noise  sequence  (or  as  close  to  white  as  possible)  is  produced  at  its  output.  In  this 
block-  processing  mode,  a  new  set  of  coefficients  (AR  model  estimates)  is  derived  per 
data  segment.  An  alternative  approach  is  to  continuously  update  the  coefficients,  util¬ 
izing  the  unbroken  stream  of  data  and  taking  into  account  each  new  sample.  This 
requires  the  existence  of  an  updating  mechanism  with  a  decaying  memory,  whereby 
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"new"  samples  are  emphasized  over  "old".  The  resulting  filters  are  known  as  time- 
recursive  or  adaptive  and  can  be  easily  derived  from  block  structures  of  the  given-data 
type. 

2.5.1  The  LMS  adaptive  filter 

The  LMS  algorithm  uses  an  implementation  of  the  method  of  steepest  des¬ 
cent  to  provide  a  time-updating  mechanism  for  the  weights  of  the  direct  tapped-delay 
line  (TDL)  realization  of  the  linear  combiner  which  can  be  used  either  as  a  predictor 
or  noise-canceller  (Figure  2.4). 

The  "current"  weight  vector  a(n)  (Figure  2.4.b)  is  set  equal  to  the  "past" 
weight  vector  a(n-l)  plus  a  change  proportional  to  the  negative  gradient  of  the 
mean-square  error  surface  which  is  a  concave  paraboloid  with  a  unique  minimum. 

a(n)  =a(n-l) +/i(-V(n-l))  (12) 

where  p  is  the  factor  that  controls  stability  and  rate  of  convergence. 

Since  the  true  gradient  of  the  true  MMS  error  is  not  known,  it  is  replaced  by 
the  gradient  of  the  instantaneous  error  power  U*(n-l)|2.  This  "estimated"  gradient  is 
given  by  Widrow  '1975’: 

V(n-l)  =-2£/(n— l)x'(n)  (13) 

where 

£^n)  =  d(n)  -arx(n) 
xr(n)  =  [z(n-l),....,  x(n-p)} 

and  forms  the  basis  for  the  Widrow-Hoff  LMS  algorithm  of  adaptive  weight  updating: 

a(n)  =a(n— 1)  +  2/if(n-l)x'(n-l)  (14) 

It  has  been  shown  by  Widrow  [1976]  that  the  expected  value  of  the  weight 
vector  converges  to  the  true  Wiener  solution.  It  is  also  shown  that  the  algorithm  will 
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converge  in  the  mean  and  remain  stable  as  long  as 

1A»1»  >»  >0  (15) 

where  Xm„  is  the  largest  eigenvalue  of  the  autocorrelation  matrix  R.  In  addition,  the 
time  constant  for  the  "learning  curve"  corresponding  to  the  natural  mode  of  the 
signal  is  related  to  the  »'*  eigenvalue  of  R  by 

r,=l/4/iX,  (16) 

Morover,  the  ratio  of  the  excess  noise  power  due  to  weight  misadjustment  over  the 
estimation  error  power  is  given  by 

M=///rR=/i  p  X.„  *i(i-)o(  (17) 

Relations  (16)  and  (17)  portray  the  major  drawback  of  the  LMS.  When  the 
matrix  R  has  a  large  eigenvalue  spread,  convergence  is  controlled  by  the  smallest 
eigenvalue  and  misadjustment  by  the  largest  one  The  result  is  slow  adaptation  cou- 
pled  with  high  misadjustment  noise. 

It  is  convenient  to  introduce  a  new  parameter  aLUS  [Griffiths,  1976]  such  that 

V  =  qlms/pE  o 

where  E'0  is  the  power  of  the  reference  sequence  x(n).  Then  (15)  becomes 

0<<XUIS<PE  oAbu  (18) 

Since  pE'0  >Xmu  it  follows  that  convergence  is  always  achieved  for  0  <aLMS  <  1. 

In  the  algorithm,  £g  is  replaced  by  the  estimate 

Eo(n)  —(l~°iMs)Eo(n—- l)  +auwslx(n)l2  -  (19) 

2.5.2  Summary  of  the  LMS  algorithm 

The  scalar,  complex  joint-process  gradient  transversal  filter  (LMS)  algorithm 
is  summarized  as  follows  (see  Figure  2.4a) 
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Initialization  (k=0,l,...,p) 


/*(—!)  -  0  (20a) 

El{~ l)  =  eLMS,  tLMS  -  El(®)  (20b) 

*(-*)  =  x  (-P-1)  =  o  (20c) 

*/N)  =0  (20d) 

Time  update  (n>0  and 


/.(")  -  /.("-D  -  (20e) 

eptn)  “  £/*("M«-*)  +  An)  (20f) 

EKn)  =  (!  -Qiws) -^»(n-l)  +  QLMS^x{nP-  (20g) 

Slightly  simplyfying  (20a)-(20g),  the  one-step  forward  linear  predictor  com¬ 
plex  gradient  transversal  filter  algorithm  can  be  written  as  (see  Figure  (2.4b)  ) 

Initialization  (k=4,2,...,p) 


at(-i)  =  0  (21a) 

=  tLus,  €lms  -  E'M  (21b) 

*(“*)  =  x (  p  1 )  =  0  (21c) 

M"1)  =  0  (21d) 

Time  update  (n>0  and  Jfc=l,2,...,p) 

a*(n)  =  «*(n-l)  -  -  ep(n-l)x*(n-*-l)  ( 21e ) 
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er(n)  “  £“*(»)*("-*)  .  «♦  ■  1  (2 If) 

k  — 

A'M  -  (1  ~qlus)  El{n-l)  +  aiW5|j(n)l2.  (2lg) 


Remarks 

Oceanic  reverberation  is  characterized  by  the  rapid  onset  of  components 
differing  widely  in  intensity.  The  presence  of  such  components  will  most  likely  create 
this  condition  of  high  eigenvalue  spread  and  at  the  same  time  the  need  for  extremely 
fast  adaptation.  It  has  been  shown  in  frequency  tracking  [Hodgkiss  and  Presley,  1981’ 
and  adaptive  equilization  [Satorius,  1979]  applications  that  adaptive  forms  of  the  lat¬ 
tice  structure  possess  significant  performance  advantages  over  the  LMS,  resulting 
mainly  from  the  orthogonalization  and  decoupling  property  between  lattice  stages.  In 
effect,  each  stage  of  the  lattice  is  asigned  to  a  single  mode  and  adapts  independently 
of  all  other  stages.  It  will  be  the  intent  of  Simulation  1  to  confirm  this  in  a  noise¬ 
cancelling  application. 


2.5.3  The  gradient  lattice  filter  (GRL) 

The  same  concept  of  an  estimated  error  gradient  may  be  used  to  update  the 
reflection  coefficients  of  a  lattice  filter.  Again,  the  true  gradient  with  respect  to  the  i, 
of  the  sum  of  the  mean  square  error  residuals  £{U,-(n)l2  +lr,-(n)l2}  is  replaced  by  the 
noisy  gradient  of  the  sum  of  the  instantaneous  square  error  residuals  le,(n)l2 +lr,(n)l2 
The  noisy  estimate  which  can  be  obtained  directly  from  (8),  can  be  used  to  estimate 
the  "next”  k{  as  follows:  [Makhoul,  1978,  Griffiths,  1977] 


Mn)  =  + 


ti(n~ l)r,-i(n— 2)  +  c.^n-l^n-l) 

A(") 


1<*'<P  ~ 1 


(22a) 


where  Dt(n)  is  computed  recursively  from 


A  (n)  =(1-aCRi)A(«-l)  +le,_)(n-l)l2  +lr,_j(n— 2)l2,  l<«'<p-l 


(22b) 


41 


The  product  ttGRLDi(n)  where  0<acffi<l,  represents  an  exponentially  averaged  esti¬ 
mate  of  the  power  level  at  the  stage  of  the  lattice.  The  normalization  by  the  factor 
Di(n)  is  intended  to  maintain  the  same  time  constant  of  adaptation  at  each  stage  of 
the  lattice.  Thus,  the  overall  convergence  rate  should  be  unaffected  by  the  eigenvalue 
spread  of  R. 

The  optimum  set  of  the  "joint-process"  coefficients  kd  may  be  updated  adap¬ 
tively  as  follows. 


k/(n) 


kf(n-l) 


e,(n— l)r,-(n— 1) 
+  G,-(n) 


(23a) 


where 


G,(n)  * *f  tr,-(n  — l)K  1<*<P  (23b) 

Only  preliminary  results  are  available  concerning  the  convergence  properties  and  the 
misadjustment  noise  of  the  gradient  lattice  [Griffiths  and  Medaugh,  1978]. 


2.5.4  Summary  of  the  GRL  algorithm 

The  scalar,  complex,  joint-process  gradient  lattice  algorithm  (GRL)  is  sum¬ 
marized  as  follows  (see  Figure  2.3b)  with  k<  «  k-  *=  Jb/*). 

Initialization  (i  =0,1,.. .,p) 


M-i)  =  o  ,  .Vo 

(24a) 

1)  =  Ef(~ 1)  =  eJCGL  >  (JCGL  =  £o(0) 

(24b) 

£|'(— 2)  =  IjCGL  . 

(24c) 

«.(-!)  =  r.(-l)  =  0 

(24d) 

r,(-2)  =  0  ,  «Vp 

(24c) 

tf-1)  =  0 

(24f) 

42 


e^-1)  “  0 

(24g) 

Time  update  (n  >  0) 

e.(")  =  M«)  -  *{*) 

(24h) 

ei,(n)  m  d(n) 

(24i) 

Order  update  (i  =0,1,...,  p) 

M«)  -  *.(»  i)  2ojcgl 

EU(n- 2)  +  E/_,(«-l) 

(24j) 

•  le,(n-l)r,^(n-2)  +  r,*(n-l)e,_1(n-l)l  ,  «> 0 

«.(«)  *  +  ^(n)r,_J(n-l)  ,  »V0 

(24k) 

ri(n)  =  r.-i(n-l)  +  *.*(«)«, n(n)  ,  «V0 

(241) 

£;L,(n-l)  =  (1  -o,ccJ^(n-2) 

(24m) 

+  qjccl  1  •'i— 1(«  — l)l2  ,  «V0 

i?-i  ( n  )  =  (1  —  QJCGl)Ef-i  (n_ 1) 

(24n) 

+  aJccJei-i(n)l2  .  «yo 

^")  -  *'("  !)  £^}  e/(n— l)ri*(n— 1) 

(24o) 

ei‘(n)  =  e/-i(n)  +  *r(n)ri(n) 

(24p) 

Predictors  (6/’>(n)  =  «$(n),  1  =0,1, ...,p) 

a/')(n)  =  *,(n) 

(24q) 

a/°(»)  =  «/,-0(")  +  *,(n 

(24r) 
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2.6  Deterministic  least-squares  lattice  structures 

The  optimality  of  the  filters  discussed  in  the  previous  section  rests  on  the 
assumptions  of  stationarity  and  in  the  given  data  case,  ergodicity.  During  transition 
periods  when  the  statistics  of  the  signal  change  drastically,  all  such  solutions  are 
clearly  suboptimal.  In  addition,  the  misadjustment  noise  of  the  gradient  methods  is 
often  a  limiting  factor  on  the  algorithm  convergence  speed.  An  alternative  to  those 
stochastic  approximation  methods  is  the  deterministic  least-squares  approach.  When 
(WSS)  stationarity  does  prevail,  the  solutions  derived  through  this  approach  perform 
identically  to  their  stochastic  (time-average)  counterparts.  However,  they  enjoy  an 
important  advantage  in  that  they  continue  to  be  optimum  (in  a  deterministic  least- 
squares  sense)  even  during  abrupt  changes  in  the  signal  statistics.  This  fact,  together 
with  the  exact  time-update  property  possessed  by  one  particular  realization,  contri¬ 
butes  to  transient  behavior  superior  to  the  gradient  solutions  and  makes  the  least- 
squares  structures  an  excellent  candidate  for  our  choice  of  adaptive  algorithm. 

The  direct  form  implementation  of  the  exact  least  squares  algorithm  was 
first  presented  by  Morf  et  al  [1977],  The  lattice  recursions  were  first  reported  in 
[Morf, Lee, Nichols  and  \iera  19(7,  and  [Morf, Viera, Lee  1977].  A  comprehensive  presen¬ 
tation  of  the  concept,  including  Hilbert-space  approach  derivations,  joint  process 
extensions  and  computer  simulations  demonstrating  the  tracking  capability  of  these 
algorithms,  is  given  by  Lee  [1980],  Satorius  and  Pack  [1979]  and  Falconer  and  Ljung 
[1978]  reported  superior  performance  of  the  exact  least  squares  algorithms  in  a 
channel-equilization  application  as  did  Hodgkiss  and  Presley  [1982]  in  simulated  fre¬ 
quency  tracking  situations  and  Nehorai  and  Morf  [1982]  in  a  line-enhancement  con¬ 
text.  Completely  absent  from  the  literature  are  reports  on  the  performance  of  these 
structures  in  their  noise-cancelling  configuration. 
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2.6.1  The  DLS  solution  to  filtering  and  prediction 

The  deterministic  least-squares  (DLS)  problem  is  formulated  as  follows:  Given 
two  dicrete  time  series 

{  *(")>  <*(») )  .  W,  <  n  <Nj 

of  dimension  m  and  /  respectively,  we  form  the  estimate  of  d(n) 

d(n)  =  &,lp)(”)r(n-;) 

=  F/'(n)x 

and  with 


£p..,/(")  =  rf(n)  ~  j(n) 

we  seek  to  minimize 

=  <r££p,.,/(n)£  *.-,/(«) 

* 

where  i,f  will  be  specified  later  (N,  <  i  <Nf  ,  i  <  f  <  N,  .) 

The  deterministic  least-squares  (DLS)  solution  to  the  filtering  problem  is 
obtained  through  the  same  arguments  used  in  the  stochastic  case  and  is  given  by. 

Fft( «\/)R,, ,(«,/)  -R*.,(i,/)  (25a) 

where 

=  Xp(,\/)X  "(.,/) 

with 


XP(.J) 


*(«)  •  •  •  *(/) 
x(t-f)  ■  ■  ■  z(/-p) 


R*.,(«\/)  =  d|f:/|X"(t\/) 
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Fj  -  (  FW’  FW- ....  Fp) } 

The  minimum  least-squares  error  is 

E/(«,/)  =min{fj^»,/)  (25b) 

=  ^.oo(«,/)-FprR2.p(«,/) 

Similarly,  the  DLS  one-step  predictor  is  given  by: 

A„r(«,/)R,, ,(«,/)  =  [  E;(.\/),  0,  0  ]  (26) 

where 

A'/(*>/)  =  (  Im,  ■■■,  } 

Ep(«,/)  —FlO0(i,f)  —  ApQ,  p(i ,f ) 

Q*,J>  =  ^*(n)X|^-l:.-p| 

ll 

Note  that  the  exact  form  of  Xp(iJ)  depends  on  the  choice  of  »  and  /.  The 
following  choices  are  most  commonly  made,  resulting  in  sample  correlation  matrices  R 
with  different  symmetry  and/or  shift  invariance  properties.  Such  properties  are  cru¬ 
cial  to  the  development  of  algorithmic  implementations  of  the  above  solutions. 


(1)  Prewindowed  case 


f-N, 

(2)  Non  windowed  case 

i=Ni+P,  f  =  Nf 

(3)  Pre-  and  postwindowed  case 

«  =  Nf,  f  =  Nf  +  p 


Remarks 

An  interesting  parallelism  exists  between  stationarity  conditions  in  the  sto¬ 
chastic  formulation  and  the  various  windowing  schemes  in  the  deterministic  least- 
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squares  development.  For  example,  pre-  and  post-windowing  results  in  a  sample  corre¬ 
lation  matrix  R  which  is  Toeplitz,  in  direct  analogy  with  wide-sense  stationarity.  The 
Toeplitz  recursion  is  equally  applicable  here  and  is  known  as  the  "autocorrelation" 
method  of  linear  prediction  in  the  speech  processing  literature  The  non-windowed  case 
leads  to  an  R  which  is  symmetric  but  not  Toeplitz  and  is  known  as  the  "covariance" 
method.  (Clearly  such  characterizations  are  unjustified  from  a  statistical  point  of 
view.) 

The  prewindowed  case  is  of  particular  interest  It  leads  to  a  matrix  R  which 
is  not  Toeplitz  but  close  to  it  and  has  shift  properties  which  lead  to  an  algorithmic 
solution  with  exact  time-update.  The  statistical  analogue  would  be  a  process  which  is 
not  wide-sense  stationary  but  close"  to  it,  in  a  sense  not  explicity  defined  in  standard 
probability  theory.  This  creates  the  possibility  of  treating  processes  with  various 
degrees  of  stationarity  through  the  same  efficient  algorithmic  solutions  once  the 
appropriate  structure  for  the  true  autocorrelation  matrix  has  been  established. 


2.6.2  The  "  pre-windowed"  DLS  lattice 

In  this  case,  the  sample  correlation  matrix  R  of  order  p  is  of  the  form 

where, 


X,.*  = 


*(0)  •• 
0 

0  0 


*{p) 


:(N) 


*(0)  •  '  '  x(N-p) 

Note  that  Rp  N  is  not  Toeplitz  but  it  consists  of  a  product  of  tw-o  (upper-  and 
lower-triangular)  Toeplitz  matrices  and  might  be  thought  of  as  being  " near-Toeplitz" 

In  fact,  it  falls  within  the  class  of  the  so-called  a— stationary  matrices  which  can  be 
expressed  as 


R-t  +  5>,.u,.u/*, 


ox  =  ±1 


(27) 


47 


where  T  is  Toeplitz  and  the  {U,}  are  upper  triangular  Toeplitz  matrices  [Friend- 
lander,  1979],  These  matrices  are  also  known  as  "low  shift  rank”  matrices  [Morf, 
1974;.  For  such  matrices,  it  is  possible  to  derive  generalized  Levinson-type  recursions. 
Specifically  in  the  pre-windowed  case,  R,,„  is  stationary  with  a  =  2  and  satisfies  the 
following  recursive  identities: 


order-  and  time-update 


order-update 


time-update 


where 


R.,(A0 


R.-ocW  Q ,,P{N) 

Q,t,p{N)  R.^N-i) 


R,.P^(N)  V,  P(N) 


R*<AN)  -R*.P(JV-1)  +  X|W_J):A,|X(^_J,. 


:W| 


(28a) 


(28b) 


(28c) 


y*AN)  =  Ei(f»-p)x^_1:,_p) 

n  ^ 

and  all  other  quantities  as  previously  defined.  These  identities  form  the  basis  for  the 
exact  order-  and  time-update  recursions  for  the  pre-windowed  DLS  prediction  lattice. 
The  DLS  joint-process  lattice  can  be  obtained  through  an  embedding  technique, 
whereby  the  joint  sample  correlation  matrix  R 2iP(N)  is  defined 


as 


where 


R,,p  =ZP(N)Z»(N) 


Z  ,{N)  = 


<*(°)  •  •  •  d{p)  ■  •  •  d(N) 

*(0)  •  ’  '  x(p)  ■  •  •  x(N) 

0  . 

0  0  x(0)  ■  ■  ■  x(N-p) 


Then,  the  filtering  normal  equations  can  be  re-expressed  in  terms  of 


as  follows: 
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F£Rlt,(N)  -  [  E*(N),  0 . 0  J  (29) 

where 

HW  -  [  I,,  F,r(AO  ] 

R,,p  has  shift  properties  similar  to  RItP.  The  two  sets  of  recursions  can  be 
combined  to  concurrently  solve  the  forward  prediction,  backward  prediction  and  filter¬ 
ing  problems,  giving  rise  to  the  complex,  multichannel  DLS  joint-process  lattice.  The 
needed  algebraic  manipulations  are  presented  in  Lee  [l980j  and  Satorius  and  Pack 
[1979], 

Exponential  weighting 

An  adaptive  version  of  the  DLS  lattice  can  be  easily  derived  by  introducing 
an  exponential  "window"  in  the  error  norm.  Specifically,  one  may  choose  to  minimize 

% 

where  X  is  a  constant  <1,  so  that  past  errors  will  have  a  smaller  influence  on  the  esti¬ 
mate.  This  constant  "fade"  factor  is  easily  incorporated  into  the  lattice  recursions 
and  allows  the  algorithm  to  compensate  for  time  variations  in  the  input  signals.  The 
"effective”  number  of  points  used  in  this  case  may  be  defined  in  terms  of  a  "time  con- 
stant"  as  follows: 

X(*-)=exp— { (Af-*)/JV„  ] 

Therefore, 

W.j-minlW.-ji} 

In  our  implementation  X  =1  —aLsl. 

2.6.3  Summary  of  the  LSL  algorithm 

The  scalar,  complex,  joint-process,  exponentially  weighted  prewdndowed  least 
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squares  lattice  algorithm  is  summarized  as  follows  (see  Figure  2.3b):  (The  vector  ver 
sion  is  presented  in  Appendix  B). 

Initialization  (i  =0,  1,...  ,  p) 


r,(-l)  =0  ,  iVp 

(30a) 

1)  —  €lsl  >  clsi  =  0  001  and  iVp 

(30b) 

l)  =  0  ,  15^0 

(30c) 

=  0  ,  0<k<i-l  ,  iVO,  and  i 

(30d) 

Tf.-i(— 1)  =0  ,  «Vp 

(30e) 

O 

II 

? 

(30f) 

rime  update  (n  >  0) 

e.(n)  =  r,(n)  =  z(n) 

(30g) 

E„(n)  =  £J(n)  =  (l  —  Qisi)  ET,{n— 1)  +  |z(n)|2 

(30h) 

l-i(n)  =  0 

(30i) 

«£,(«)  =  d(n) 

(30j) 

Order  update  (i  =0,1,. ..,p) 

A,(n)  *=  (1  —  o iSL)  A,(n— 1) 

(30k) 

e1-H(n)r1',(n-l) 

1  -  7,-^n-l)  ’ 

*<(")  =  A,*(n)  /  EU(n)  ,  «> 0 

(301) 

*.r(n)  =  A,(n)  /  E/L^n-l)  ,  i*0 

(30m) 
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«.(")  “  <<h(»)  +  *.rMr,-,(n-l)  ,  |V0 

(30n) 

r.(n)  “  »7h('»-1)  +  **(")«, -t(n)  ,  «V0 

(30o) 

£/(")  “  EUin)  -  lA,.(n)|2  /  £/_(n-l)  ,  ,yo 

(30p) 

£/(")  -  ~  1  A,(n)l2  /  EU(n)  ,  ,y 0 

(30q) 

%_,(«)  -=  7,--b(b)  +  lr,_^(n)I2  /  E^(n)  ,  .yo 

(30r) 

V(«)  =  (1  -  *i5i)A/(n— 1)  -  ^(")r/(") 

1  -  7<- 1(") 

(30s) 

Jb/(n)  -  - 

#(.») 

(30t) 

ef( n )  “  «/-i(n)  +  »)»■((«) 

(30u) 

Predictors 

a/’)(n)  =  *,r(n) 

(30v) 

-  *,?(") 

(30w) 

o/'Hn)  =  «/•-*)(«)  +  Jfcf(„)  */4-*>(n-l) 

(30x) 

1</<I-1 

*/f,(n)  =  +  A/( n )  a/'H)(n) 

(30y) 

The  likelihood  variable 

The  parameter  which  controls  the  adaptation  response  of  the  LSL  is  7 

defined  by 

(31) 

It  can  be  shown  [Lee,  1980j  that 

0  <  Ix.p  <  1. 

7  has  an  interesting  statistical  interpetation  as  the  likelihood 

of  deviation  of  the  p 
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succesive  samples  from  a  joint  Gaussian  distribution  parameterized  by  R,  p.  [Morf  and 
Lee,  1978].  It  can  be  expected  to  be  ~  0  except  when  the  signal  statistics  undergo 
drastic  changes,  in  which  case  it  will  approach  unity. 

Note  that  in  the  LSL  recursions  (Eq.  30s)  the  "  crosscorrelation"  quantities 
A,-  are  weighted  by  a  gain  factor  equal  to 

1 

1  -  i 

which  replaces  the  "smoothed"  power  estimates  of  the  stochastic  gradient  algorithms. 
Clearly,  when  7, --4  ~  1  the  algorithm  will  undergo  extremely  fast,  almost  instantane¬ 
ous,  adaptation.  It  is  the  presence  of  this  optimum  gain  factor  which  sets  the  "exact" 
least  squares  lattice  apart  from  the  suboptimum  gradient  approximations  and  results 
in  its  superior  transient  response  (see  Simulation  2.) 

On  the  choice  of  adaptation  coefficient 

The  proper  choice  of  adaptation  coefficient  can  be  made  by  considering  the 
more  general  question  of  resolution  vs  statistical  reliability.  If  WSS  holds  for  all  time, 
then  a  —0  is  the  optimum  choice  and  the  time-invariant  Wiener-Hopf  solution  is 
obtained.  In  cases  of  practical  interest,  local  WSS  if  more  often  encountered  and  the 
desired  resolution"  is  controlled  by  the  time  scale  of  variation  of  the  signal  statistics. 
By  that  it  is  meant  that  one  must  refrain  from  averaging  over  a  number  of  distinct 
statistical  events.  Therefore  the  lower  limit  of  the  adaptation  coefficients  (or  the 
upper  limit  of  the  allowable  averaging  interval)  is  set  by  the  stationarity  properties  of 
the  input  process.  In  the  gradient  descent  algorithms  (LMS,GRL)  the  lower  bound  of 
a  is  imposed  by  the  misadjustment  noise,  caused  by  excessively  large  steps  in  the  gra¬ 
dient  search  process  (e.g.  equation  (17)).  In  the  case  of  the  exact  deterministic  least- 
squares  solution  (LSL),  there  is  no  misadjustment  error.  The  lower  bound  of  ocLSL  is 
strictly  a  question  of  statistical  reliability. 
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A  complete  statistical  analysis  of  the  lattice  structure  has  not  been  pursued 
in  this  thesis.  However,  insight  may  be  gained  by  considering  the  variance  of  the  fol¬ 
lowing  crosscorrelation  estimate  between  two  bandwidth-limited  white  noise  processes 
x(n)  and  y(n): 

R>t(k)  =  T7— 

in 0 

where, 

B  ~  qlsl) 

The  variance  of  the  estimate  is  given  by  [Bendat  and  Piersol,  1971': 

Vflr&,(*)]  =  /?„(*)) 

The  above  indicates  that  large,  spurious  crosscorrelation  values  can  be  gen¬ 
erated  between  two  uncorrelated  processes  (/?„(*)  =0),  for  sufficiently  small  Ntjr 
and/or  large  autocorrelations  within  the  two  processes.  Although  direct  comparison 
between  this  simple  result  and  the  behavior  of  the  lattice  parameters  for  arbitrary 
signals  is  not  possible,  the  LSL  joint-process  filter  has  been  observed  to  obey  this  gen¬ 
eral  rule  Specifically,  the  filter  has  shown  a  tendency  towards  cancellation  of  the  pri¬ 
mary  channel  when  the  primary  and  reference  channels  consisted  of  independent  white 
noise  sequences,  for  aLSL  >0.02.  This  empirically  derived  value  will  be  used  as  the 
upper  bound  for  throughout  this  thesis. 

2.7  Computer  simulations 

The  exposition  of  the  theoretical  foundation  of  the  LMS,  GRL  and  LSL  filters 
in  the  preceding  sections  of  this  chapter  enables  us  to  make  certain  conjectures 
regarding  their  operational  characteristics.  First,  one  may  expect  the  LMS  to  be  suc- 
ceptible  to  the  "eigenvalue  spread"  condition.  On  the  other  hand,  a  lattice  filter  (GRL 
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or  LSL)  should  not  be  affected  owing  to  their  orthogonalization  (stage-by-stage  decou¬ 
pling)  property.  Secondly,  the  stochastic  approximation  solutions  (LMS  and  GRL)  are 
clearly  suboptimum  during  transient  conditions  (i.e.  departures  from  WSS).  The  exact 
deterministic  LSL  solution,  on  the  other  hand,  remains  optimum  (in  a  least-squares 
sense,  not  in  a  stochastic  MMSE  sense)  during  such  transient  states  and  may  be 
expected  to  perform  better.  The  simulations  of  this  section  were  intended  to  provide 
support  for  these  mathematical  conjectures.  They  are  a  representative  subset  of  a 
large  number  of  simulations  performed  during  the  course  of  this  research  project  [e.g. 
Alexandrou  and  Hodgkiss,  1982d,e],  Simulations  1  and  2  were  designed  to  bring  out 
the  operational  characteristics  of  the  three  chosen  structures  (LMS,  GRL,  LSL)  under 
transient  conditions  and  in  the  presence  of  the  "eigenvalue  spread"  and  "scale 
change  effects,  respectively;  the  third  to  demonstrate  the  potential  advantage  of  the 
multichannel  lattice  filter  over  a  cascade  of  single  channel  filters,  when  multiple  refer¬ 
ence  signals  are  available.  The  scalar,  complex  versions  of  LMS,  GRL  and  LSL  as  well 
as  the  vector,  complex  LSL  have  been  programmed  in  Ratfor  [Kernighan,  1975,  and 
used  throughout  this  work. 

For  the  sake  of  clarity,  simple  complex  exponentials  in  additive  white  noise 
were  used,  which  undergo  frequency  and/or  magnitude  "steps"  at  prescribed  points 
during  the  iteration.  However,  the  simulated  transients  do  represent  conditions  which 
are  likely  to  be  encountered  in  real  reverberation  signals.  Therefore,  these  results  are 
pertinent  to  the  performance  of  the  chosen  algorithms  in  a  reverberation  cancellation 
context.  These  noise-cancelling  simulations  have  not  appeared  previously  and  are  a 
major  contribution  of  this  thesis. 

2.7.1  Simulation  1  :  Eigenvalue  spread 

The  objective  here  is  to  examine  the  performance  of  the  LMS,  GRL  and  LSL 
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in  a  noise-cancelling  setting,  where  the  reference  signal  components  have  widely 
separated  power  levels.  Both  main  (primary)  and  reference  channels  consist  of  a  pair 
of  complex  exponentials  corrupted  by  white  noise.  The  reference  channel  signals 
undergo  a  frequency  "step"  at  iteration  number  1536.  Table  2.1,  below,  fully  describes 
the  simulation. 


Simulation  1 

before  step 

after  step 

main 

reference 

main 

reference 

CJi 

7 r 

JL 

JL 

l 

8 

4 

8 

8 

SNR^db)  . 

20 

10 

20 

10 

(jJrt 

6tt 

7  IT 

67T 

w  2 

8 

8 

8 

8 

SNR^db) 

20 

20 

20 

20 

Table  2.1  Summary  of  Simulation  1 


The  main  and  reference  channels  are  displayed  in  Figure  2.5a  and  2.5b 
respectively,  in  the  form  of  a  Range-Doppler-Map  (RDM)  created  by  taking  successive, 
50%  overlapped  FFTs  of  the  two  time-senes.  A  filter  order  of  p=3  was  used  and  the 
adaptation  parameters  were  chosen  such  that  the  response  of  all  three  algorithms  to 
the  high  power  sinusoid  was  similar.  The  joint-process  error  outputs  produced  by 
LMS.GRL  and  LSL  are  shown  in  Figures  2.5c, d,e.  It  is  evident  that  the  LMS  is  lagging 
substantially  in  the  cancellation  of  the  lower  power  sinusoid  compared  to  the  two  lat¬ 
tice  filters.  The  LSL  performs  marginally  better  than  the  GRL  which  is  consistent 
with  previous  findings  in  a  frequency  tracking  application  [Hodgkiss,  1981.] 

Remarks 

The  eigenvalue  spread  condition  is  likely  to  exist  in  a  real  composite  rever- 
beration  signal  since  the  surface,  volume  and  bottom  components  are  normally 
characterized  by  large  intensity  disparities.  The  LMS  is  therefore  a  poor  candidate  for 
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our  purposes. 

2.7.2  Simulation  2  :  Scale  change 

The  fundamental  difference  between  the  stochastic  approximation  realiza¬ 
tions  (e.g.  LMS,  GRL)  and  the  deterministic  least-squares  filter  (LSL)  is  the  fact  that 
the  latter  continues  to  be  optimum  within  its  design  even  during  abrupt  changes  in 
the  characteristics  of  the  signal(s).  The  LMS  and  GRL  are,  in  such  cases,  clearly 
suboptimum  in  either  the  stochastic  (MMSE)  or  the  deterministic  least-squares  sense. 
This  simulation  was  intended  to  bring  out  the  specific  consequences  of  this  fact.  The 
complex  exponentials  in  the  reference  channel  perform,  in  addition  to  the  frequency 
step,  undergo  instantaneous  magnitude  "jumps"  of  10,20  and  30  db.  The  simulation  is 
described  in  Table  2.2. 


Simulation  2 

before  step 

after  step 

main 

reference 

main 

reference 

rj. 

7T 

7T 

7T 

7T 

8 

4 

8 

8 

SNR^db) 

20 

10 

20 

20,30,40 

dir 

7  77 

6  7T 

67T 

8 

8 

8 

8 

SNR^db) 

20 

10 

20 

20,30,40 

Table  2.2  Summary  of  Simulation  2 


The  filter  outputs  for  these  three  cases  are  shown  in  Figures  2.6  and  2.7a. 
The  magnitude  plots  were  artificially  ”  clipped”  for  display  purposes. 

The  GRL  proved  to  be  extremely  sensitive  to  this  "scale  change"  effect  and 
became  increasingly  unstable,  failing  completely  to  recover  from  the  transient  in  the 
40  db  case.  The  LMS  is  similarly  succeptible  although  it  appears  to  be  considerably 
more  robust.  The  LSL  was  able  to  adapt  immediately  to  all  magnitude  jumps.  The 
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likelihood  variable  can  be  seen  (Figure  2.7b)  to  approach  unity  as  the  severity  of  the 
discontinuity  inceases.  The  crosschannel  coefficients  (tj)  of  the  GRL  and  the  LSL  exi- 
bit  a  similar  behavior  (Figure  2.8). 

It  is  likely  that  the  specific  reason  for  the  instability  of  the  gradient  filters 
lies  in  the  "smoothed"  power  estimate  (Eqs.  (19),  (22b),  (23b))  used  to  normalize  the 
coefficients  of  these  filters.  This  low-pass  filtered  power  estimate  does  not  change  fast 
enough  to  compensate  for  large  instantaneous  power  jumps. 

Remarks 

The  "scale  change"  condition  is  also  likely  to  be  in  effect  for  real  reverbera¬ 
tion  signals.  For  instance,  the  onset  of  surface  reverberation  in  a  volume  reverbera¬ 
tion  backround  can  readily  amount  to  an  abrupt  increase  in  intensity  of  the  order  of 
the  simulated  magnitude  jumps. 

2  7  3  Simulation  3:  Multichannel  vs  cascade 

In  situations  where  spatially  adaptive  beamforming  is  in  order  and  direct 
control  of  the  spatial  transfer  function  is  desirable,  the  need  for  the  true  vector  (mul¬ 
tichannel)  filter  is  obvious.  If,  however,  the  primary  channel  interference  simply  con¬ 
sists  of  a  summation  of  components  each  of  which  can  be  "tapped"  individually,  one 
might  argue  that  a  valid  alternative  to  the  multichannel  structure  is  a  cascade  of 
scalar  (single-channel)  filters  where  a  single  interference  component  is  used  as  the 
reference  channel  at  each  stage  of  the  cascade  (Figure  2.9).  This  simulation  will  show 
that  although  this  may  in  fact  be  the  case  when  the  interference  components  are 
mutually  uncorrelated,  in  general  the  multichannel  filter  is  the  better  choice  because 
of  its  ability  to  compensate  for  interference  components  correlated  between  all  refer¬ 
ence  channels  (and  not  present  in  the  primary  channel).  In  order  to  clarify  this  point, 
consider  the  case  of  a  scalar  primary  process  d(n)  and  a  twochannel  reference  signal 
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consider  the  case  of  a  scalar  primary  process  d{n)  and  a  two-channel  reference  signal 
x(n), 


The  matrix  form  of  the  solution  is  (see  eq.  2.2) 

FrR„  -R* 

For  the  sake  of  simplicity  consider  only  the  O'*  order  solution  (no  time  lags  involved.) 
For  that  case,  we  have: 


R„=£{ 


x 

x 


i(n) 

2(") 


Cl’(") 


z2'(n)  ]  } 


Rx1xl{n<n)  R,il2{n,n) 
■R*2,/n’n)  Rx^lJ^n  in) 


R* 


£{</(n)*,(n)}  £{rf(n)z2(n)}  ]  -  (  R*I(nln) 


Rixj,n,n)  ] 


Therefore, 


Fr=[/f  /!) 


or, 


[  ft  n 


Rxlxl(n,n)  Rti,£n,n) 

Rxi.1{n,n)  i?JaIj(n,n) 


=  [  Rix,(n,n) 


Rixz{n,n)  ! 


ftRxx  *,(”.«)  + /l-RxjxJ”  >”)  =  Riix(n  ,n) 

ftRxxx£n,n)  +  fiRx^^n  ,n)  =RlUi(n,n) 


Therefore,  the  two-channel  solution  has  access  to  the  cross-reference  channel  com¬ 
ponents  R,il2  and  7?^  and  has  the  option  of  removing  them  from  either  reference 
channel  if  they  are  a  major  source  of  output  error.  This  fact  can  be  useful  in  practice 
as  it  can  potentially  aleviate  the  effect  of  additive  noise  in  the  reference  channels. 
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described  in  Section  2.3. 

In  order  to  confirm  this,  two  reference  channels  were  created  each  consisting 
of  a  single  complex  sinusoid  in  white  noise  "stepped"  in  frequency  at  iteration  point 
1536  to  match  the  primary  input  signals.  The  primary  channel  of  simulation  1  was 
also  used  here.  A  large  white  noise  component  (  SAT?.  -  lOrfi  )  was  added  to  both 
reference  channels  to  play  the  role  of  the  cross-channel  interference.  The  simulation  is 
fully  described  in  Table  2.3. 


Simulation  3 

before  step 

after  step 

main 

refl 

ref2 

main 

refl 

ref2 

W1 

JL 

8 

T 

4 

8 

00 

SNRi(db) 

20 

10 

20 

10 

6tt 

8 

8 

(hr 

8 

67T 

8 

SNR^db) 

20 

10 

20 

10 

SNRw(db) 

10 

r  10 

10 

10 

Table  2.3  Summary  of  Simulation  3. 


The  theoretical  evidence  (Equation  5a)  suggests  that  the  presence  of  this 
large  white  noise  component  will  render  noise-cancelling  impossible  on  a  single  channel 
basis.  The  multichannel  structure,  on  the  other  hand  should  be  able  to  compensate  for 
such  cross-channel  noise  and  proceed  with  its  noise-cancelling  task.  The  simulation 
results  provide  support  for  this  assertion.  The  single-channel  LSL  with  reference  chan¬ 
nel  1  is  unable  to  perform  noise-cancelling  (Figure  2.10a, b)  and  therefore  cascade  pro¬ 
cessing  of  one  channel  at  a  time  is  not  possible  in  this  setting.  On  the  other  hand,  the 
vector  (m  =2)  LSL  is  successful  in  eliminating  both  the  cross-channel  white  noise  com¬ 
ponent  and  the  sinusoidal  interference  from  the  primary  input  (Figure  2.10c,d,e). 


59 


Remarks 

This  scenario  will  arise  when  multiple  preformed  reference  beams  are  avail¬ 
able.  When  the  beams  overlap,  any  component(s)  entering  through  the  overlap  region 
(more  specifically  the  part  of  it  which  is  does  not  intersect  the  main  beam)  will  serve 
as  the  cross-channel  interference  components  and  will  have  a  deleterious  effect  on  can¬ 
cellation  via  cascade  processing. 

2.8  Chapter  summary 

The  adaptive  filtering  problem  was  set  in  the  proper  theoretical  framework. 
Adaptive  beamforming  was  presented  in  the  context  of  ANC.  The  origin  of  the  lattice 
structure  was  examined  and  the  relative  advantages  of  the  various  implementations 
were  considered.  The  basic  differences  between  the  stochastic  approximation  and  the 
exact  deterministic  least-squares  methods  was  established  Through  computer  simula¬ 
tions,  the  lattice  filters  were  shown  to  be  superior  to  direct  TDL  forms  in  a  noise¬ 
cancelling  setting.  The  deterministic  least-squares  lattice  was  shown  to  remain  stable 
under  drastic  magnitude  transients  that  drive  the  gradient  algorithms  to  instability. 
Finally ,  the  ability  of  the  multichannel  structure  to  cancel  components  correlated 
between  reference  channels  and  thus,  its  superior  performance  over  a  cascade  of 
single-channel  filters,  was  demonstrated. 
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Figure  2.1  The  noise  cancelling  concept,  (a)  General  structure,  (b)  Widrow’s  linear  model. 

Here  m,  and  mt  represent  uncorrelated  additive  noise  components  in  the  two 
channels.  The  H  mismatch"  between  the  primary  and  reference  channels, 
expressed  by  H(z),  must  be  compensated  for  by  the  adaptive  filter  transfer  func¬ 
tion  A(z). 
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Figure  2.2  Adaptive  beamforming  via  ANC.  (a)  Fixed,  preformed  reference  beam.  An  external 
steering  mechanism  is  required  to  track  the  changing  directional  properties  of  the 
interference,  (b)  Constrained  adaptive  beamformer.  The  multichannel  adaptive 
filter  weights  the  constrained  elements  to  form  a  spatially  adaptive  reference 
beam. 
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Figure  2.3  The  lattice  filter,  (a)  Prediction  error  (whitening)  filter;  k‘  and  k[  are  the  forward 
and  backward  reflection  (PARCOR)  coefficients;  e,(n)  and  r,(n)  are  the  order 
forward  and  beckward  prediction  error  sequences,  (b)  Joint-process  filter;  the 
r,(n)  constitute  an  orthogonal  basis  for  the  reference  process  x(n)  and  are 
weighted  by  the  cross-channel  coefficients  kf  to  produce  an  estimate  of  the  pri¬ 
mary  process  d(n). 
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Figure  2.4  Transversal  tapped-dclay-ltne  (TDL)  filter .  (a)  joint-process  configuration,  (b) 
one-step  forward  predictor. 


Figure  2.5  Simulation  1:  Eigenvalue  spread,  (a)  Primary  channel  input,  (b)  reference  channel  input,  (c)  LMS  output,  (d) 
GRL  output,  (e)  LSL  output.  Not  that  the  LMS  lags  significantly  in  cancelling  the  weak  sinusoid. 
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Figure  2. 


simulation  2i  scale  chonge 
LMS:  joint  process  error  output 
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simulation  2i  scale  change 
GRL:  joint  process  error  output 
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fl  Simulation  2;  Scale  change  (a)  LMS  output  (10,20,30  dB  jumps)  (b)  GRL  output  (10,20,30  dB  jumps).  Both 
structures  become  increasingly  unstable  as  the  severity  of  the  scale  change  increases.  (The  observed  clipping 
was  imposed  to  facilitate  display). 
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L5L:  joint  process  error  output 
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Figure  2.7  *  Se«/e  cAanje.  (a)  LSL  output  (10,10,30  dB  jumps)  (b)  LSL  "  likelihood  parameter"  (7)  (10  20  30  dB 
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GRL:  cross -channel  coefficient 
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LSL:  cross -channel  coefficient 
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Figure  2.8  Simti/afion  2:  Scale  change .  (a)  GRL  cross-channel  coefficient  (if)  and  (b)  LSL  cross-channel  coefficient  (if). 

The  instability  induced  by  the  power  transients  is  evident  in  the  GRL  coefficient,  especially  for  the  30  dB  jump. 

The  LSL  coefficient  exhibits  no  signs  of  instability. 
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Figure  2.9  Multichannel  vs  cascade,  (a)  Cascade  of  scalar  filters,  (b)  multichannel  joint-process  filter. 
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Figure  2.10  Simulation  8:  Multichannel  va  cascade .  (a)  Single  channel  LSL  joint-process  error  output  and  (b)  crosa-channel 
coefficient  (k*).  The  scalar  filter  is  unable  to  compensate  for  the  presence  of  a  large  uncorrelated  noise  com¬ 
ponent  in  the  reference  channel.  On  the  other  hand,  the  two-channel  LSL  joint-process  error  output  (c)  indi¬ 
cates  that  the  multichannel  structure  can  remove  the  additive  noise  component  common  to  both  reference 
channels.  This  is  confirmed  by  the  behavior  of  the  two  cross-channel  coefficients  (d,e). 
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CHAPTER  m 


BOUNDARY  REVERBERATION  REJECTION 
IN  A  SHALLOW  WATER  ENVIRONMENT 

In  this  chapter,  the  adaptive  beamforming  concept  will  be  applied  to  real 
reverberation  data  from  a  shallow-water  active  sonar  experiment.  The  shallow-water 
environment  creates  some  unique  difficulties.  As  opposed  to  deep-water  settings  where 
the  acoustic  paths  are  few  and  identifiable,  shallow-water  reverberation  returns  arrive 
in  rapid  succession,  imposing  stringent  adaptation  requirements  on  the  algorithm  and 
complicating  the  performance  evaluation  process. 

Two  types  of  problems  will  be  considered  here.  The  first  is  the  detection  of 
acoustic  signals  embedded  in  a  reverberation  backround.  This  is  the  problem  which 
may  be  encountered,  for  instance,  in  a  fisheries  sonar  application.  The  performance  of 
a  horizontally  directed  sonar  deployed  from  a  surface  vessel  is  often  reverberation- 
limited.  Particularly  difficult  to  detect  would  be  low-Doppler,  distributed  reflectors 
such  as  fish  schools  moving  slowly  and/or  at  an  angle  to  the  sonar  beam.  It  will  be 
shown  that  the  detectability  of  such  echo  returns  can  be  greatly  enhanced  through 
adaptive  filtering.  The  second  problem  is  selective  reverberation  cancellation  (SRC).  In 
this  case,  signal  and  noise  are  constituent  components  of  the  received  composite  rever¬ 
beration  process.  Volume  reverberation  is  normally  assumed  to  be  the  component  of 
interest  with  boundary  reverberation  being  the  interference. 

In  Chapter  2  we  have  shown  that  the  correct  choice  of  adaptive  processor  is 
of  critical  importance  for  succesful  reverberation  cancellation.  However,  this  choice  is 
not  the  only  factor  that  merits  consideration.  Equally  important  is  the  decision  con¬ 
cerning  the  filter  configuration  to  be  adopted  in  a  particular  experiment.  In  general, 
this  decision  .will  be  influenced  by  the  experimental  geometry,  the  degree  of  control 
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over  the  spatial  characteristics  of  the  system  which  is  afforded  by  the  data  acquisition 
system,  the  properties  of  "signal"  and  "interference"  and  the  objectives  of  a  particu¬ 
lar  deployment.  The  most  appropriate  structure  for  our  purposes  is  a  spatially  con¬ 
strained  adaptive  beamformer  based  on  the  adaptive  noise-cancelling  (ANC)  principle. 
Given  the  distinctly  different  scattering  mechanisms  responsible  for  the  three  rever¬ 
beration  types,  the  independence  assumptions  (see  Section  2.3)  are  readily  satisfied.  No 
further  assumptions  need  to  be  made  about  the  relative  intensity  and  duration  of  sig¬ 
nal  and  interference  and  no  modifications  of  the  algorithm  are  necessary  in  order  to 
impose  the  constraints.  The  question  remains  as  to  the  exact  manner  in  which  spatial 
constraints  are  imposed.  Spatial  prefiltering  can  range  from  complete  conventional 
beamforming  designed  to  encompass  the  interference  and  exclude  the  "signal"  com¬ 
ponent  from  the  reference  channel  to  simple  element-to-element  subtraction  followed 
by  adaptive  beamforming  controlled  by  a  multichannel  algorithm.  Both  fixed,  pre¬ 
formed  and  spatially  constrained  adaptive  reference  beams  will  be  used  in  this  chapter 
with  attention  to  their  relative  advantages.  The  possibility  of  applying  the  results  on 
the  (vertical)  spatial  correlation  characteristics  of  boundary  and  volume  reverberation 
presented  in  Section  1.4,  in  conjunction  with  direct  ANC,  will  also  be  considered.  The 
Least-Squares  Lattice  (LSL)  in  its  scalar  and  multichannel  forms  will  be  exclusively 

used  as  the  adaptive  processor,  based  on  theoretical  considerations  and  simulation 
results  presented  in  Chapter  2. 

3.1  The  Dabob  Bay  environment 

The  experiment  was  conducted  in  the  deepest  part  of  Dabob  Bay,  Washing¬ 
ton  (maximum  depth  ~200m).  The  physical,  biological  and  acoustical  properties  of 
Dabob  Bay  have  been  studied  intensively  [Helton,  1976].  It  is  characterized  by  large 
fluctuations  of  temperature  and  salinity.  In  addition  to  the  usual  seasonal  variations 
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in  the  top  ~30 m,  cellular  inhomogeneities  exist  in  the  deeper  regions,  caused  by  intru¬ 
sions  of  Pacific  Ocean  water.  Water  parcels  have  been  detected  with  distinctly 
different  values  of  temperature  and  salinity,  with  thickness  of  ~40  m  and  horizontal 
extent  up  to  ~1000  m  [Ebbesmeyer,  1973),  which  may  remain  intact  for  several  weeks 
before  they  dissipate.  These  transient  water  masses  cause  wide  fluctuations  in  the 
sound  velocity  structure  and  create  the  need  for  up-todate  sound  velocity  profiling 
during  acoustic  experiments.  Coring  device  sampling  [Burns,  1962]  has  shown  that  the 
deep  regions  of  Dabob  Bay  are  covered  by  silt  and  mud  with  some  sand  and  fine 
gravel.  A  harder,  highly  reflective,  layer  consisting  of  glacial  till  (a  dense  mixture  of 
sand,  gravel  and  silt)  begins  1— 15m  below  the  mudline.  A  maximum  bottom  back- 
scattering  coefficient  of  —20 db  was  measured  at  ZOkHz.  Correspondingly,  the  surface 
backscattering  coefficient  was  found  to  be  in  the  range  of  -35  to  -25</6.  Intensive 
volume  sampling  and  volume  scattering  measurements  over  a  period  of  two  years  has 
revealed  multiple  scattering  layers  with  substantial  seasonal  variations  [Anderson, 
1981).  The  maximum  measured  volume  scattering  coefficient  is  ~-60<f6. 

3.1.1  Experiment  Description 

The  sonar  system  moved  with  constant  velocity  along  straight  line  paths.  It 
carried  a  transducer  array  whose  elements  were  combined  by  a  programmable  conven¬ 
tional  beamformer  capable  of  producing  several  transmit/receive  beam  sets.  The 
overall  system  geometry  is  depicted  in  Figure  3.1.  The  transmit  beam  pattern  choice 
included  a  "wide"  (60°  X  60  °)  and  a  "narrow"  (15°  X  15°)  beam.  Two  receive 
beam  sets  were  available.  The  first  (beam  set  0)  consisted  of  9  preformed  (15  °  X  15  ° ) 
beams  steered  in  the  directions  indicated  in  Table  3.1.  In  this  coordinate  system 
elevation  runs  from  -90°  (up)  to  90°  (down)  and  bearing  from  -180°  to  180  °  relative 
to  broadside  (0  ° ). 


73 


Beam  number 

1 

2 

3 

4 

5 

6 

7 

8 

9 

Elevation 

o- 

-7.5  * 

-7.5* 

-7.5* 

-7.5* 

7.5* 

7.5* 

7.5* 

7.5* 

Bearing 

0* 

15- 

7.5* 

-7.5* 

-15* 

15* 

7.5* 

-7.5* 

-15’ 

Table  3.1  Steering  directions  of  preformed  beam  set 


It  may  be  noticed  that  beams  3, 4, 7, 8  overlap  partially  with  the  center  beam  (1).  This 
beam  arrangement  has  good  potential  as  a  noise-cancelling  configuration.  The  second 
receive  set  (beam  set  1)  contains  the  outputs  of  9  row  sums  comprising  the  array  face 
and  constitutes  a  relatively  "raw"  data  set  well-suited  for  multichannel  adaptive 
beamforming  operations.  The  sonar  system  transmitted  180m«  pulses  and  recorded 
digitally  the  complex-basebanded  reverberation  return.  Narrowband  calibration 
returns  arriving  broadside  to  the  array  are  occasionally  present  which  are  useful  "con¬ 
trol”  signals  valuable  in  assessing  the  cancellation  results.  CTD  profiles  were  recorded 
before  each  experimental  run. 

3.2  Adaptive  "  whitening"  of  reverberation  data 

Before  we  embark  on  dual-channel  ANC  operations  it  would  be  of  interest  to 
examine  the  effect  of  single-channel  processing  with  the  prediction-error  (whitening) 
filter  on  real  reverberation  data.  Two  potential  benefits  come  to  mind.  First,  time 
evolving  "high  resolution"  spectral  estimates  can  be  obtained  or,  equivalently,  the 
appropriate  AR  model  for  reverberation  can  be  identified.  Second,  by  adjusting  the 
adaptation  coefficient,  the  sensitivity  of  the  filter  to  transient  signals  off  the  main 
reverberation  frequency  "ridge"  can  be  substantially  reduced.  As  a  result,  these  sig¬ 
nals  may  be  allowed  to  "leak"  in  the  error  output  with  considerably  improved  SNR. 

The  composite  reverberation  spectrum  has  a  complex  and  highly  variable 
structure.  Boundary  reverberation  has,  typically,  an  abrupt  onset.  Morover,  surface 
reverberation  is  often  characterized  by  Doppler  shifts  distinctly  different  from  the 
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volume  return  [Igarashi  and  Stern,  1970;  Swarts,  1971].  It  follows  that  effective  whi¬ 
tening  of  the  reverberation  ridge  cannot  be  accomplished  with  a  fixed  notch  filter.  On 
the  other  hand,  the  frequency  tracking  capability  of  the  adaptive  filter  should  allow  it 
to  be  more  successful  to  that  end. 

3.2.1  Data  description 

Ping  EXP3.1  was  chosen  to  undergo  adaptive  "prewhitening".  During  this 
run,  the  sonar  depth  was  100  m  and  receive  beam  set  0  (9  preformed  beams)  was  in 
effect.  Only  the  output  of  the  center  beam  (beam  1  in  Figure  3.1b)  is  considered  here. 
The  sound  velocity  profile  measured  prior  to  this  run  and  the  associated  ray-tracing 
diagram  are  displayed  in  Figure  3.2.  The  received  reverberation  is  highly  composite 
with  both  boundaries  contributing  to  the  return.  Because  no  real  signals  are  present  in 
this  data  set,  two  synthetic  echoes  (Tl  and  T2)  were  inserted  in  the  vicinity  of  the 
reverberation  ridge  in  order  to  demonstrate  the  degree  of  signal  enhancement  attained 
through  adaptive  whitening.  They  are  described  in  Table  3.2  below. 


Signal  echo 

Tl 

T2 

Onset  time 

1.5  s 

2.5  s 

Duration 

180  ms 

180  ms 

SNR 

—  0 

-0 

Velocity 

0  m/s 

0.375  m/s 

Table  3.2  Sythetic  echo  description 


The  input  waveform  is  displayed  in  Figure  3.3.  The  time  series  plot  was  made  follow¬ 
ing  low-pass  filtering  and  decimation  by  a  factor  of  8  in  order  to  facilitate  the  display 
of  the  time-varying  features  of  the  data.  The  Range  Doppler  Map  (RDM)  plot  is 
intended  to  show  the  time-varying  character  of  the  ping  in  the  frequency  domain.  It 
was  constructed  via  successive  128-point  FFTs  of  windowed  (Kaiser-Bessel,  <*=2.5) 
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data  segments  with  50°$  overlap.  The  (dB)  magnitude  of  each  transform  corresponds 
to  a  single  line  in  the  pseudo-3d  plot.  The  Time-Varying-Gain  (TVG)  applied  to  this 
ping  was  designed  to  compensate  for  the  attenuation  factor  due  to  spherical  spreading 
The  fact  the  cylindrical  spreading  (~r“*)  is  in  effect  for  most  of  the  return  in 
the  shallow  bay  explains  the  observed  rising  level  with  range.  The  main  reverberation 
ridge  occupies  the  near-DC  region,  following  complex  basebanding.  The  secondary 
ridge  at  un  =0.25  is  not  physically  meaningful  and  could  possibly  be  the  result  of  alias¬ 
ing  suffered  during  the  data  acquisition  phase.  The  synthetic  signal  echoes  are  well- 
concealed  by  reverberation  in  both  the  time  and  frequency  domain. 

3.2.2  Processing  results 

The  AR-model  order  for  this  ping  was  determined  by  making  use  of  the  AR- 
cutoff  property  of  the  lattice  filter  (see  section  2.4.1).  The  reflection  coefficients  of  ord- 
ers  (>  7)  were  found  to  be  virtually  equal  to  zero  and  therefore  p  =  6  was  accepted  as 
the  appropriate  filter  order.  This  is  consistent  with  the  findings  of  Haykin  et.  al.[1982] 
in  their  investigation  of  radar  clutter,  the  atmospheric  counterpart  to  reverberation 
and  of  Hansen  [1984]  in  his  study  of  the  oceanic  reverberation  spectrum.  Both  have 
concluded  that  reverberation  (clutter)  can  be  effectively  modelled  as  an  AR  process  of 
relatively  low  order. 

Three  whitening  runs  were  made  ,  each  with  the  scalar  LSL  prediction-error 
filter  of  order  p  =6  and  adaptation  coefficients  (aWL)  of  0.0002,  0.002  and  0.02  respec¬ 
tively.  The  results,  displayed  in  Figure  3.4  include  the  time  evolving  "high  resolution" 
spectral  estimate  and  the  RDM  of  the  prediction  error  output  for  these  values  of 
For  aLSL  =0.0002,  the  spectral  estimate  fails  to  resolve  T2  which,  as  a  result,  is  not 
whitened  and  is  preserved  at  the  output  with  a  much  improved  SNR.  For  ocLSL  =0.002, 
T2  is  partially  resolved  and  correspondingly  partially  removed  from  the  output.  For 
aLSL  =  0*02,  the  spectral  estimate  includes  many  of  the  details  of  the  original  RDM 
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and  to  a  large  extent,  tracks  its  time-varying  behavior.  Signal  T2  is  completely 
resolved  and  compensated  for.  The  zero-Doppler  signal  (Tl)  was  cancelled  every  time, 
which  is  to  be  expected  since  it  resides  in  the  frequency  null  induced  by  the  reverbera¬ 
tion  process. 

Remarks 

A  clear  trade-off  is  evident  regarding  the  choice  of  adaptation  coefficient.  If 
the  objective  is  enhanced  detectability  of  a  transient  low-Doppler  echo  in  the  vicinity 
of  the  reverberation  ridge,  then  a  relatively  small  is  in  order.  In  this  case,  the 
effective  averaging  may  be  excessive  for  the  spectral  modeling  objective  and  may 
result  in  the  "smearing”  of  certain  time-varying  spectral  features.  A  larger  af gfr.  on 
the  other  hand,  will  result  in  a  more  realistic  spectral  estimate  at  the  expense  of  com¬ 
plete  elimination  of  all  signals  from  the  output.  It  should  be  mentioned  that  signal 
duration  is  also  important  with  regard  to  this  trade-off.  Gradual,  persistent  echoes 
from  spatially  distributed  scatterers  are  more  likely  to  be  whitened  than  transient 
returns  from  discrete  reflectors.  Morover,  increased  signal  cancellation  will  result  as 
the  proximity  of  the  echo  to  the  reverberation  ridge  is  increased.  In  view  of  these  con¬ 
cerns,  the  adaptive  whitening  operation  is  of  limited  utility  in  the  context  of  echo 
detection. 

3.3  ANC  with  fixed,  preformed  reference  beams 

A  more  promising  alternative,  for  signal  detection  in  a  reverberation  back- 
round,  is  dual-channel  ANC  processing  with  preformed,  fixed  reference  beams  partially 
overlapping  with  the  main  beam.  Signal  components  common  between  the  main  and 
the  reference  beams  will  be  subject  to  cancellation.  On  the  other  hand,  echoes  arriving 
through  the  angular  region  covered  exclusively  by  the  main  beam,  will  be  protected 
from  cancellation.  It  follows  that  such  echo  returns  will  be  enhanced  by  an  ANC 
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operation  between  the  main  and  the  reference  beam.  In  order  to  demonstrate  this, 
the  same  ping  (EXP3.1)  will  be  processed  with  the  multichannel,  joint-process  LSL 
filter,  the  (multichannel)  reference  input  consisting  of  the  four  interior  beams  of  the 
multibeam  set. 

The  joint-process  error  output  for  p  =3  and  aLSL  =0.02  is  shown  in  Figure 
3.5.  The  reverberation  backround  is  reduced  substantially,  rendering  signals  Tl  and 
T2  clearly  visible  in  both  the  time  and  frequency  domains.  Also  displayed  is  a  dB- 
cancellation  curve  which  is  the  low-pass  filtered  ratio  of  the  input  signal  magnitude 
over  the  output  error  magnitude  and  shows  the  overall  character  of  cancellation 
achieved.  It  is  noteworthy  that  the  zero-Doppler  echo  (Tl)  is  equally  well  preserved  as 
the  other  low-Doppler  echo  (T2).  In  addition,  signal  duration  is  no  longer  an  important 
factor  for  signal  preservation.  The  sole  requirement  is  that  the  echo  arrival  angles  fall 
within  the  protected  angular  region  of  the  main  beam.  Therefore,  even  persistent, 
zero-Doppler  echoes  from  distributed  scatterers  can  potentially  be  extracted  through 
this  process. 

Remarks 

1*  It  should  be  mentioned  that  in  this  case  the  signal  echoes  are  completely 

absent  from  all  reference  beams.  Should  a  signal  component  be  allowed  to 
enter  one  of  the  reference  channels,  ANC  theory  predicts  potentially  serious 
signal  cancellation  at  the  output.  Therefore,  in  designing  reference  beams, 
emphasis  should  be  placed  in  suppressing  their  sidelobe  response  in  the  direc¬ 
tions  of  high  main  beam  response. 

ii.  This  signal  processing  technique  is  likely  to  be  useful  in  a  horizontal  fisheries 

sonar  application.  The  Doppler  structure  of  a  fish  school  is  determined  by 
varying  degrees  of  systematic  and  random  motion,  including  rhythmical 
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swimming  motions,  cruising,  short-term  accelerations  and  erratic  changes  in 
speed  and  direction  [Holliday,  1974).  The  aspect  between  the  swimming  direc¬ 
tion  of  the  school  and  the  main  beam  axis  is  also  a  factor.  It  follows  that  the 
received  waveform  will  be  characterized  by  varying  degrees  of  net  Doppler 
shift  and  with  the  exception  of  fish  swimming  at  high  speeds  along  the  beam 
axis,  fish  detection  with  a  Doppler  sonar  is  reverberation  limited.  Reverbera¬ 
tion  interference  is  particularly  severe  in  shallow-water  applications  such  as 
the  tracking  of  migrating  species  in  rivers  [Pincock  and  Easton,  1978;  Hen- 
dershot  and  Acker,  1984].  The  results  of  this  section  indicate  that  ANC  with 

appropriately  formed  reference  beams  offers  the  potential  for  improved  fish- 

« 

count  estimation  in  such  reverberation-limited  environments. 

m.  The  question  remains  as  to  the  composition  of  the  residual  reverberation  sig¬ 
nal.  It  seems  unlikely  that  the  particular  beam  set  used  is  adequate  for  selec¬ 
tive  reverberation  cancellation.  Since  the  reference  beams  do  not  encompass 
all  possible  directions  of  arrival  of  boundary  reverberation  components  enter¬ 
ing  the  mam  beam,  one  may  expect  the  reverberation  present  at  the  filter 
output  to  be  still  highly  composite.  The  following  REVGEN  simulation  was 
intended  to  resolve  this  question 

3.3.1  REVGEN  simulation  (RS3.1) 

Two  simulation  runs  were  performed.  In  the  first  run,  boundary  as  well  as 
volume  scatterers  were  placed  in  the  medium.  In  the  second  run,  only  volume  scatter- 
era  were  included  in  order  to  create  a  "desired"  signal  in  the  SRC  scenario.  System 
parameters  were  duplicated  to  the  extent  possible.  Simplified  "pencil"  beams  were 
used  and  an  isovelocity  medium  was  assumed.  A  summary  of  the  simulation  parame¬ 
ters  is  given  in  Table  3.3  below. 
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Simulation  parameter 

Symbol 

Value 

Pulse  length 

r 

180ms 

Shell  thickness 

T 

15  m 

No.  of  volume  scatterers  per  slice 

N* 

1  ’»i 

100 

No.  of  surface  scatterers  per  slice 

K 

50 

No.  of  bottom  scatterers  per  slice 

K 

50 

Layer  height 

H 

175 

Sound  velocity 

c 

1500  m/s 

Volume  reverberation  coefficient 

-70  dB 

Surface  reverberation  coefficient 

-30  dB 

Bottom  reverberation  coefficient 

*b 

-20  dB 

Table  3.3  REVGEN  simulation  RS3.1 


The  approximate  number  of  scatterers  per  scattering  volume  and  area  for  volume  and 
surface  (bottom)  reverberation  respectively,  are  given  by 


and 


N., 


N'„0ct 
4  kT 


n;jct 

4jiT 


where  0  is  the  equivalent  horizontal  receive  beam  width  in  radians.  The  above  com¬ 
putations  are  based  on  a  wedge  approximation  for  the  receive  beam  patterns  which 
is  valid  at  the  far  ranges,  following  beam  intersection  with  the  boundaries.  Here 
0  =  15  *  =  0.262rai.  Substituting,  we  have 


and 


N„  =  0.375 NJ,  =  37.5 


N,f  =0.375 N;,  =18.75 

Therefore,  a  sufficiently  large  number  of  scatterers  contributed  for  the  reverberation 
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return  to  "converge"  to  a  Gaussian  distribution. 

No  signal  echoes  were  included  in  this  simulation.  The  composite  reverbera¬ 
tion  data  were  processed  identically  to  the  real  case  and  the  results  are  displayed  in 
Figure  3.6.  It  can  be  seen  that  the  ANC  output  (dotted  line)  falls  between  the  input 
(composite)  reverberation  process  and  the  "desired"  volume  reverberation  signal. 
Therefore,  the  output  is  indeed  of  composite  nature.  The  dB-cancellation  curve  shows 
similar  overall  cancellation  as  in  the  real  data  case. 

3.4  An  Application  of  Constrained  Adaptive  Beamforming 

In  this  section  we  focus  on  SRC.  Volume  reverberation  is  the  "signal"  of 
interest.  The  "interference"  consists  of  surface  and  bottom  reverberation  entering  ini¬ 
tially  through  the  sidelobes  and  eventually  through  the  main  lobe  of  the  receive  beam. 
We  will  examine  the  possibility  of  recovering  the  volume  component,  normally  weak 
and  completely  masked  by  the  boundary  return,  through  adaptive  processing. 

The  feasibility  of  ANC  with  fixed,  preformed  beams  depends  on  the  sidelobe 
structure  of  the  particular  transmit  and  receive  beam  patterns.  Two  different  refer¬ 
ence  beam  sets  are  in  general  required  for  side  lobe  and  main  lobe  interference  cancel¬ 
lation.  In  constructing  a  side  lobe  cancellation  beam,  great  care  must  be  exercised  to 
closely  match  the  side  lobe  features  of  the  main  beam.  This  may  prove  difficult  for 
beams  with  complex  side  lobe  structure.  The  main  lobe  cancellation  beam,  on  the 
other  hand,  must  have  a  null  in  the  direction  of  interest  in  order  to  protect  the 
volume  signal  from  cancellation.  As  the  interference  arrival  angle  changes  with  range, 
the  fixed  null  shape  may  be  a  serious  limitation  to  cancellation  in  the  far  ranges.  Mor- 
over,  an  external  decision  is  needed  on  the  appropriate  time  for  switching  from  side 
lobe  cancellation  to  main  lobe  cancellation  mode.  The  inevitable  transient  response  of 
the  algorithm  to  such  a  discontinuous  "step"  presents  an  additional  problem. 
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As  an  alternative,  one  may  allow  the  adaptive  algorithm  to  directly  control 
the  spatial  transfer  function  of  the  system,  provided  that  spatially  separated  receive 
elements  are  independently  available  and  to,  in  effect,  create  its  own  adaptive  refer¬ 
ence  (cancellation)  beam.  With  appropriately  constrained  reference  elements,  this 
configuration  has  the  potential  of  cancelling  both  side  lobe  and  main  lobe  interference 
with  a  graceful  transition  between  the  two  modes.  In  addition,  it  enjoys  a  significant 
advantage  in  operational  flexibility  as  it  needs  little  or  no  a  priori  information  about 
the  shape  of  the  main  receive  beam.  For  these  reasons  this  Constrained  Adaptive 
Beamforming  operation  is  especially  promising  for  SRC. 

3.4.1  Data  description 


The  experimental  data  set  offers  an  opportunity  for  application  of  this  con¬ 
cept.  Ping  EXP3.2  is  a  good  case  in  point.  The  same  Mwide"  transmit  beam  set  is  in 
effect  and  beam  set  1  was  used  to  receive.  The  latter  consists  of  the  individually 
recorded  outputs  of  9  array  row  sums.  These  rows  will  be  used  as  the  reference  ele¬ 
ments  in  the  adaptive  beamforming  process.  The  basic  experimental  geometry  is  unal¬ 
tered  with  the  exception  of  sonar  depth  which  now  is  30  m.  The  velocity  profile  and 
ray  tracing  diagram  for  this  run  are  displayed  in  Figure  3.7. 

A  REVGEN  simulation  (RS3.2)  was  carried  out  for  this  scenario.  The  simula¬ 
tion  parameters  of  Table  3.1  were  used  and  the  receive  elements  were  approximated 


by  rectangular  pistons  at 


\ 

2 


spacing. 


The  main  beam  was  constructed  by  adding  the  outputs  of  the  9  rows.  The 
simulated  beam  pattern  is  shown  in  Figure  3.8a.  In  addition,  8  constrained  reference 
elements  were  created  by  pairwise  subtraction  among  the  9  rows,  in  the  manner  dep¬ 
icted  in  Figure  2.2b.  The  beam  pattern  of  a  single  constrained  element,  characterized 
by  a  central  null,  is  the  dotted  curve  in  Figure  3.8b. 
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The  real  and  simulated  (composite)  reverberation  signal  entering  the  main 
beam  are  displayed  in  parallel  in  Figure  3.9.  The  first  plot  is  the  regular  RDM  display. 
The  second  plot,  an  RDM  with  ~  90%  overlap,  provides  a  more  detailed  view  of  the 
near  ranges  where  boundary  reverberation  enters  through  the  sidelobes.  The  data  has 
been  frequency-shifted  to  compensate  for  the  Doppler  shift  due  to  sonar  motion. 
Clearly  visible  is  the  typical  "  hook"  pattern  of  negative  Doppler  features  in  the  initial 
~  0,5  3  of  the  pinS-  This  Pattern  is  entirely  due  to  high-angle  boundary  returns  which 
are  imparted  a  lesser  net  Doppler  shift  by  the  velocity  of  the  sonar.  With  increasing 
range,  as  the  lower  sidelobes  and  eventually  the  main  lobe  become  the  boundary 
reverberation  outlets,  the  pattern  merges  into  the  main  ridge.  This  artifact  is  useful 
for  identifying  sidelobe  returns  and  establishing  their  cancellation  following  adaptive 
processing.  The  final  plot  is  the  regular  time-series  diplay. 

Although  the  real  and  simulated  data  have  similar  overall  structure,  several 
differences  are  apparent.  The  real  data  are  characterized  by  a  somewhat  wider 
"hook"  which  may  be  an  indication  that  the  real  beam  pattern  has  an  appreciable 
side  lobe  response  at  a  higher  angle  than  the  simulated  beam  or  that  the  actual  velo 
city  of  the  sonar  was  underestimated  in  the  simulation.  In  addition,  the  simulated 
reverberation  is  better  sustained  at  long  ranges.  Two  reasons  may  be  cited  for  this. 
First,  no  decaying  function  of  grazing  angle  was  introduced  in  the  simulation  for  boun¬ 
dary  reverberation.  This  decay  with  grazing  angle  is  a  well-documented  fact  [Urick, 
1975).  The  second  reason  may  be  inferred  from  the  raytracing  diagram  (Figure  3.7) 
and  the  time  of  arrival  vs  angle  of  arrival  plot  (Figure  3.10).  Because  of  the  well- 
developed  thermocline,  a  shadow  zone  is  formed  in  the  near  surface  past  a  range  of  ~ 
500  m.  The  surface  contribution  disappears  abruptly  at  that  point.  In  the  simulated 
case,  an  isovelocity  profile  is  used  and  the  surface  continues  to  contribute  at  long 
ranges.  Finally,  a  high-Doppler  calibration  signal  can  be  observed  in  the  real  data 
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RDM  at  It  arrives  broadside  to  the  array.  This  signal  was  not  included  in  the 

simulation. 

3.4.2  Processing  with  a  single  constrained  reference  element 

A  single  constrained  reference  element,  created  by  subtracting  one  row  ele¬ 
ment  from  another,  has  a  beam-pattern  (Figure  3.8b)  which  at  first  glance  appears  to 
be  well-suited  for  SRC.  Excluding  its  broadside  null,  it  is  nearly  omnidirectional  in  the 
vertical,  thus  covering  all  possible  directions  of  arrival  for  boundary  reverberation. 
One  may  expect  it  to  be  effective  in  cancelling  the  sidelobe-born  boundary  returns  and 
to  some  extent  main  lobe  interference,  while  protecting  much  of  the  volume  signal  in 
its  null. 

Ping  EXP3.2  was  processed  in  this  manner,  using  the  scalar  LSL  with  p  =3 
and  at lsl  =0.02.  The  single  reference  channel  consisted  of  rows  (5-6).  The  results  for 
both  the  real  and  simulated  data  are  displayed  in  Figure  3.11.  Cancellation  is  gen¬ 
erally  low  in  the  near  ranges.  In  addition,  the  simulated  case  displays  a  significant 
increase  in  backround  noise  level  for  t  <  0.5*.  Moderate  cancellation  was  achieved  in 
both  cases  for  intermediate  ranges  (l.O  <  t  <2.0)  while  the  behavior  in  the  far  ranges 
is  distinctly  different  for  real  and  simulated  data.  For  the  latter,  cancellation  contin¬ 
ues  to  improve  with  range  while  little  or  no  cancellation  is  observed  for  the  former. 

The  low  overall  cancellation  in  the  near  ranges,  is  a  manifestation  of  the 
ANC  theoretical  prediction  regarding  the  effect  of  uncorrelated  noise  components  in 
either  the  primary  or  the  reference  channel.  This  effect,  which  was  outlined  in  Section 
2.2,  comes  into  play  here  as  follows:  The  transmit  beam  ensonifies  the  entire  angular 
region  from  -60*  to  450*.  The  reference  beam,  being  nearly  omnidirectional  in  the 
vertical,  receives  reverberation  from  the  entire  range  with  the  exception  of  the  near¬ 
broadside.  The  main  beam,  on  the  other  hand,  has  a  definite  side  lobe  response  which 
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operates  on  the  high  angle  returns.  Each  time  boundary  reverberation  arrives  through 
the  nulls  between  side  lobes  it  is  excluded  from  the  primary  channel.  The  same  com¬ 
ponents  enter  the  reference  channels  without  being  attenuated  and  in  effect  cause  the 
uncorrelated-to-correlated  noise  density  ratio  (M^z)  in  Eq.  2.5a  )  to  increase,  resulting 
in  reduced  signal-to-noise  ratio  at  the  filter  output.  This  explains  the  observed  rise  of 
the  noise  "floor".  The  somewhat  better  cancellation  apparent  for  the  high  angle 
(±60  )  returns  in  the  real  vs  the  simulated  data,  may  be  an  artifact  of  the  RDM 
displays  caused  by  dynamic  range  differences,  unavoidable  in  view  of  the  necessarily 
incomplete  model  of  physical  reality  which  is  used  in  the  simulation.  The  particular 
side  lobe  structure  of  the  real  main  beam,  which  may  be  substantially  different  from 
the  simulated  beam,  could  also  explain  this  difference. 

The  difference  in  far-range  cancellation  performance  is  due  to  the  grazing 
angle  effect.  In  the  real  case,  there  is  simply  little  boundary  reverberation  left  in  the 
far  ranges  and  therefore  the  cancellation  potential  is  low.  This  effect  was  not  included 
in  the  simulation  and  as  a  result  significant  cancellation  was  achieved  there.  Interest¬ 
ingly,  there  is  a  clear  trend  of  improving  cancellation  performance  with  range,  which 
is  consistent  with  another  ANC  theoretical  prediction,  this  one  regarding  the  noise 
spectrum  at  the  filter  output  (relation  2.5c).  Specifically,  the  output  noise  spectrum 
was  found  to  be  proportional  to  the  signal-to-  correlated  noise  ratio  at  the  primary 
input.  In  our  experimental  geometry,  this  ratio  should  decrease  with  range  as  the 
boundary  returns  begin  to  enter  through  the  main  lobe  instead  of  the  side  lobes. 
Therefore,  one  may  expect  the  output  noise  power  to  decrease  and  cancellation  to 
improve.  This  concept  may  be  tested  through  the  simulated  data.  Once  more,  two 
separate  simulation  runs  were  performed  with  only  volume  and  only  boundary  scatter- 
ers  present,  respectively.  Thus,  "signal"  and  "interference"  were  isolated.  We  have 
calculated  the  volume  to  boundary  reverberation  power  ratio  at  the  output  of  the 
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main  beam  and  plotted  it  together  with  the  dB-cancellation  curve  in  Figure  3.12b.  It 
is  seen  that  the  two  curves  are  nearly  mirror  images  of  each  other.  The  decreasing 
trend  in  the  main  channel  SNR  corresponds  well  to  the  increasing  cancellation  trend. 
In  addition,  many  of  the  large-scale  fluctuations  seem  to  obey  the  same  relationship. 
Figure  3.12a  clearly  shows  that,  although  cancellation  improves  with  range,  this 
configuration  is  inadequate  for  SRC. 

Two  additional  observations  can  be  made,  concerning  the  real  data  results. 
First,  the  high-Doppler  "control"  signal  was  virtually  unaffected  by  this  processing. 
This  is  an  indication  that  the  constraining  scheme  is  effective.  Second,  a  large  zero- 
Doppler  feature  was  revealed  (marked)  at  t~  1.0a.  This  is  probably  a  volume  rever¬ 
beration  feature  the  origin  of  which  will  be  speculated  upon  in  a  later  section. 

Remarks 

i.  Processing  with  a  single  constrained  reference  beam  would  be  more  effective 

in  the  side  lobe  cancellation  mode,  if  the  transmit  beam  were  identical  to  the 
main  receive  beam.  In  that  case,  the  "uncorrelated  noise"  effect  would  not 
exist  and  near-range  cancellation  would  improve. 

Regarding  the  constrained  reference  cancellation  beam,  a  choice  must  be 
made  about  the  width  of  thr  broadside  null.  It  would  be  desirable  to  have  a 
wide  null  initially,  when  interference  arrives  through  the  sidelobes  and  a  pro¬ 
gressively  narrower  null  as  the  main  beam  begins  to  intersect  the  boundaries. 

m.  An  alternative  explanation  for  the  improved  cancellation  observed  in  the  far 
ranges  is  in  terms  of  the  increased  spatial  coherence  of  boundary  reverbera¬ 
tion  predicted  by  the  point-scattering  model  (Section  1.4.) 
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3.4.3  Processing  with  a  spatially  adaptive  cancellation  beam 

The  results  of  Simulation  3  indicate  that  the  multichannel  joint-process  filter 
can  compensate  for  noise  components  not  present  in  the  primary  channel,  which  are 
correlated  between  the  reference  channels.  Therefore,  if  the  multichannel  LSL  filter 
were  to  be  used  in  the  present  case,  with  the  8  constrained  reference  elements 
comprising  the  vector  reference  input,  one  may  expect  it  to  be  more  effective  than  the 
single  channel  filter  by  aleviating  the  near  range  "uncorrelated  noise"  effect. 
Equivalently,  the  multichannel  adaptive  algorithm  should,  in  theory,  form  an  adaptive 
reference  (cancellation)  beam  with  the  tendency  to  emulate  the  side  lobe  features  of 
the  main  beam.  This  adaptive  reference  beam  will  still  have  the  central  null  according 
to  the  beam  pattern  product  theorem  [Urick,  1975]  and  therefore,  the  basic  constrain¬ 
ing  scheme  will  still  be  in  effect. 

The  results  for  the  same  ping  (EXP3.2)  processed  through  the  8-channel  LSL 
filter  are  displayed  in  Figure  3.13.  Overall,  cancellation  is  substantially  improved  over 
the  single  reference  element  case.  As  evidenced  by  the  virtual  elimination  of  the 
hook  pattern  from  the  RDM  plots,  the  side  lobe  boundary  interference  has  been 
successfully  removed.  Somewhat  surprising  is  the  dramatic  improvement  in  far  range 
cancellation.  As  a  result,  in  the  real  data  case,  volume  feature  Vl  is  even  more  clearly 
in  evidence  at  t  ~  1.0  s  and  an  additional  feature  (V2)  was  exposed  at  t  ~  1.6  s. 

Figure  3.14a  shows  that  the  multichannel  constrained  adaptive  beamformer 
is  generally  very  successful  towards  achieving  SRC.  This  is  particularly  true  for 
^  —  L0  6.  Insight  may  be  gained  on  the  adaptive  operation  performed  by  this  struc¬ 
ture,  by  considering  the  adaptive  cancellation  beam  it  creates.  Vertical  "cuts"  of  this 
beam,  at  zero  bearing,  were  calculated  from  each  set  of  8  cross-channel  reflection 
coefficients  at  ~0.024  a  intervals  and  are  displayed  in  a  "waterfall"  plot  (Figure 
3.14b).  Note  that  the  adaptive  cancellation  beam  still  has  the  central  adaptive  null, 
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as  expected.  Morover,  it  strives  to  approximate  the  sidelobe  structure  of  the  main 
beam  in  the  early  ranges  of  high-angle  boundary  returns.  At  longer  ranges  ((  >  1.0  «) 
when  the  side  lobe  activity  subsides,  the  adaptive  cancellation  beam  concetrates  on 
the  main  lobe  boundary  returns  to  the  extent  allowed  by  the  broadside  null. 

3.4.4  REVGEN  simulation:  symmetric  case 

In  order  to  futher  examine  the  nature  of  the  adaptive  cancellation  beam,  the 
same  processing  was  carried  out  in  a  more  benign  symmetric  experimental  geometry. 
In  this  REVGEN  simulation  (RS3.3),  the  overall  depth  was  set  to  300  m  and  the  sonar 
was  positioned  at  150  m.  The  surface  and  bottom  backscattering  coefficients  were 
both  set  equal  to  -30  db.  In  the  sequence  of  plots  presented  in  Figure  3.15,  it  is 
apparent  that  the  symmetrical  arrangement  is  beneficial  to  the  SRC  performance  in 
the  near  ranges.  The  adaptive  cancellation  beam  is  itself  symmetric  and  its  behavior 
easier  to  interpret.  In  its  side  lobe  cancellation  mode  (f  =0.5)  it  succesfully  emulates 
the  side  lobe  structure  of  the  main  beam.  As  the  main  lobe  begins  to  intersect  with 
the  boundaries,  the  adaptive  filter  is  seen  to  adaptively  alter  the  effective  width  of  the 
broadside  null,  making  it  progressively  narrower  with  increasing  range  (plots  for 
t  =1.0  and  t  =2.0  in  Figure  3.15c).  This  is  intuitively  satisfying  and  explains  the 
dramatic  improvement  in  far-range  cancellation  performance  achieved  by  the  adap¬ 
tive  beamformer  over  single  channel  (i.e.  fixed  null  shape)  ANC. 

3.4.5  Volume  M  feature"  interpretation 

A  number  of  hypotheses  were  considered  regarding  the  origin  of  the  volume 
features  (Vl  and  V2)  exposed  by  the  adaptive  beamforming  operation  on  the  real 
data.  First,  the  possibility  that  they  are  acoustic  echoes  from  schools  of  fish  or  cellu¬ 
lar  water  masses  within  the  angular  range  of  the  null  was  debated.  This  hypothesis 
was  rejected  based  on  the  processing  results  for  ping  EXP3.3  which  preceded  ping 
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EXP3.2  by  ~  104  seconds.  During  this  time,  assuming  they  were  moving  slowly  rela¬ 
tive  to  the  sonar,  one  whould  expect  the  echoes  to  be  at  discernably  different  ranges 

for  the  two  pings.  This,  in  fact  was  not  the  case;  the  features  appeared  at  precisely 
the  same  ranges. 

Alternatively,  the  features  may  be  explained  by  the  presence  of  a  strong, 
well-defined  scattering  layer  at  a  shallow  depth.  Such  a  layer,  consisting  mainly  of 
pre-spawning  populations  of  Pacific  herring,  has  been  consistently  observed  in  the 
spring  months  [Friedl,  1970;  Anderson,  1981;  Frost,  1985],  If,  due  to  refraction  effects, 
the  near  broadside  rays  residing  in  the  null  region  intersect  this  layer  at  appropriate 
ranges,  they  could  possibly  acount  for  the  features.  The  high-resolution  ray  tracing 
plot  created  to  provide  support  for  this  hypothesis,  exposed  additional  evidence  (Fig¬ 
ure  3.16.)  Specifically,  two  caustics  are  clearly  visible  at  ranges  which  correspond  well 
with  the  feature  position.  The  presence  of  the  caustics  in  the  strong  shallow  scattering 

layer  can  probably  account  for  the  intense  scattering  which  produces  the  volume 
features. 

3.5  Vertical  reverberation  correlation  and  ANC 

As  a  final  step,  the  theoretical  model  predictions  presented  in  Section  1.4  on 
the  vertical  correlation  characteristics' of  reverberation  will  be  tested  in  the  context  of 
direct  ANC.  Row  1  will  constitute  the  main  channel.  Given  the  width  of  the  transmit¬ 
ting  beam  pattern,  one  would  expect  the  volume  component  to  be  decorrelated 
between  row  1  and  row  8.  A  twochannel  reference  beam  will  be  used,  consisting  of 
rows  8  and  9.  The  processing  results  for  real  and  simulated  data  are  shown  in  Figure 
3.17.  The  behavior  of  the  REVGEN  data  is  consistent  with  the  model  predictions. 
Cancellation  is  virtually  zero  in  the  early  ranges  as  both  volume  and  boundary  rever¬ 
beration  have  low  vertical  coherence.  The  cancellation  improves  gradually  and  eventu- 
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ally  SRC  is  achieved  as  the  boundary  return  becomes  increasingly  coherent.  Cancella¬ 
tion  improves  with  range  for  the  real  data  as  well,  although  the  trend  is  not  quite  as 
clear.  This  may  be  construed  as  evidence  of  a  degree  of  departure  from  the  spatial 
correlation  structure  predicted  by  the  simple  point-reflector  model. 

3.6  Chapter  summary 

A  number  of  adaptive  filtering  configurations  were  applied  to  the  problems  of 
signal  detection  in  a  reverberation  backround  and  of  selective  reverberation  cancella¬ 
tion.  Adaptive  prewhitening,  ANC  with  fixed  preformed  reference  beams  and  con¬ 
strained  adaptive  beamforming  were  used  to  treat  real  and  simulated  shallow-water 
reverberation  data.  It  was  shown  that,  although  signal  enhancement  is  possible  with 
other  methods,  only  (spatially)  adaptive  beamforming,  implemented  through  a  mul¬ 
tichannel  joint-process  structure,  is  capable  of  achieving  the  SRC  objective.  Certain 
model  theoretical  predictions  on  ANC  performance  limitations  were  shown  to  be  appli¬ 
cable  as  were  the  point-scattering  model  predictions  of  the  spatial  correlation  charac¬ 
teristics  of  reverberation.  REVGEN  simulations  were  performed  to  parallel  the  experi¬ 
mental  data  and  generally  good  correspondence  was  observed  between  the  real  and 
simulated  processing  results. 
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Figure  3.1 


Shallow-water  experimental  geometry.  The  surface,  volume  and  bottom  contribute 
to  the  backscattered  return. 
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Figure  3.2  Raytracc:  ping  EXP3.1.  Velocity  profile  and  corresponding  raytrace  diagram  for 
real  data  (a)  and  REVGEN  simulated  scenario  (b).  Both  surface  and  bottom 
reverberation  is  produced  at  all  ranges. 
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Figure  3.3  Ping EXP 8.1:  Low-pass  time  series  plot  (a)  and  Range-Doppler-Map  (b).  The  two 
synthetic  echos  Tl  and  T2  are  masked  by  reverberation. 
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Figure  3.4  Adapttve  prewhitening  results.  Adaptive  spectral  estimate  and  Range-Doppler-Map  of  prediction  error  output  for 
aLSi  =0.0002  (a,b);  aLSL  =0.002  (c,d)  and  aLSL  =0.02  (e,f).  As  aLSL  increases  the  resolution  of  the  spectral 
estimate  improves  and  the  low-Doppler  echo  vanishes. 
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Figure  3.5 


7l\  jrtrmei  refertnce  btams •  The  R^nge-Doppler-Map  (a)  and  the  tin 
stznlh  Tl  andT^I  ‘w*  LSL  J°mt'process  error  outPut,  show  the  synthetic  ec 

sides  The  d^cJ”  ,°f  e  gFeat  y/  ®nhLanced  «  ^e  reverberation  backround  su 
sides.  The  dB-cancellation  curve  (c)  shows  cancellation  to  improve  with  ran-e 
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Figure  3.6  REVGEN  simulation  of  EXP8,1.  The  LSL  joint-process  error  output  (dotted  line 
in  (a))  remains  much  higher  than  the  ”  desired”  volume  signal,  indicating  that 
fixed,  preformed  reference  beams  are  not  adequate  for  SRC  despite  the  fact  that 
substantial  cancellation  is  achieved  (b). 
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Figure  3.7  Raytracc :  Ping  EXP  3.2.  Velocity  profile  and  corresponding  raytrace  diagram  for 
real  data  (a)  and  REVGEN  simulated  scenario  (b).  Note  that  in  the  real  data 
case  a  shadow  zone  is  formed  in  the  near  surface  and  no  surface  reverberation  is 
produced  at  the  far  ranges.  In  the  simulated  case,  on  the  other  hand,  surface 
reverberation  continues  to  be  produced  at  the  far  ranges. 
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Figure  3.8  EXP  3.2  beam  patterns,  (a)  primary  beam  pattern,  produced  by  adding  all  9 
rows,  (b)  single  vertical  0  *  bearing  cuts  of  the  primary  beam  (solid  line)  and  of  a 
single  constrained  beam  pattern  produced  by  subtracting  two  rows  (dotted  line.) 
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Figure  3*9  EXP  8,2  data .  Real  (LHS)  and  REVGEN  (RHS)  primary  channel:  (a)  standard 
RDM,  (b)  expanded  near-range  plot  (note  the  typical  "hook"  pattern  caused  by 
the  high-angle  boundary  reverberation  returns)  and  (c)  time-series  plot. 
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Figure  3.10  Angle  vs  time  of  arrival .  (a)  Real  (EXP3.2),  and  (b)  REVGEN  (RS3.2)  data. 

These  plots  were  intended  to  illustrate  the  origin  of  the  ”  uncorrelated  noise" 
effect.  In  the  early  ranges,  the  constrained  reference  beam  (dotted  line)  receives 
reverberation  from  angular  directions  where  the  primary  beam  (solid  line)  has  a 
null  response. 


100 


Figure  3.11  ANG  with  single  constrained  reference  element.  Ping  EXP3.2  (LHS)  and 
REVGEN  simulation  RS3.2  (RHS):  (a)  Standard  RDM,  (b)  expanded  RDM.  Sub¬ 
stantial  reverberation  cancellation  is  seen  for  both  real  and  simulated  data  in 
the  far  ranges  but  near-range  performance  is  marginal.  Note  increase  in  noise 
level  in  the  simulated  results,  (c)  Output  time-series  and  (d)  dB-  cancellation 
curve.  A  volume  "signal"  (Vl)  is  exposed  in  the  real  data  output.  Cancellation 
is  seen  to  improve  with  range  for  simulated  data,  level  out  for  real  data. 
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Figure  3.12  SRC  performance  of  single  constrained  reference  beam  (RS8.2).  Composite  (top 
solid  line),  pure  volume  (bottom  solid  line)  and  joint-process  error  output  for 
single  constrained  reference  beam  (dotted  line.)  Clearly,  the  SRC  objective  is 
not  achieved,  (b)  dB-cancellation  curve  and  Signal-to-interference  ratio  in  the 
main  channel.  The  two  curves  are  nearly  mirror  images  of  each  other. 
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Figure  3.13  Adaptive  beamforming  with  8  constrained  reference  elements.  Ping  EXP3  2 
(LHS)  and  REVGEN  simulation  RS3.2  (RHS):  (a)  Standard  RDM,  (b)  expanded 
RDM.  Cancellation  improves  dramatically  for  both  real  and  simulated  data  in 
all  ranges,  (c)  Output  time-series  and  dB-cancellation  curve.  An  additional 
volume  "signal"  (V2)  is  revealed. 
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adaptive  cancellation  beam  (3  reference  elements) 


Figure  3.14  SRC  performance  of  spatially  adaptive  reference  beam  (RS3.2).  (a)  Composite 
(top  solid  line),  pure  volume  (bottom  solid  line)  and  output  of  8-channel  beam- 
former  (dotted  line).  With  the  exception  of  a  near-range  region  (0.5—  l.Oaec), 
SRC  is  achieved,  (b)  Waterfall"  plot  of  the  adaptive  reference  (cancellation) 
beam  produced  by  the  algorithm.  The  adaptive  beam  strives  to  emulate  the 
sidelobe  structure  of  the  primary  beam  in  the  near  ranges. 
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Figure  3.15  SRC  performance  of  spatially  adaptive  reference  beam  (RSS.3): 

case,  (a)  SRC  performance  improves  in  the  near  ranges  and  (b) 
adaptive  cancellation  beam  is  produced. 
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Figure  3.15  (cont.)  Vertical  "cuts”  of  the  adaptive  beam  pattern  at  t  =0.5 sec  (c),  t  =  1.0 sec  (d)  and  t  =2.0 sec  (e)  show 
three  stages  of  the  adaptation  process.  Initially,  the  reference  beam  (dotted  line)  emulates  the 
sidelobe  structure  of  the  primary  beam  (solid  line)  thus  eliminating  the  near-range  "uncorrelated 
noise"  effect.  As  the  high-angle  sidelobe  returns  diminish  and  boundary  reverberation  begins  to  enter 
through  the  main  lobe,  the  reference  beam  progressively  narrows  its  center  null  to  compensate  for  the 
near-grazing  boundary  returns. 


105 


dB  MAGNITUDE 


106 


Figure 


3.16 


RANGE  (m) 

i?ter?reta‘ion-  The  two  volume  "signals"  exposed  in  the  dat£ 
f  2  by  thC  adaptlve  beamforming  operation  (a),  correpond  well  in  rang< 
with  two  caustics  resulting  from  the  sharp  velocity  gradient  in  the  near-surface 
The  raytrace  diagram  (b)  covers  an  angular  range  of  ±3  *  which  is  covered  b) 
the  null  of  the  reference  channel.  Therefore,  the  strong  volume  reverberating 
returns  originating  at  the  caustic  regions  are  protected  from  cancellation  anc 
preserved  at  the  filter  output. 
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Figure  3.17  Vert, cal  correlation  and  ANC.  Real  (EXP3.2)  (a)  and  simulated  (RS3.2)  data  (b) 
The  main  and  reference  channels  are  the  outputs  of  rows  1  and  9  respectively 
Cancellation  improves  with  range  for  the  simulated  data  as  predicted  by  the 
point-scattering  model.  The  real  data  results  are  not  as  conclusive. 


CHAPTER  IV 


SURFACE  REVERBERATION  INTERFERENCE 
IN  UPPER  OCEAN  VELOCIMETRY 

Doppler  sonar  systems  offer  significant  operational  advantages  over  mechani¬ 
cal  flow  meter  devices  and  are  becoming  increasingly  popular  within  the  oceano¬ 
graphic  community  as  remote  velocity  sensing  tools.  Typically,  they  are  narrow-beam 
sonars  which  transmit  short  pulses  of  sound  at  carrier  frequencies  ranging  from 
several  kHz  to  ~  10  MHz.  The  Doppler  shift  of  the  volume  backscatter  has  informa¬ 
tion  about  the  relative  radial  velocity  between  the  platform  and  the  scatterers. 
Assuming  that  the  primarily  biological  scatterers  are  passive  tracers  of  water  motion 
and  that  platform  motion  is  insignificant  on  the  velocity  scales  of  interest,  or  indepen¬ 
dently  measured  and  compensated  for,  this  estimate  may  be  interpreted  as  water 
mass  velocity.  Doppler  sonars  are  most  frequently  categorized  into  pulse-to-pulse 
coherent  and  incoherent  types.  In  a  coherent  system  the  velocity  is  estimated  from  the 
pulse  to  pulse  variations  of  the  echo  phase.  In  an  incoherent  system,  each  return  is 
considered  to  be  an  independent  realization  (sample  function)  of  a  random  process  and 
the  velocity  determination  problen  is  essentially  one  of  spectral  estimation.  A 
classification  of  Doppler  sonar  types  is  given  by  Pinkel  [1980]. 

Sidelobe  interference  from  the  surface,  bottom  or  a  strong  scattering  layer  is 
a  definite  concern  during  a  Doppler  sonar  experiment;  it  can  distort  the  volume  velo¬ 
city  information  or  completely  mask  the  volume  reverberation  component.  This  prob¬ 
lem  has  been  cited  explicitly  by  Pinkel  [1980],  Lhermitte  [1983]  and  Zedel  [1985]  and  is 
very  likely  to  arise  during  a  near-surface  deployment  of  a  horizonally  directed  Doppler 
sonar.  The  near-surface  region  and  the  mixed  layer  in  particular,  is  characterized  by 
dynamic  velocity  features  which  are  related  to  the  air-sea  interaction  process  and  may 
represent  important  momentum  transfer  mechanisms  [e.g.  Pollard,  1970;  Weller  et.al., 
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1985],  For  these  reasons,  it  is  a  region  of  great  interest  to  the  physical  oceanographer. 
Currently,  limited  sampling  of  the  mixed  layer  velocity  field  can  be  obtained  by 
downward-looking  Doppler  sonars  [e.g.  Regier,  1982].  A  horizontal  sonar  would  be  con¬ 
siderably  more  efficient  for  that  purpose  if  not  for  the  surface  reverberation  interfer¬ 
ence  problem  which  is  made  more  acute  by  the  upward  refraction  typically  encoun¬ 
tered  in  the  mixed  layer. 

This  situation  is  a  potential  candidate  for  SRC.  Succesful  application  of  the 
adaptive  techniques  described  in  the  previous  sections  may  significantly  extend  the 
useful  operational  range  of  a  horizontal  Doppler  sonar.  In  this  chapter,  real  and  simu¬ 
lated  data  from  a  pulse-to-pulse  incoherent  sonar  will  be  used  to  examine  the  poten¬ 
tial  of  recovering  the  volume  velocity  information  by  eliminating  the  contamination 
due  to  sidelobe-born  surface  interference. 

4.1  Experiment  description 

The  data  were  collected  in  May  1983  off  the  coast  of  San  Diego,  from 
Research  Platform  FLIP  and  through  the  pulse-to-pulse  incoherent  sonar  system  of 
the  Internal  Wave  Group  (SIO).  A  six-panel  array  was  mounted  at  a  depth  of  38 
meters  on  the  port  side  of  FLIP’S  hull,  its  normal  axis  forming  ~45  0  angle  with  the 
fore-aft  direction.  Each  trapezoidal  panel  consisted  of  210  transducer  elements  which 
were  combined  to  produce  a  narrow  (~  2  * )  beam.  A  maximum  of  3  receive  channels 
was  available  at  75  khz.  Panel  A  was  the  transmitting  array  and  panels  B,C,D  were 
used  to  receive.  The  experimental  geometry  and  the  panel  arrangement  are  depicted 
in  Figure  4.1.  The  beam  pattern  of  a  single  panel,  calculated  by  approximating  the 
hexagonal  transducer  elements  by  circular  pistons,  is  shown  in  Figure  4.2.  It  was  found 
to  be  in  good  agreement  with  the  available  calibration  measurements.  Because  the 
element  spacing  is  somewhat  larger  than  X/2,  spatial  aliasing  occurs  giving  rise  to 
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three  families  of  grating  lobes  which  are  ~15  dB  down  from  the  main  lobe.  The  upper 
two  grating  lobes  (elevation  ^37  #  and  bearing  ~  ±  26  * )  were  expected  to  let  an 
amount  of  surface  reverberation  in  the  receive  channels  during  the  experiment. 

Two  data  sets  were  collected,  one  with  the  array  pointed  horizontally  (run 
12A)  and  the  other  with  the  array  rotated  up  ~  3  *  (run  11A).  A  pulse  length  of  20 
ms  was  used.  Each  transmission  was  followed  by  a  2  sec  "listen  interval"  and  the 
complex  basebanded  signal,  sampled  at  500  Hz,  was  digitally  recorded.  A  step  gain  of 
40  db  was  applied  at  ~  0.3  sec  following  each  transmission  in  order  to  preserve  the 
backscattered  return  at  long  ranges. 

During  data  collection,  the  wind  speed  was  ~25  kn  and  the  sea  surface 
displayed  a  large  swell  with  "chop"  and  "whitecaps.”  The  estimated  wave  height  was 
~  4  m.  The  diagram  in  Figure  4.3  shows  the  orientation  of  FLIP  relative  to  the  wind 
direction  as  recorded  during  data  collection.  It  is  evident  that  the  main  lobe  axis  of 
the  port  horizontal  sonar  nearly  coincides  with  the  direction  of  the  wind. 

4.2  Velocity  estimation  procedures 

In  the  case  of  a  pulse- to- pulse  incoherent  Doppler  sonar,  the  velocity  informa¬ 
tion  is  embedded  in  the  spectrum  of  the  reverberation  return.  Normally,  the  entire 
spectrum  cannot  be  resolved  in  the  spatial  scales  of  interest.  However,  much  useful 
information  resides  in  the  first  and  second  moments  of  the  spectrum  (mean  and  spec¬ 
tral  width,  respectively).  Because  the  autocorrelation  and  the  spectrum  are  Fourier 
transform  pairs,  these  moments  may  be  estimated  in  the  time-domain  based  on  a 
well-known  transform  property.  Namely,  the  spectral  moments  are  equal  to  the 
derivatives  of  the  autocorrelation  function  evaluated  at  zero  lag.  The  derivatives  may 
be  determined  by  evaluating  the  autocorrelation  function  at  one  additional  non-zero 
lag.  This  forms  the  basis  for  the  Spectra!  Moment  Estimation  (SME)  method 
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[Rummler,  1968;  Miller  and  Rochwarger,  1972]  which  is  widely  used  in  oceanic  Doppler 
sonar  applications  [Pinkel,  1980]  and  will  be  employed  here.  Specifically,  the  first  spec¬ 
tral  moment  (mean  frequency)  is  given  by: 


m  = 


i 

2jt  Ar 


tan-1 


Im 

Re  R(t,Ar) 


(4.1) 


where  Ar  is  a  suitably  "small"  time  lag  and  R  is  the  reverberation  autocorrelation 
function.  The  velocity  can  then  be  estimated  through: 
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where  e  is  the  speed  of  sound  in  the  ocean  and  fe  is  the  carrier  frequency. 


(4.2) 


Remarks 

i.  The  SME  method  will  be  effective  when  the  Doppler  spectrum  is  unimodal.  If 
multiple  spectral  peaks  exist,  reflecting  distinct  velocity  "events",  this 
method  will  produce  a  meaningless  "average"  estimate.  It  is  likely  that  high 
resolution  AR  (or  ARMA)  spectral  estimation  techniques  may  offer  significant 
advantages  over  SME  by  resolving  the  multiple  Doppler  features.  For  exam¬ 
ple,  Hansen  [1984]  has  shown  that  velocity  anomalies  in  the  volume  signal, 
presumably  caused  by  fast-swimming  fish,  can  be  detected  through  AR 
modeling.  It  may  be  that  surface  reverberation  can  be  equally  well  detected 
by  such  a  spectral  modeling  scheme  and  the  volume  and  surface  velocity 
components  thus  isolated. 

ii.  Although  several  different  autocorrelation  estimators  exist,  the  work  of  Han¬ 
sen  [1984]  indicates  that  they  result  in  essentially  equivalent  velocity  esti¬ 
mates  in  terms  of  estimate  variance  [Theriault,  1981].  Here,  the  following 
estimate  will  be  used: 
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where  t  is  the  ping  (sample  function)  index. 

4.3  Data  description 

Relative  intensity  profiles  from  a  representative  ~6  min  long  data  segment 
are  displayed  in  Figure  4.4.  Each  profile  is  a  20  ping  average.  Two  surface  reverbera¬ 
tion  features  are  clearly  visible  in  this  plot.  First,  a  near-range  (~  65  m)  peak  is 
observed  which  corresponds  well  in  range  with  the  distance  to  the  surface  in  the  direc¬ 
tion  of  the  upper  grating  lobes.  Evidently,  the  surface  return  is  strong  enough  to 
overwhelm  the  volume  signal  by  ~20  dB  despite  the  reduced  response  of  the  sidelobes. 
As  the  sidelobe  surface  hit  ends,  the  volume  component  reappears  until  a  grazing  main 
lobe  surface  hit,  probably  the  result  of  upward  refraction,  occurs  at  ~400-600  m  and 
dominates  the  return.  The  data  are  significantly  different  for  the  two  runs  in  the  far 
ranges.  It  appears  that  for  the  horizontal  deployment  (runl2a)  FLIP’S  motion  causes 
the  beam  to  periodically  miss  the  surface.  This  seems  to  occur  with  a  ~  1  min  period. 
The  "glitch”  at  ~300  m  is  a  remnant  of  the  step  gain  correction. 

Expanded  intensity  plots  of  the  near-ranges  (prior  to  the  gain  change) 
together  with  the  corresponding  velocity  profiles  calculated  through  SME  are 
displayed  in  Figure  4.5.  Striking  velocity  "anomalies"  are  observed  in  the  vicinity  of 
the  intensity  features.  They  are  mostly  negative  which  is  consistent  with  the 
downwind  orientation  of  FLIP  during  the  experiment  and  display  a  complex,  variable 
character.  The  volume  velocity  estimates  are  of  the  order  of  several  cm  /a  as  expected. 
A  substantial  increase  in  the  variance  of  the  velocity  estimate  is  observed  in  the  200- 
300  m  range  which  can  be  attributed  to  low  SNR. 
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4.3.1  The  Doppler  shift  of  surface  reverberation 

Doppler  shifts  from  the  ocean  surface  have  been  observed  in  experiments 
involving  both  electromagnetic  and  acoustic  wave  transmission.  [Crombie  and  Watts, 
1967;  Pidgeon,  1968;  Valenzuela  and  Laing,  1970,  Igarashi  and  Stearns,  1971;  Swarts, 
1972].  Most  frequently,  the  mean  Doppler  frequency  of  surface  reverberation  (clutter) 
corresponds  to  the  phase  velocity  of  surface  waves  that  satisfy  the  Bragg  resonant 
condition  (i.e.  water  waves  of  wavelength  half  the  transmitted  wavelength  for  grazing 
incidence).  This  ’’diffraction  grating"  effect  was  clearly  demonstrated  in  a  laboratory 
experiment  by  Liebermann  [1963],  By  insonifying  a  randomly  rough  surface,  he  showed 
that  monochromatic  acoustic  radiation  is  preferentially  scattered  from  travelling  sur¬ 
face  waves  with  the  appropriate  wavelength,  resulting  in  Doppler  shifts  characteristic 
of  the  wave  phase  velocity.  Other  possible  Doppler  shift  mechanisms  cited  in  the 
literature  are  orbital  velocities  associated  with  gravity  waves  and  "  whitecaps"  travel¬ 
ling  in  the  direction  of  the  wind  [Hicks  et.  al.,  I960]. 

4.3.2  Surface  "roughness”  vs  bubble  layer:  A  hypothesis  based  on  obser¬ 
vation 

The  single  (180  ping  average)  intensity  and  velocity  estimates  plotted 
together  in  Figure  4.7  reveal  an  interesting  discrepancy:  The  surface  velocity  "ano¬ 
maly"  does  not  coincide  in  range  with  the  surface  reverberation  intensity  peak. 
Specifically,  the  velocity  feature  develops  during  the  second  half  of  the  intensity 
feature.  This  may  indicate  the  presence  of  two  distinctly  different  scattering  mechan¬ 
isms  at  or  near  the  surface.  We  offer  the  following  hypothesis  based  on  the  geometry 
of  the  sidelobe  hit  (Figure  4.6)  and  the  postulated  presence  of  a  layer  of  bubbles 
entrained  by  turbulence  just  below  the  surface.  Such  a  bubble  layer  has  been  con¬ 
sistently  detected  in  moderate  to  high  sea  states  [e.g.  Thorpe,  1981]  and  can  cause 
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intense  scattering.  If  this  layer  undergoes  a  random  (chaotic)  motion  it  can  be 
expected  to  contribute  mainly  to  the  broadening  of  the  Doppler  spectrum  rather  than 
result  in  a  net  Doppler  shift.  At  the  same  time,  a  surface  insonified  "patch"  develops 
which  may  cause  net  Doppler  shifts  associated  with  capillary  wave  velocity  or  whi- 
tecap  drift.  Initially  (aa’  in  Figure  4.6)  the  zero-mean  spectral  broadening  due  to  the 
bubble  layer  return  may  obscure  the  surface  velocity  information.  In  the  latter  stages 
of  the  sidelobe  hit  (bb’  in  Figure  4.6),  on  the  other  hand,  the  volume  component  of  the 
combined  return  diminishes  and  the  surface  return  dominates.  In  that  region,  the  SME 
velocity  estimate  is  in  the  range  of  15  to  30  cm/s. 

The  wavelength  of  the  capillary  waves  which  selectively  scatter  the  75  kHz 
(X  «2em)  incident  acoustic  radiation  will  be,  according  to  the  diffraction  grating  con¬ 
cept: 


2  CO80 


1.252  cm. 


where  0  is  the  grazing  angle  (here  -37  # ).  The  corresponding  capillary  wave  phase 
velocity  is  given,  to  a  linear  approximation,  by: 
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where  T  is  the  surface  tension,  p  is  the  water  density  and  g  is  the  acceleration  of 

gravity.  Taking  T  «=  73. 82dyn /cm  and  p  =1.025  f/em!  [Neuman  and  Pierson,  1966]  we 
have, 


t>£  ~  23.60  cm  /». 

This  value  is  in  good  agreement  with  the  surface  velocity  estimate  and  provides  sup¬ 
port  for  the  "  diffraction  grating"  effect. 

Further  support  for  the  present  hypothesis  is  provided  by  the  velocity  esti¬ 
mate  in  the  ranges  of  the  grazing  main  lobe  surface  hit.  In  such  a  low  grazing  angle 
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scenario,  a  subsurface  layer  of  bubbles  has  often  been  suggested  as  the  dominant 
scattering  mechanism.  This  and  the  fact  that  no  large  velocity  "anomalies"  were 
found  in  that  region  are  consistent  with  our  hypothesis. 

4.4  Sidelobe  interference  cancellation 

An  attempt  was  made  to  utilize  the  limited  beamforming  potential  of  this 
system  to  create  a  reference  beam  suitable  for  surface  reverberation  cancellation.  A 
fixed  constrained  beam  may  be  produced  by  subtracting  two  neighboring  panels  (e.g.  B 
and  C).  Let  this  be  known  as  configuration  Cl.  The  "  uncorrelated  noise"  effect  cited 
in  Section  3.4.1  will  not  be  present  here  because  the  transmit  and  receive  beams  are 
identical.  However,  a  potential  difficulty  lies  in  the  large  separation  distance  between 
the  centers  of  any  two  panels  which  is  d  =  0.526m  or  d  =  26.3X.  Panel-topanel  sub¬ 
traction  will  cause  severe  spatial  "aliasing"  resulting  in  a  finely  structured  "grating 
null"  pattern  which  may  seriously  distort  the  sidelobe  structure  of  the  reference  beam. 
This  may  in  turn  destroy  the  ANC  potential  of  this  particular  configuration. 

4.4.1  Spatial  aliasing  in  constrained  reference  beams 

In  order  to  quantify  the  grating  null  effect,  consider  two  identical  omnidirec¬ 
tional  transducers  spaced  a  distance  d  apart  and  suppose  that  a  signal  x(f)  impinges 
on  the  two-element  array  at  an  angle  0  relative  to  broadside.  The  phase  delay 
between  the  two  elements  will  be 

<t>  =  2?r(<//X)sin0 

and  the  output  of  the  constrained  array  formed  by  subtracting  one  element  from  the 
other  will  be 

xc  =  z(<)  —  z(t)c>*. 

In  addition  to  the  obvious  broadside  null  corresponding  to  0=0*  (representing  the 
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desired  constraint),  additional  "grating"  nulls  will  appear  for 

d/Xsin0«n,  n-1,2,... 

In  our  case,  i/\ « 0.526/0.02  *26.3  and  the  first  grating  null  will  appear  at 
6  «sin-»(l/26.3)  ~2.18  *.  Therefore,  according  to  the  beam  pattern  multiplication 
theorem  [Urick,  1975],  one  may  expect  a  large  number  of  closely  spaced  grating  nulls 
in  the  reference  beam  formed  by  subtracting  the  elements  of  two  panels. 

The  effect  of  panel  separation  on  the  resulting  reference  beam  is  illustrated  in 
Figure  4.8.  The  first  plot  (4.8.a)  is  of  a  vertical  "cut"  of  the  beam  pattern  produced 
by  omnidirectional  elements  separated  by  ~0.526  m,  the  distance  between  two  negh- 
boring  panels.  Note  the  finely  structured  grating  null  pattern.  The  second  plot  (4.8.b) 
displays  a  0  bearing  cut  of  the  resulting  reference  beam  and  shows  the  broadside  null 
in  the  main  lobe,  representing  the  desired  constraint.  The  third  plot  (4.8.c)  in  the 
sequence  is  of  a  -26  cut  and  shows  the  original  grating  lobe  to  be  disected  by  one  of 
the  grating  nulls.  Since  the  preservation  of  this  sidelobe,  the  main  interference  outlet, 
is  crucial  for  succesful  ANC  the  reference  beam  created  by  panel-khpanel  subtraction 
is  a  poor  one  for  our  purpose. 

The  situation  can  be  remedied  by  effectively  bringing  the  two  arrays  closer 
together.  For  instance,  by  allowing  sections  of  the  arrays  to  overlap,  two  offset-phase- 
center  beams  may  be  constructed  which  are  sufficiently  close  to  prevent  the  formation 
of  deleterious  grating  nulls.  The  series  of  plots  in  Figure  4.8.d,e,f  demonstrates  the 

effect  of  reducing  the  distance  between  the  array  phase  centers  to  The  single 

wide  broadside  null  is  seen  to  virtually  eliminate  the  main  lobe  while  leaving  the  origi¬ 
nal  grating  lobe  undisturbed  in  the  reference  beam.  This  configuration  (C2)  is  nearly 
ideal  for  sidelobe  interference  cancellation  and  will  be  tested  through  a  REVGEN 
simulation. 


117 


4.4.2  Processing  results 

Segments  of  the  real  FLIP  data  from  both  experimental  runs  were  processed 
with  the  joint-process  LSL  filter  (j>=3,  0lSL  =0.02).  The  results  (Figure  4.9)  are  con¬ 
sistent  with  the  predictions  of  the  previous  section  regarding  the  experimental  panel 
arrangement  (Cl).  Although  substantial  cancellation  (~10  dB)  is  achieved  in  the  range 
of  the  surface  sidelobe  hit,  a  large  surface  reverberation  component  remains  in  the 
ANC  output.  Lesser  (~  3  dB)  cancellation  of  the  volume  signal  is  observed  owing  to 
the  fact  that  the  main  lobe  is  not  completely  eliminated  from  the  reference  beam  by 
the  narrow  broadside  null.  Cancellation  improves  once  more  (to  ~10  dB)  following 
the  onset  of  the  main  lobe  grazing  surface  hit  but  it  is  unlikely  that  SRC  is  achieved; 
this  is  a  situation  analogous  to  the  one  described  in  Section  3.4.1  ,  where  the  fixed  null 
width  proved  to  be  a  serious  limitation  in  maintaining  a  high  interference-tosignal 
ratio  in  the  reference  channel  and  consequently  achieving  SRC. 

4.4.3  REVGEN  simulation 

Two  simulations  were  performed.  In  the  first  simulation  (RS4.1)  the  experi¬ 
mental  panel  arrangement  (Cl)  was  used.  In  the  second  (RS4.2)  the  "  ideal" 
configuration  (C2)  was  implemented.  In  both  simulations  the  experimental  geometry 
was  slightly  altered  in  order  to  contain  the  computational  burden  imposed  by  the  nar¬ 
row  beam  patterns.  Specifically,  a  very  high  density  of  scatterers  is  required  to 
achieve  near-Gaussian  reverberation  statistics.  By  moving  the  arrays  deeper,  to  240 
m,  and  thus  allowing  the  scattering  volume  to  expand  prior  to  the  initial  surface  hit,  a 
lower  density  is  required  and  substantial  computational  savings  are  gained.  This 
alteration  does  not  constitute  a  serious  departure  from  the  experimental  setting  with 
regard  to  our  objective.  The  volume  and  surface  layers  were  assigned  distinctly 
different  velocities  in  order  to  facilitate  the  performance  evaluation  process.  The  simu- 
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lation  parameters  are  summarized  in  Table  4.1 


Simulation  parameter 

Symbol 

Value 

Carrier  frequency 

ft 

75.0  kHz 

Pulse  length 

T 

20ms 

Sampling  rate 

77 

500  samples/s 

Shell  slice  thickness 

T 

5  m 

No.  of  volume  scatterers  per  slice 

~R!7~ 

3000 

No.  of  surface  scatterers  per  slice 

1000 

Layer  height 

H 

260  m 

Sound  velocity 

c 

1500  m/s 

Volume  reverberation  coefficient 

-70  db 

Volume  velocity 

vf 

(5,1)  cm/s 

Surface  reverberation  coefficient 

-20  dB 

Surface  velocity 

(200,20)  cm/s 

Table4.1  REVGEN  simulation  (RS4.1) 


In  addition,  the  calculated  FLIP  panel  beam  pattern  Avas  inserted  in  the  simulation. 

Intensity  and  velocity  profiles  from  the  simulated  data  are  displayed  in  Fig¬ 
ure  4.10  a,b.  Note  that  the  intensity  and  velocity  surface  "events"  coincide  in  range 
in  the  absence  of  a  simulated  subsurface  zero-mean-velocity  layer.  The  processing 
results  for  Simulation  RS4.1  are  shown  in  Figure  4.10  c,d.  Marginal  intensity  cancella¬ 
tion  is  obtained,  similar  to  the  real  data  results.  Virtually  no  improvement  in  the  velo 
city  estimate  is  seen.  The  results  for  Simulation  RS4.2,  shown  in  Figure  4.10  e,f,  are 
dramatically  different.  No  evidence  of  the  surface  reverberation  intensity  feature 
remains  and  nearly  complete  recovery  of  the  volume  velocity  information  is  achieved. 

Remarks 

i.  The  "ideal"  beam  resulting  from  configuration  C2  is  not  suggested  as  the 

best  possible  practical  implementation  of  an  adequate  reference  beam  for  this 
application.  External  transducer  elements  emulating  the  sidelobe  structure 
of  the  main  beam  should  be  equally  effective. 
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ii.  No  simulations  were  performed  for  the  main  lobe  grazing  surface  hit.  The 

results  of  Section  3.4  apply  equally  well  here.  SRC  can  be  achieved  and  there¬ 
fore  volume  velocity  recovered  only  via  spatially  adaptive  beamforming 
implemented  through  a  multichannel  joint-process  structure. 

4.5  Chapter  summary 

Data  from  a  horizontal  pulse-to-pulse  incoherent  Doppler  sonar  system  were 
used  to  estimate  the  velocity  features  of  surface  reverberation.  A  hypothesis  involving 
the  presence  of  a  thin  layer  of  bubbles  just  below  the  surface  was  found  to  be  con¬ 
sistent  with  the  observed  intensity/velocity  structure  of  the  surface  return.  Spatial 
aliasing  in  constrained  reference  beams  was  found  to  be  detrimental  to  SRC. 
REVGEN  simulations  of  this  experiment  were  used  to  show  that  SRC  and  volume 
velocity  recovery  are  possible  if  the  spatial  aliasing  problem  is  aleviated. 
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Figure  4.1  FLIP  experimental  geometry.  Panel  A  was  used  to  transmit  and  panels  B  C  D  to 
receive.  Surface  reverberation  entering  through  the  sidelobe  contaminates  the 
volume  reverberation  return. 
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FLIP  sonar  simulation 

single  pone]  beom  pattern 


Figure  4.2  Beam  pattern  of  a  single  FLIP  panel.  The  upper  two  grating  lobes  at  -37  *  eleva¬ 
tion  and  ±  26  *  bearing  are  the  surface  reverberation  outlets. 
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Figure  4.3  FLIP  orientation  during  data  collection.  The  main  lobe  axis  nearly  coincides  with 
the  wind  direction. 
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Figure  4.4  Representative  intensity  profiles ;  (a)  run  11A,  (b)  run  i2A.  Each  profile  is  a  20 
ping  average.  The  periodic  intensity  fluctuations  in  the  far  ranges  reflect  the 
motion  of  FLIP. 
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Figure  4,6  Near-range  intensity  and  velocity  profiles',  (a,c)  run  1JA,  (b,d)  run  12A.  A  definite  velocity  "event"  can  be  seen 
associated  with  the  surface  reverberation  intensity  feature. 
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Figure  4.8  The  near-range  intensity/ velocity  "  anomaly” ;  (a)  run  11A,  (b)  run  12A.  The  velo 
city  feature  is  seen  to  develop  during  the  second  half  of  the  intensity  feature. 
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Figure  4.7  Illustration  of  the  "  bubble  layer"  hypothesis.  Towards  the  end  of  the  sidelobe  sur¬ 
face  hit  the  volume  part  of  the  combined  volume/surface  scattering  region 
diminishes,  allowing  the  surface  velocity  component  to  emerge. 


dB  magnitude  dB  magnitude  dB  magnitude 

t 

-60.0  -  30.0  0.0  -  60.0  -  30.0  0.0  -  60.0  -  30.0 


elevation  (deg) 


elevation  (deg) 


Figure  4.8  Effect  of  panel  separation  on  the  reference  beam.  The  equivalent  omnidirectional  beam  pattern  for  configuration 
Cl  (a)  is  characterized  by  a  multitude  of  grating  nulls;  in  addition  to  the  desired  broadside  null  (b),  more  nulls 
are  obtained  in  the  reference  beam,  one  of  which  disects  the  original  grating  lobe  (c).  The  corresponding  plots 
for  configuration  C2  show  the  single  wide  null  (d)  to  completely  eliminate  the  main  lobe  from  the  reference 
channel  (e)  while  leaving  the  grating  lobe  unaffected  (f). 
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Figure  4.9  ANC  results  for  real  FLIP  data.  (a,b)  Run  11A,  (c,d)  run  12A.  Although  &  degree  of  cancellation  is  obtained, 

SRC  is  clearly  not  achieved  with  the  experimental  configuration  (Cl). 
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Figure  4,10  ANC  results  for  simulated  data.  Intensity  and  velocity  profiles  of  simulated  data  (a,b),  data  processed  with 
configuration  Cl  (c,d)  and  data  processed  with  "  ideal"  configuration  C2  (e,f). 
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CHAPTER  V 


SEA  BEAM  SIDELOBE  INTERFERENCE  CANCELLATION 

Sea  Beam  is  a  system  developed  by  the  General  Instruments  Corporation 
which  uses  time  of  arrival  information  obtained  through  a  multibeam  echo-sounder  to 
produce  high  resolution  contoured  charts  of  the  ocean  floor.  It  transmits  7  ms  pulses 
at  a  carrier  frequency  of  12.158  kHz  with  adjustable  repetition  rate.  It  employs  a  20- 
element  transmit  array  mounted  along  the  ship’s  keel  and  a  40-hydrophone  receive 
array  which  lies  athwartships.  The  resulting  transmit  beam  pattern  is  wide  (54  0 ) 
athwartships  and  narrow  (2  2/3  0  )  in  the  fore-aft  direction.  The  opposite  is  true  for 
the  receive  beam  (2  2/3°  and  20°,  respectively).  A  total  of  16  receive  beams  is 
formed  by  electronically  steering  at  athwartships  intervals  of  2  2/3*.  An  illustration 
of  the  overall  system  geometry  is  given  in  Figure  5.1a.  The  calculated  transmit  and 
receive  beam  patterns  are  displayed  in  Figures  5.1b  and  5.1c,  respectively. 

Sea  Beam  is  now  widely  in  use  within  the  oceanographic  community  and  has 
contributed  greatly  in  advancing  the  study  of  sea  floor  morphology  by  making  large- 
scale,  detailed  surveys  possible.  It  is  generally  considered  to  be  a  highly  useful  and 
reliable  scientific  tool.  However,  years  of  operational  experience  with  the  system 
installed  on  the  R/V  Thomas  Washington  of  the  Scripps  Institution  of  Oceanography 
(SIO)  have  revealed  a  number  of  artifacts  in  the  Sea  Beam  output  which  may  be 
misinterpreted  as  seafloor  features.  Recently,  de  Moustier  [1985a],  using  a  data 
acquisition  system  developed  by  the  Marine  Physical  Laboratory  (MPL)  of  SIO,  col¬ 
lected  the  acoustic  envelope  from  all  16  receive  channels.  Through  these  data  he  was 
able  to  identify  the  underlying  causes  of  artifacts  such  as  the  "  omega”  and  "  tunnel" 
effects.  He  concluded  that  the  bathymetric  artifacts  result  from  echo  detection  and 
processing  errors  associated  with  certain  characteristics  of  the  sea  floor  such  as  a  sud¬ 
den  change  in  slope,  or  caused  by  interference  entering  through  the  sidelobes  of  the 


130 


131 


receive  beams. 

5.1  Sea  Beam  sidelobe  interference 

Sidelobe  interference  becomes  a  serious  problem  when  the  sea  floor  surveyed 
is  relatively  flat  and  highly  reflective.  In  that  case,  the  strong  near-specular  return 
entering  the  main  lobe  of  the  downlooking  beam  "leaks"  through  the  sidelobes  of  all 
other  beams  (Figure  5.2a)  and  is  present  in  all  receive  channels  at  levels  comparable 
to  the  main  lobe  backscattered  returns.  The  sidelobe  interference  "ridge,"  character¬ 
ized  by  its  simultaneous  arrival  in  all  receive  channels,  can  be  clearly  observed  in  the 
typical  Sea  Beam  acoustic  record  shown  in  Figure  5.2b.  In  certain  operational  states 
the  system  will  track  the  sidelobe  response  and  interpret  it  as  a  synchronous  arrival 
from  all  beams.  The  resulting  contoured  output  is  a  concave-up  half  cylinder  (Figure 
5.3).  This  is  known  as  the  "tunnel  effect." 

Perhaps  a  more  serious  consequence  of  the  sidelobe  interference  is  that  it 
masks  the  backscattered  acoustic  return  in  the  interior  (near-vertical)  beams.  As  a 
result,  it  seriously  inhibits  the  process  of  extracting  additional  information,  beyond 
bathymetry,  from  the  Sea  Beam  acoustic  signals.  For  example,  the  grazing  angle 
dependence  of  .the  backscattered  intensity  may  be  an  important  indicator  of  the  sea 
floor  roughness  characteristics.  In  attempting  to  estimate  this  function  from  Sea 
Beam  data,  de  Moustier  [1985b]  noted  that  the  interior  beams  were  rendered  useless 
by  sidelobe  interference.  Unless  the  interference  is  removed,  any  potentially  impor¬ 
tant  information  possesed  by  the  interior  beams  will  remain  unexplored. 

Here,  we  offer  a  possible  remedy  for  the  sidelobe  interference  problem 
through  a  simple  application  of  Adaptive  Noise  Cancelling. 
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5.2  Sidelobe  interference  cancellation 

The  same  mechanism  responsible  for  creating  the  sidelobe  interference  prob¬ 
lem  offers  the  potential  for  its  removal  through  ANC.  Specifically,  the  downlooking 
beam  is  a  natural  reference  channel  which  could  be  used  to  operate  sequentially  on 
the  other  15  beams.  Because  the  interference-to-signal  ratio  is  high  in  both  the  main 
and  reference  channels,  a  substantial  improvement  in  signal-to-interference  ratio  (SIR) 
is  predicted  by  ANC  theory  (Equations  2.5a,b). 

Application  of  this  concept  on  real  Sea  Beam  data  is  not  currently  possible 
because  only  the  envelope  of  the  acoustic  return  is  recorded.  A  modification  of  the 
MPL  data  acquisition  system  is  under  way  to  allow  the  collection  of  amplitude  and 
phase  information  through  quadrature  sampling  of  all  16  beams. 

5.2.1  REVGEN  simulation 

Preliminary  conclusions  on  the  feasibility  of  this  technique  may  be  reached 
through  a  REVGEN  simulation.  A  simulation  run  (RS5.1)  was  performed  which  retains 

all  known  Sea  Beam  system  parameters  including  the  calculated  transmit  and  receive 
beam  patterns.  A  summary  of  this  REVGEN  run  (RS5.1)  is  given  in  Table  5.1. 


Simulation  parameter 

Symbol 

Value 

Carrier  frequency 

Tc 

12.168  kHz 

Pulse  length 

T 

7ms 

Sampling  rate 

/. 

1000  samples/s 

Shell  slice  thickness 

T 

5  m 

No.  of  volume  scatterers  per  slice 

300 

No.  of  bottom  scatterers  per  slice 

K 

300 

Layer  height 

H 

4000  m 

Sound  velocity 

c 

1600  m/s 

Volume  reverberation  coefficient 

-$0  db 

Bottom  reverberation  coefficient 

*b 

-50  db 

Table  5.1  REVGEN  simulation  RS5.1 
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A  total  of  4  receive  beams  (steered  at  0‘,  3*,  6*  and  9*)  was  deemed  sufficient  for 
the  purpose  of  this  simulation.  The  synthetic  Sea  Beam  record  is  displayed  in  Figure 
5-4. a.  The  sidelobe  interference  pattern  is  clearly  visible  and  bears  close  resemblance 
to  the  corresponding  real  data  feature. 

The  3  outer  beams  (3  ’ ,  6 ' ,  9  * )  were  processed  through  the  joint-process 
least-squares  filter  (p=3,  aLSL  =0.02)  with  the  0‘  beam  serving  as  the  reference  chan¬ 
nel.  The  results  (Figure  5.4.b)  are  promising.  The  sidelobe  interference  has  been 
succesfully  removed  with  no  evidence  of  distortion  in  the  remaining  signal. 

Remarks 

i.  The  emergence  of  the  sidelobe  interference  "spike"  from  a  backround  of 
volume  reverberation  represents  an  intensity  "jump"  of  ~20-40  dB.  Accord¬ 
ing  to  Simulation  3  (scale  change)  of  Chapter  2,  the  stochastic  approximation 
algorithms  (i.e.  LMS  and  GRL)  would  be  unsuccessful  in  this  setting.  During 
exploratory  cancellation  attempts  with  the  LMS  and  GRL  on  the  simulated 
Sea  Beam  data,  the  two  filters  were  driven  to  instability  as  predicted  by 
Simulation  3. 

ii.  If  similar  cancellation  performance  is  obtained  for  real  Sea  Beam  data,  the 
"tunnel  effect"  can  be  eliminated.  However,  additional  simulations  are 
needed  to  determine  the  extent  to  which  recovery  of  the  interior  beam  back- 
scattered  return  can  be  achieved. 

iii.  An  alternative  operation  which  may  be  effective  in  removing  the  Sea  Beam 
sidelobe  interference  involves  the  placement  of  appropriate  spatial  nulls  in  all 
receive  beams.  Essentially,  an  entirely  new  beamforming  scheme  must  be 
implemented;  it  would  require  tapping  the  outputs  of  all  40  receiver  elements 
and  hence  further  modification  of  the  current  Sea  Beam  data  acquisition  sys- 
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tern. 

5.3  Chapter  summary 

It  was  demonstrated  through  a  REVGEN  simulation  that  the  same  mechan¬ 
ism  responsible  for  creating  the  Sea  Beam  sidelobe  interference  problem  offers  the 
potential  for  its  removal  through  ANC,  pending  a  minor  modification  in  the  Sea  Beam 
data  acquisition  system. 


Figure  5.1  Sea  Beam  experimental  geometry  and  beam  patterns .  (a)  A  schematic  diagram  showing  the  relationship  between 
the  transmit  beam  and  the  sixteen  receive  beams.  The  calculated  transmit  and  receive  beam  patterns  are 
displayed  in  (b)  and  (c)  respectively. 


d-B  MAGNITUDE 
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Figure  5.2  Sea  Beam  txdelobe  interference,  (a)  Generating  mechanism:  The  near-specular 
return  entering  the  main  lobe  of  the  downlooking  beam  (solid  line),  "leaks" 
through  the  sidelobe  structure  of  an  outer  beam  (dotted  line),  (b)  A  typical  Sea 
Beam  record  displaying  the  characteristic  sidelobe  interference  "  ridge"  (from  de 
Moustier,  1985a). 


Figure  6.3  The  "  tunnel  effect .  As  a  result  of  sidelobe  interference,  the  contoured  output  of  Sea  Beam  resembles  a  trough 
approximately  centered  on  the  ship’s  track  (from  de  Moustier  1985a). 
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SEA  BEAM  Simulation 


Figure  5.4  Sea  Beam  sidelobe  interference  cancellation,  (a)  REVGEN  simulated  data  and  (b) 
ANC  processing  results.  By  using  the  downlooking  (<T)  beam  as  the  reference 
channel,  the  sidelobe  interference  ridge  is  successfully  removed  from  the  outer 
beams  of  this  simulated  Sea  Beam  record. 


CONCLUSIONS 


The  central  theme  of  this  thesis  has  been  the  application  of  the  adaptive 
beamforming  concept  to  the  problem  of  selective  reverberation  cancellation.  We  have 
exposed  relevant  reverberation  properties,  evaluated  a  representative  cross-section  of 
adaptive  algorithms  and  demonstrated  experimentally  the  feasibility  of  selective 
reverberation  cancellation  through  real  and  simulated  data  from  three  active  sonar 
experiments. 

The  spatial  correlation  properties  of  reverberation  predicted  by  the  point¬ 
scattering  model  were  shown  to  be  directly  applicable  in  a  reverberation  cancellation 
context.  Specifically,  the  predicted  disparities  in  vertical  correlation  between  the 
boundary  and  volume  returns  provide  a  natural  mechanism  for  the  separation  of  the 
two  components  through  adaptive  noise  cancelling. 

Based  on  theoretical  considerations  substantiated  by  computer  simulations,  it 
was  shown  that  lattice  joint-process  filters  have  superior  convergence  properties  over 
direct  tapped-delay-line  implementations  in  the  presence  of  power  disparities  in  the 
reference  channel.  In  addition,  deterministic  least-squares  solutions  were  shown  to  per¬ 
form  better  than  stochastic  approximation  (gradient)  solutions  during  drastic  power 
transients  in  the  reference  channel. 

Real  and  REVGEN-simulated  composite  reverberation  data  from  a  shallow- 
water  experiment  were  used  to  determine  the  degree  of  cancellation  obtainable 
through  adaptive  processing.  It  was  shown  that  significant  SNR  enhancement  can  be 
achieved,  for  echoes  embedded  in  a  reverberation  backround,  through  ANC  with  fixed, 
preformed  reference  beams.  Morover,  it  was  shown  that  the  weak  volume  reverbera¬ 
tion  component  can  be  extracted  from  a  composite  process  through  boundary  rever¬ 
beration  cancellation  obtained  via  constrained  adaptive  beamforming. 
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Data  from  a  near-surface  deployment  of  a  Doppler  sonar  revealed  an 
interesting  "anomaly"  in  tbe  intensity/velocity  structure  of  surface  reverberation 
which  led  to  a  hypothesis  involving  the  presence  of  a  thin  layer  of  bubbles  immediately 
below  the  surface.  In  addition,  it  was  demonstrated  through  a  REVGEN  simulation 
that  recovery  of  the  volume  velocity  signal  is  possible  through  ANC  if  appropriate 
reference  beams  are  constructed. 

A  REVGEN  simulation  of  the  Sea  Beam  bathymetric  system  was  used  to 
demonstrate  that  the  sidelobe  interference  problem  giving  rise  to  the  "tunnel  effect" 
artifact  can  be  aleviated  through  adaptive  processing. 

There  are  several  directions  in  which  the  results  developed  here  can  be 
extended.  Experimental  measurements  of  the  spatial  correlation  characteristics  of  oce¬ 
anic  reverberation  under  a  variety  of  environmental  conditions,  would  be  very  useful 
in  assessing  the  reverberation  cancellation  potential  of  a  given  experimental  arrange¬ 
ment.  In  addition,  it  would  be  of  interest  to  pursue  further  the  promising  simulation 
results  of  Chapter  III.  If  an  appropriate  reference  beam  can  be  constructed,  it  may 
indeed  be  possible  to  recover  volume  velocity  information  by  selective  removal  of  the 
surface  reverberation  component.  Morover,  significant  potential  benefits  may  be 
gamed  from  the  application  of  the  sidelobe  interference  cancellation  scheme  suggested 
in  Chapter  I\  to  real  Sea  Beam  data,  pending  a  minor  modification  in  the  Sea  Beam 
data  acquisition  system  to  allow  the  collection  of  both  amplitude  and  phase  informa¬ 
tion.  Finally,  recently  developed  algorithms  based  on  the  more  general  ARMA  model 
(Morf  and  Lee,  1978;  Lee,  1980]  offer  the  potential  for  better  representation  of  the 
reverberation  process  and  may  lead  to  improved  cancellation  performance. 


APPENDIX 


A.  Levinson’s  algorithm 


This  algorithm  is  an  efficient  solution  of  the  linear  prediction  problem  (Equa¬ 
tion  2.4),  when  the  autocorrelation  matrix  R&  is  Toeplitz.  We  begin  with, 


[An.  A <», . Rf,  =  [£p',0,....,0] 


(A.l) 


The  unknowns  are  the  {AW}  and  E‘.  The  aim  is  to  determine  E$M  and  { }  in  a 
way  that  takes  maximum  advantage  of  the  previous  computations  made  to  find  Ep 
and  { A^} . 

The  method  starts  by  trying  an  "obvious"  solution,  assuming  that  adding  a 
zero  to  the  previous  solution  may  work.  It  would  work  if  in  the  resulting  equation 


(A.2) 


the  term 


is  zero.  Since  this  will  not  happen  in  general,  we  have  to  find  a  way  of  forcing  ap  to 
zero.  For  this  we  introduce  the  auxiliary  (reversed)  equations 


(A.3) 


Next,  we  form  a  weighted  combination  of  A.2  and  A.3 


from  which  it  is  easy  to  see  that  choosing 


k?  =  ^p[E;r 

gives  a  solution  of  the  extended  equation,  i.e. 
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Similarly,  with  the  help  of 


we  can  update  the  { B^}  as 


v— wr 


. -nuu-i  -  w . si'u,i +k>\im 

The  recursions  for  E'  and  Ep  are 

Rt>«  =e;-op\e;]^p 
e;m  =e;  -pp[E;)-lap 

This  concludes  the  multichannel  prediction-error  Levinson  algorithm. 


B  Joint  multichannel  complex  least  squares  lattice  (JMCLSL) 

The  multichannel  least  squares  lattice  algorithm  is  summarized  as  follows: 

Initialization  (i  =0,  1,...  ,  p) 


r,(-l)  =  0  ,  .Vp 


Etf-l)  =  «Im  ,  «Vp 


A.(-i)  =  o  ,  «yo 

l)  =  0  ,  «Vp 

Ktf.-l)  =  0 


(Ba) 

(Bb) 

(Be) 

(Bd) 

(Be) 


Time  update  (n  >  0) 
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e#(«)  -  r#(n)  -  x(n)  (Bf) 

E^(n)  =  E'(n)  =  (l  —  Qjclsl )  EJ(n  — l)  +  x(n)xw(n)  (Bg) 

7-.(")  =  0  (Bh) 

e4(n)  =  d(n)  (Bj) 

Order  update  (i  =0,1,. ..,p) 

Lattice 

^•(n)  =  (1  —  Qjclsl)  A,(n— 1)  (Bj) 

e,_,(n)r,^(n-l) 

1  -  7,-c(n-l)  ’  ** 

K,?(n)  =  Ali(n)E^i{n)  ,  tVO  (Bfc) 

K^n)  =  A,-(n)E{*^(n— 1)  ,  «yo  (BJ) 

e.(n)  *  «.*— i(n)  +  K/(n)r,_,(n— 1)  ,  «yo  (Bm) 

r,(n)  =  r,_,(n-l)  +  K*(n)e,_1(n)  ,  yo  (Bn) 

E/(n)  =  E*_,(n)  -  A,(n)Et^(n-l)A/,(n)  ,  yo  (Bo) 

E/(«)  -  EL(n-l)  -  Afln)Ea(n)A,(«)  ,  ,yo  (Bp) 

=  7,-c(n)  +  r,^1)(n)E1r_|(n)r,-_1(n)  ,  yo  (Bq) 

M«)  =  (1  -  aJCLSL)Af{n-l)  -  (gr) 
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Ktf»)  =  A/t»)BT(«)  (Bs) 

e.^n)  =  e/_,(n)  +  K/(n)r,(n)  (gt) 

Predictors 

A,f’)(n)  -  K/(n)  (Bu) 

B(0(n)  =  K'(n)  (Bv) 

Ai’)(n)  =  Ai’^n)  +  Kfln)  B^l(n-l)  (gw) 

l<ifc<»— 1 

Bi°(»)  =  B^)(»-l)  +  K/(n)  Ai‘‘-*)(n)  (gx) 
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